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Dedication of the New NBS Laboratories 


On November 15, 1966, Secretary of Commerce 
John T. Connor dedicated the new ultra-modern facili- 
ties of the National Bureau of Standards at Gaithers- 
burg, Md. This $120 million laboratory complex serves 
as the nerve center of the Nation’s measurement and 
standardization activities. With three Institutes de- 
voted to Basic Standards, Materials Research, and 
Applied Technology, the Bureau’s new home has been 
designed to meet the needs of today and the challenges 
of the future in the physical sciences. 

Dedication ceremonies were held in a courtyard 
between the Bureau’s new library and one of the 
laboratory buildings. In the center of the courtyard 
is an apple tree descended from the tree under which 
Sir Isaac Newton was sitting when, according to the 
story, an apple fell on his head, leading him to conceive 
the law of universal gravitation. 

Surrounding this tree was an audience of 2500, 
including distinguished guests from government, indus- 
try, and the world scientific community. The Reverend 
Edward G. Latch, Chaplain of the House of Representa- 
tives, gave the invocation and the benediction. _Intro- 
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Eleven stories high, the Administration Building at the National Bureau of Standards’ new site near Gaithersburg, Md., 


above the other buildings in the laboratory complex. 
not requiring laboratory space. 


duced by Dr. A. VY. Astin, Director of the National 
Bureau of Standards, Secretary Connor said: 

“It is no accident of history that the establishment of 
NBS in 1901 coincided with the beginning of America’s 
transformation into the industrial giant of the twentieth 
century. . Today the demands made upon the 
Bureau and the nature of its program reflect the com- 
plexity of our atomic and space-age world.” 

The Secretary began his remarks with a message 
from President Johnson. The President noted that 
throughout its 65 years’ existence the Bureau had 
advanced the frontiers of measurement in pace with 
the demands of science and industry; he extended his 
congratulations and best wishes for continuing achieve- 
ment. 

Secretary Connor stated that at age 65 the Bureau 
had significant accomplishments behind it and a future 
promising even greater contributions to the Nation. 
He noted that the new facilities are a symbol of the 
Bureau’s future, but that the germ of future achieve- 
ment lies in the members of the Bureau’s staff. 

In connection with the dedication, the Secretary 
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The building houses the Director and his staff as well as other NBS activities 








The new NBS site is 20 miles northwest of the old Washington 
site, at the Darnestown exit of Interstate Highway 70S. It borders 
and is west of 70S and Gaithersburg, Md. The plan view shows 
the layout of existing buildings. 


sponsored a Symposium on Technology and World Trade 
which was held at the Bureau during the two days 
following the dedication. It featured as speakers 
experts from industry) government, and education of 
many countries, attracting an attendance of nearly 
500. 

Following the dedication ceremonies, the invited 
guests were given a tour of the installation. On 
Saturday, November 19, the general public was invited 
to visit the new facilities and to view over 100 labora- 
tories and demonstrations. In an eight-hour period 
over 20,000 people roamed the site. Some of the 
things they saw were a 100-ft linear accelerator which 
produces one of the world’s most intense electron 
beams; a recently completed high-flux research reactor; 
a 12-million pound capacity hydraulic testing machine; 
the NBS Museum containing historical instruments, 
equipment, standards, and materials pertaining to the 
Bureau's activities during its 65-year history; a vault 
containing the national measurement standards — or 
parts symbolizing those standards; a demonstration of 
glass blowing techniques used to make highly unusual, 
and precise apparatus; and many of. the Bureau's 
laboratories in the fields of metrology, electricity, 
mechanics, heat, atomic physics, radiation physics, 
analytical chemistry, physical chemistry, polymers, 


metallurgy, and building research. 

By the close of 1966, transfer of the Bureau’s Wash- 
ington (D.C.) laboratories to the Gaithersburg site will 
be largely completed and most of the 2700 staff mem- 
bers will be relocated in 15 major buildings at the site. 
Of about 300 persons remaining at the Washington site, 


most will be relocated in new special-purpose labora- 
tories now under construction at Gaithersburg. In 
addition, the Bureau will continue to maintain a staff 
of about 600 in Boulder, Colo., for work in radio stand- 
ards and cryogenics, and a staff of about 300 at the 
Clearinghouse for Federal Scientific and Technical 
Information at Springfield, Va., as well as small groups 
of employees at widely scattered field stations. 

Although the Bureau is located 20 miles northwest of 
its old site, the mailing address remains: National 
Bureau of Standards, Washington, D.C. 20234. The 
new telephone number is (Area Code 301) 921-1000. 
To reach the Bureau by car, use the Darnestown exit 
on Interstate Highway 70S. 
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Symmetrical Bending of Thin Circular Elastic 
Plates on Equally Spaced Point Supports 


A. F. Kirstein and R. M. Woolley 
Institute for Basic Standards, National Bureau of Standards, Washington, D.C. 20234 


(August 23, 1966) 


\ special application of Bassali’s solution for transverse flexure of thin elastic plates supported 
at several points is presented for the case of symmetrical bending. Equations for moments, shearing 
forces, and stresses are developed which may be useful for design purposes. The experimental 
results although limited in quantity are in good agreement with the theoretical predictions. 
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1. Introduction 


The determination of bending moments, twisting moments, and shearing forces in a thin 
circular elastic plate subjected to symmetrical bending is a problem which is often encountered 
in the design and analysis of structural elements or systems. This study deals,with the solution 
of this problem for a thin circular elastic plate supported at points equally spaced along a concen- 
tric support circle and subjected to a transverse load which is symmetrically distributed over a 
concentric circular area. These structures may be typified as end closures, bulkheads, and dia- 
phragms. Usually the analysis of such a structure is simplified by the introduction of engineering 
approximations which pertain to the particular structural system under examination. However, 
this paper presents a special application of a more general solution developed by Bassali [1] ' 
which more closely represents the conditions realized in practical structures and obviates the 
necessity for some of the approximations. 

This specialized treatment of Bassali’s theory provides the equations necessary to calculate 
moments, shearing, forces, and associated stresses anywhere within the plate. It also shows 
that the expressions for maximum bending stresses at the center are independent of the angular 
orientation and the number of supports. Further reduction of these expressions ‘result in the 
Grashof [2] equations. 

A comparison between a limited amount of experimental results and the theoretical pre- 
dictions of tangential strains along the concentric support circle show good agreement. This 
good agreement for strains along with that for deflections [3] tend to substantiate the theory. 


2. List of Symbols 


-, @ polar coordinates 

c radius of the plate 

a radius of the concentric support circle 

b radius of the loaded area (region 1) of the plate 
h thickness of the plate 


' Figures in brackets indicate the literature reierences at the end of this paper 


] 





=rlc 
=alc n'=n+2 
= b/c qg =b'lc ; 
number of support points (m = 3) 
=2z/m, polar angle subtended by support points 
s =sa, polar angle subtended by the sth support point (s=1, 2, 3, . 
s= O— 8, 
modulus of elasticity 
Poisson’s ratio 
= (3+ v)/(v—1) 
intensity of load 
transverse load intensity over 0 


r <b (region 1) 


=rs 
transverse load intensity over b <r <c (region 2) 


total load on the plate defined by eg (1.02) 
moment 

shearing force 

stress 


strain. 
Subscripts: 


1, 2 refer to regions | and 2, respectively 
r, t refer to radial and tangential, respectively 
rt refers to twisting. 


3. Background 


Bassali obtained the solution for the problem of flexure of a thin circular elastic plate sup- 
ported at several arbitrarily located interior points and transversely loaded over a circular area 
eccentrically located with respect to the center of the plate. Bassali considers the intensity of 


the transverse load over the circular area indicated as region | in figure 1 to be 
pi=por"™-2, (n > 2) (1.01) 


and the intensity over region 2 to be po=0. This, of course, shows that the load is distributed 
symmetrically with respect to the center of region 1, and the value of n defines the load distribution 
(n=2 represents uniform load distribution). The total load on the plate is given by 


bn 
Po = 27po — (1.02) 
n 


The boundary of the plate is considered to be free as the plate is supported at interior points, but 
no special treatment is required when some or all of the supports lie on the boundary. 

In a previous paper [3] a special application of this more general solution was used to obtain 
a method of determining the deflection of a plate subjected to symmetrical bending with the load 
uniformly distributed over region 1. A comparison of theoretical and experimental results indi- 
cated that the theory adequately predicted the deflection of the plates over the range of ‘geometries 
tested. Further observations indicated that the theory accounted for a constraining effect on the 
deflection which for the most part was due to the annular portion of the plate overhanging the 
support circle. It was also noted that this constraining effect decreased as the number of supports 
were increased, but did not vanish when the number of supports became sufficiently large to 
produce a support condition equivalent to that of a simple continuous line support. 
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FIGURE 1. Concentric arrangement of the circular plate. 


With regard to the preceding discussion it was deemed desirable to extend the previous 
work [3] to encompass moments, shearing forces, and stresses, for the purpose of presenting 
equations which may be adapted to the design and analysis of structures of this type. 


4. Moments and Shearing Forces for Symmetrical Bending 


It can be shown that Bassali’s solution, reduced to the case of symmetrica! bending, yields 


the following expressions for the radial and tangential bending moments in region |] 
Pol+v)[.,,] , Pol—v) 
i a eee 4 
: 8armk | loamxKp? 


Wj Se ka = ja, (1.04) 
8amk | l6amkp? : 


respectively, where 





(1 — t?)(1 — p?t?) 
(1—2pt cos gs+p?t?) 


m 
r=|S {In (1—2pt cos gs+ p*t?)—« In (p?—2pt cos gs + t?)— 
l 


s 


2mt* 2mkq” 2mk (p" 7 
+ + m= ——*— + 2mk In q+ (2-1) 1.05) 
(+p) m on mk in q ~ q" ( > 





. -) G-fey7 ~aru—p, 
P iain a —? “OS et 22 — 
Q; b {(x 1) In (1 — 2pt cos gs + p*t?) + (—2pt cos @, +p? 


s=1 





— p*)(1 — #)(1 — pt?) bir I 2__ 2) 4mxKp"™ 
(1 — p*y See a —_ af + mix — 1ip?— 2) * (1.06) 


(= 2pt cos ¢s+ p?t?)?  (p?—2pt cos gs +t? n'q” 





For region. 2, 


loam«p? 


M;,=- 


Pol + v) vs| : Pol —v) | 1.07) 
. Samk 


Po TV) ; x o(1 - ' 
on guiomenees A : Lt. (1.08) 

Samk em 
respectively, where 


(l—e)1—p*t*) _| 


m 
ln (1 —2pt cos ¢,+ p*t?)—k In (p* — 2pt cos gs 4 
§ Bo 2pt COS Ys T p* t)) 


=|) 


s=1 


9 2 
2mK¢ 

+ m———— + 2m In p (1.09) 
K+ | 


and 


‘P—2t(1—p*? (1—p*)(1—t*)(1 — p*t*y 


—2pt cos gst p*t?) (1—2pt cos gt p*t?? 





K(p* — t*)? ) . . 9 eo as 
oe er —~¢ + m(k — 1)(p° — t°)— 2mxp? + 2mkq* 
(p°— Zpt COS Qs 


It is noted here that the equation for the twisting moment is the same for regions | and 2 and 


may be expressed as 


PoAl—v)f 4 {pt sin gdp? —t*) pt sin gs(1 — p*) (1 — t*) (1 — pt”) 
“n= Ss _ ‘ 9,9\9 ~~ 
(1—2pt cos ¢, + p*t?)? 


8armxp? | L =, 1 —2pt cos gs + p?t?) 





Kpt sin @5(p? — t*) : ot sin ~, ) Z 
r re —(x*— 1) arctan mtv «| (1.11) 


(p? — 2pt cos gs + t*) 1 — pt cos ¢sJ 





The 


and that the equation for the tangential shearing force is also the same for regions 1 and 


tangential shearing force is given by 


_ 2pt sin g(1 — t?) (1 — p?t?)) |: 


; lye 2pt cos gs+ t?) (1—2pt cos gst pt?)  (1—2pt cos ¢g,+p7t?? J 
(1.12) 


P, m | 2Kpt sin gs _ 2pt sin ¢s 








O=- 


darrmk | 


2 Bassali’s eq (2.44) for M,,:, was incorrect in that a factor, his P,, was omitted from the summation. No doubt this was a typographical error as his equation for 
V,, was correct, and for the genera! solution of eccentric loading M,,, must equal M,,, at the boundary between regions | and 2 
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The radial shearing forces in regions 1 and 2 as reduced from the general solution are 


G = Po z. k(p* — t?) >. (1 — 2p? + p??)t? 
" 4armk (p? — 2pt cos gs + p?t?)? (1—2pt cos gs + pt?) 








(l—&(1—p*ey | ame 
a + af ra . « 
(1—2pt cos gs + p?t??? m(x — 1) q” (1.13) 


and 


Po j = 2— 4? 1 — 292+ p2t?)t? 92 _ n242)2 
TMK i ba ae ied is. |—mi+1)}, 


— 9 
larmk | = 


Q!.,= 


« 





‘ \(p? —2pt cos gs 4 a) i— 2pt cos gs + pt?) “(1- 2pt cos ¢;+ p?t?? 
(1.14) 


respectively. 


5. The Uniformly Distributed Concentric Load 


The case of the uniformly distributed concentric load with its limiting cases of the concen- 
trated central load and the load uniformly distributed over the entire plate is of more general 
interest than the case of the symmetrically distributed load given in the preceding section. There- 
fore this investigation will deal primarily with the case of the uniformly distributed load. It should 
be noted that no special treatment is required to examine the limiting cases mentioned above, as 
they are implicit in the solution. The concentrated central load is represented by permitting q to 
shrink to zero, and the load uniformly distributed over the entire plate is realized by setting q= 1. 
It remains then only to select the appropriate equations for the region under analysis. 

The equations for moments and shearing forces given in the preceding section were derived 
for the symmetrically distributed concentric load as defined by eq (1.01). In order to obtain similar 
equations for the uniformly distributed concentric load it is a simple matter of supplying the ap- 
propriate values for n, n', and q’, viz, n=2, n'=4, and q'=q/V2. Thus, eqs (1.03) and (1.04) 
become 


/ +p 9 —vV 
u,,=— Pol +) [V,]+ Pol») 
Samk 


[Qi] 


l6am«p” 
and 


9 TV 4 —p 
Mi, = — Pol +v) ry,)—- (1 —v) 


8irmk l6amxp” 


[Qi], 


respectively, where eqs (1.05) and (1.06) become 


m 
=, 

Y, = ; 
7 
s=1 





2,2 2 (-(i—gr) 
In (1—2pt cos gs + p?t?)—k In (p? — 2pt cos gs + t*)— P 


(1 —2pt cos ys + pt?) 


9 


ss Kq? 
2mt? (1 + a 


—t mkKp" 
_ +m(1+ t?)+2mk In q+ a mx (1.17) 
q 





K+] 


{x a a _ = py — 21 — pf _ (1—p*) (1 #*) (1 — pit’? 


(1—2pt cos gs+p7t?) (1—2pt cos ¢s+ p?t?)? 





k(p? — t?)? 
(p? — 2pt cos gs + t?) 





4d 
| + mie — 1) ip? — 2) — AP] (1.18) 
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respectively. Following this same procedure, the radial and tangential bending moments in region 
2, obtained from eqs (1.07) and (1.08) are 


Po(1+v) Po(1—v) 
M ,. = — ——— [VW ] 4. —— 
‘ 8amk [Ve] l6m7mKp* [2] 
and 
M:, = Pol + v) [W.] 


SamK 


Pol Ss : v) 
loam«Kp? 


[M2], 


respectively, where eqs (1.09) and (1.10) become 


=) 9 9,9 2 9 " (] t*)(] p*t*) 
{In (1 —2pt cos g, + p*t?)— k In (p? — 2pt cos ¢, + f*) — 
— | (1—2pt cos ¢4 p*t?)} 


‘ Kq° 

2me?(1+ 55) 
272 

ee 

m (1 2mx In p| 

K+ 


he. 2 _ ( l *t?) 2t7( | »7)* ( | “\(] t~)( | >? 7)? 
0» S ) (2 1) In (1 — 2pt cos ¢, + p?t?) 4 aa it e p eo 
s=1 (1 Zpt COS gs + prt-) (| 2pt cos + p*t?)?? 


o 2 2\2 ) 
: K(p t“) | 
} 


+ m(k — 1)(p? — t?)— 2m«p* + mKq |. (1.22) 


(p° —2pt cos gst t*) 
respectively. 


Since the expressions for the twisting moment, M,,, the tangential shearing force, Q;, and the 
radial shearing force in region 2, Vr, are independent of the load distribution, eqs (1.11), (1.12), 
and (1.14) are also appropriate for the case of the uniformly distributed load. However, Q,, of 
eq (1.13) being dependent on n, becomes 
0! Po = | k(p” — ft”) (i-Zr+erer 


? . * ) , B ‘ 9,9 
“' 4armk | & \(p? —2pt cos ¢.4 (1 —2pt cos gs + p?t?) 


i-oi~-seyr  } 
| = some > > mater 
(1 —2pt cos ¢, + p?t*)) 


for uniform loading. 


6. Stress at the Center of the Plate 


An examination of eq (1.11) shows that M,,=0 as a limit when p— 0, and is independent of 
the angular orientation. Therefore Mohr’s circle becomes a point, M;,,=M,,, and o,,=o7,. It 
is noted that M,, and M,, given in eqs (1.15) and (1.16), respectively, differ only in the sign of the 
second term of each equation, and that these terms, having p? in the denominator cannot be solved 
directly when p=0. However, they can be evaluated by setting M,,=M,,. Thus, eqs (1.15) 
and (1.16) yield 


' 


(l—p*t*y?—2e(1—p*) (1—p*(1—2*\(1 —p*e*? 





(x? — 1) In (1 —2pt cos os + pt?) + — 
, f (1—2pt cos gst+p*t?) (1—2pt cos y+ pt??? 


a(p* — *)? ) * 
Bais eee imix — i)(9*— 2 
(p° — zpt cos gs + t* 


9 


m«p* 
g 


)- 0 
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from which it is obvious that the second terms of eqs (1.15) and (1.16) are equal to zero when p= 0. 
The moments at the center of the plate (9=0) may now be written as 


P g 
2mkt* | 1 — 
Po(1+ v) ( A 
i] r =| SS ee, ee Yo! © 
1,,= Mi, —— 2mk In t+ ray 





+2mk In q— mk 


and the stresses on the surface are given by 


o2(1—-£ 
3P(1+v)| 7! (1 -) 
Or, =O1,=+ : +2 In 
4th? K+] 


:... 





p=0 
As a matter of convenience for graphical representation eq (1.25) may be written as 


3Po 
4arh? 


Or, =01,=— 


[é—A] 


” (l—vp), 
where €=(1+ vp) ln q+ or (1.27) 


and A=(14+ v)(14+In #)+(1 —v)e?. (1.28) 

Figure 2 is the graphical representation of eqs (1.27) and (1.28) for v=0.3. It is apparent 
from the A curve that the stresses at the center of the plate are also affected by the portion of the 
plate that overhangs the supports. The & curve shows how gq affects the stresses at the center. 
It should be noted that when g— 0, 0; and.o;—%. It is obvious that this theory cannot be used 
to compute stress at or very near a concentrated load or reaction point. 

An examination of eq (1.25) shows that the stresses at the center of the plate are independent 
of m and 6, and it can be shown that this equation reduces to 


_ 3Po (34 
” Sirh2 © " 





T 








qort 


FIGURE 2. Factors which affect the stress at the center of the plate. 
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FIGURE 3. Distribution of maximum principal stresses in the plate. 


for the case of a simply supported, uniformly loaded plate ((=q=1). Furthermore, for the case 
of a simply supported circular plate having a central load (¢q= b/c and t= 1) eq (1.25) reduces to 





oC 


(l+v)_ =; 4c? (1+ v) 


(1.30) 


2arh? 


| ] i $3 b? earl) 


Thus, it is shown that for these limiting cases eq (1.25) reduces to the Grashof equations. 

It should be noted that the usable strength of a plate is usually greater than that implied by 
eqs (1.25), (1.26), (1.29), and (1.30) [4]. These are elastic equations, and experience shows that 
local yielding of materials at a highly stressed point does not necessarily indicate structural failure 
or excessive deflection. Therefore, predictions of structural failure based on these equations 
would require the use of an appropriate theory of failure. 


Moments and Stresses Elsewhere in the Plate. To illustrate the stress distribution over a 
segment of a plate, figure 3 presents lines of equal maximum principal stress sensitivity (o0/Po) 
over one half the period of symmetrical distribution. This figure was prepared from the results ob- 
tained from eqs (1.11), (1.19), and (1.20) for the following conditions: m= 3, t= 0.8, g=0, v= 0.3, and 
h=0.125. It can be noted from the figure that the support is located on the ray @=0°, and that 
the larger values of the maximum principal stress sensitivity extend the farthest out from the center 
along the ray 0=z/m. These rays have further distinction in that an examination of the twisting 
moment, eq (1.11), shows that M,,=0 for any values of p and t along these angular orientations. 
It follows then that the bending moments, M, and M,, along these rays produce principal stresses. 


7. Comparison of Theoretical and Experimental Results 


The test results presented herein were obtained from retesting specimen A of a previous 
investigation [3]. Since the specimens, materials, and test methods were fully described in the 
previous paper, it suffices here to indicate that the dead load method of testing was used, and the 
tangential strains were measured along the support circle (t=0.976 and m=3) by means of foil 
type electrical resistance strain gages having a gage length of 1/8 in. These gages were mounted 
on the top and bottom surface of the specimen at the angular locations indicated in figure 4. 


8 








w THEORY 


10 





STRAIN SENSITIVITY, ‘, 











° 
re) 


SUPPORT 


LOCATION ALONG SUPPORT CIRCLE, pvecrees 


FIGURE 4. Theoretical and experimental strain sensitivities along the support circle. 


tangential strains were computed for the expression 


| ¢ 
é:, = E (1, — vor )=+ = ( Mi, —yp M,). (1.31) 


Inserting the appropriate values for radial and tangential bending moment from eqs (1.19) and 
(1.20) yields 


er 


3P (1 —v?) / OD» 
= WY, | a 
2 4ammEh?k« | ) 


2p?) 


where WV» and Q»2 are given by eqs (1.21) and (1.22), respectively. 

The agreement between the experimental data, which represents two separate sets of test 
results, and the theoretical results shown in figure 4 is very good. This, along with the good 
agreement reported for deflection [3], lends credence to the use of the theory to predict the flexural 
behavior of the plates. 


8. Summary 


The special application of Bassali’s more general solution presented herein provides the equa- 
tions necessary to calculate the elastic moments, shearing forces, and stresses anywhere in a thin 
circular elastic plate supported at points equally spaced along a concentric support circle and 
subjected to symmetrical bending. However, the maximum stresses at the center of the plate 
are probably of more general interest in the design and analysis of the thin circular plate subjected 
to a uniformly distributed concentric load. It is interesting to note that the equations take on a 
very simple form for these stresses and appear to be well adapted for design purposes. It is of 


9 





further interest that these stresses, eq (1.25), are independent of angular orientation and the num- 
ber of supports. As in the case of deflection [3], the effect of the annular portion of the plate over- 
hanging the support circle is evident as illustrated by the A curve in figure 2. 

Caution should be exercised in the use of these elastic equations for design purposes as local 
yielding may not be a valid design criterion for the particular material and structure under con- 
sideration. The selection of an appropriate theory of failure is of the utmost importance in this 
case. 


The good agreement between the theoretical and experimental strains presented herein along 


with that for deflection [3] serve to substantiate the ability of the theory to predict the flexual be- 
havior of the plates. 


The authors are indebted to W. H. Pell of the National Science Foundation for his advice 
during the initial phases of this investigation, and to D. R. Tate for his valuable suggestions. 
Special mention is due L. J. Davis for his active participation in the Laboratory. 
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For the accurate determination of internal pressures in solid-filled .piston and cylinder devices 


the effective area and its change with pressure are very significant factors. 


The computation of the 


change of area with pressure from elastic theory leads however to considerable uncertainties. 

The present paper describes how ultrasonic measurements made in solid-filled piston and cylinder 
devices to 45 kilobars are combined with measurements of the lengths of the samples as function of 
pressure and with ultrasonic measurements under hydrostatic pressure to yield values for cylinder 


expansion. 


The results are compared with values obtained from modified elastic theory. 


Key Words: Compressibility, cylinder distortion, high pressure, nickel, tellurium, ultrasonics. 


1. Introduction 


Several solid filled piston and cylinder devices with 
rotatable pistons have been used in this laboratory 
to determine transition ,pressures at pressures up to 
45 kbar. The internal pressure was computed as 
the ratio of ram force on the piston to effective area 
of the cylinder. Appropriate corrections were made 
for friction between moving and stationary parts and 
for distortion of the cylinder. An estimate of the 
change of the effective area of the cylinder with inter- 
nal pressure was calculated using equations devel- 
oped from elastic theory for long thick-walled cylinders 
subjected to internal pressure. These relations were 
modified using the solution to the tangential strain 
problem for an infinite cylinder with a semi-infinite 
zone of internal pressure [1].! 

The estimated uncertainty in this calculation of the 
area is two parts per thousand at 25 kbar and four 
parts per thousand at 40 kbar provided that the elas- 
tic theory on which the calculation is based accurately 
describes the problem. This uncertainty is sufficiently 
large to warrant some effort in a direct determination 
of the change of area of the cylinder with internal 
pressure. This paper describes such a direct deter- 
mination by means of ultrasonic measurements on a 
sample contained in the cylinder. 

The ultrasonic measurements are used to deter- 
mine the compressibility of the sample as a function 
of pressure. From this the length, which the sample 
would have in a rigid cylinder at any given pressure, 
is computed. This hypothetical length will be com- 
pared with a sample length obtained from a direct 


'Figures in brackets indicate the literature references at the end of this paper. 


measurement at the same pressure. The ratio of 
these two lengths is a measure of the expansion of the 
cylinder. 


2. Definition of Expansion Factor 


We assume a mean value for the expansion of the 
cylinder over the pressurized length (fig. 1) and we 
define as the expansion factor: : 


mean expanded area of the cylinder A,,. 


a(p)= : - 
F nonexpanded area of the cylinder Ao 





a(p) is a function of pressure. The relative change 
of the area of the cylinder with pressure is 


AA _ Aer A0_ op) —1 
Ay Ao in ; 
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Deformation of sample under pressure with and without 
deformation of cylinder. 








FIGURE 1. 





We further assume that the pressure is uniform 
throughout the sample, that the end faces of the sam- 
ple are plane, and that the cylinder remains straight. 

The expansion of the cylinder does not change the 
volume V(p) which the sample assumes at pressure 
Pp, so that 


V (p)= AoL ofl — Kp) =A, Ley(p) (2) 
where L,.,(p) is the length of the sample in the ex- 
panded cylinder at pressure p, Lo is the sample length 
at zero pressure, and K is the compressibility of the 
sample (see fig. 1). Whence follows 


a(p)= Aex _ Lol] —K(p) a. L(p) 


Ao Lex(p) 





L.x(p) 3) 


According to eq (3) the distortion of the cylinder and 


the change of area could then be obtained from a. 


measurement of the actual length of the sample and.a 
computation of the uniaxial compression of the sample 
in a rigid-walled cylinder based on the compressibility 
K(p) of the sample material. 


3. Determination of Sample Length /(p) 
for Uniaxial Compression 


At low pressures, where the compressibility of the 
sample placed in a rigid cylinder is constant, its length 
L(p) is 


L(p)= Lo — Kp). (4) 
Generally K is a function of pressure and, if this 
function is known, L(p) may be computed from 


p 
Lip)=Lo1— | K(p)dp). 


( 


Although data of compressibility versus pressure have 
been reported for a number of materials, the agree- 
ment between different sources for the same material 
is often very disappointing. It is desirable to use 
materials with small compressibilities to reduce the 
uncertainty in the determination of L(p). 

The choice of materials is also narrowed down by 
the fact that the present purpose requires a material 
with low yield strength to ensure a more nearly 
hydrostatic stress pattern. 

Typical uncertainties in L(p) caused by an uncer- 
tainty of 10 percent in the compressibility of the 
material are 0.1 percent for nickel at 20 kbar and 
3 percent for the more compressible tellurium at 
40 kbar. 

In 1957 R. K. Cook published a paper [2] dealing with 
the variation of the elastic constants and static strains 
with hydrostatic pressure and their calculation from 
ultrasonic measurements. The final equation reported 


in that paper relates the sample length at the hydro- 
static pressure p with the fundamental resonance 
frequencies of the same sample measured as a function 
of pressure for longitudinal and shear waves propa- 
gating along the axis of the cylindrical sample: 


Lo 


‘ 
where Lo is the length of sample at.zero pressure 


L(p) is the length of sample at pressure p 

ho= 4 Lipo 

fi, fe are the resonance frequencies for longitudinal 
and shear waves, and 


=" 
Kaa 
where Kagq is the adiabatic bulk compressibility 
Kj, is the isothermal bulk compressibility 
ho f{—sf%) is essentially the adiabatic bulk 
o 


modulus Bag= A>). 


at 
In the derivation ot eq (6) two approximations were 
made 


> 


(1 pki) 


“(1—pK(p)) 


and 


: _ 
[ K(p)dp 


L+A [? 
— K(p)d 
3 [. PrP (9) 


In the case of very compressible materials like tel- 
lurium the first approximation leads to a value of L(p) 
that is too large by a few tenths of 1 percent at 40 kbar. 
An error of similar magnitude but opposite sign is 
introduced if A in eq (6) is neglected. The second 
approximation can be avoided by rewriting eq (6). 
With less compressible materials like, for example, 
nickel, both corrections amount to less than 0.1 percent 
in L(p). This indicates again that materials with 
low compressibility are to be preferred for such 
measurements. 

Equation (6) permits the computation of the length 
of a sample subjected to hydrostatic pressure from the 
resonance frequencies obtained over the entire pres- 
sure range. The same equation less the factor 3 in 
the denominator in front of the integral holds for a 
sample contained in a rigid cylinder and compressed 
longitudinally. 





If the measurement is made in a nonrigid cylinder 
where the length of the sample is reduced by a~! due 
to cylinder expansion, eq (6) has to be corrected by 
multiplying ho=4. Lipo with a(p)~?. Since a(p) is a 
function of pressure it has to be written under the 
integral. Equation (6) for a nonrigid cylinder then 
becomes 


be 


1+A [?a?(p)dp_ 
L(p) 


hy Di. aie 4 ~ (10) 
Fs eas 3 | 
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a*(p) does not differ appreciably from unity. To 
facilitate the solution of eq (10) we assume that a(p) 
is a linear function of p 


dba 


a(p)=1+ Bp: as 


and we use an estimated value for B in the actual 
numerical evaluation of this and related equations. 

L(p) obtained from eq (10) can be combined with an 
L.x.(p) obtained from dial gage measurements of the 
sample length as described in section 6 to yield the 
expansion factor: 


Lo l 


Qi Pp) = Lop) iy 1 + A [’ _adp _ 
h J” 9 t » 
ta Pet | $ 


3 





Instead of using ultrasonic measurements in a solid- 
filled piston cylinder device to determine L(p), a sample 
length L4(p) could be determined under hydrostatic 
conditions with higher precision and then be con- 
verted to the hypothetical length L(p) with the help 
of the relation 


LB — Ls 
LE 


Li p)— ie 
—_. 


and the equation for a(p) is then 


¥2( p) = ——— 
aaa Lex(p) 


The superscript B in these equations denotes values 
obtained under hydrostatic conditions. 

Equations (12) and (13) are not very sensitive to 
errors made in the determination of the frequencies 
or the compressibility derived from them as long as 
materials with low compressibility are used. The 
main source of error is the determination of L.,(p). 
In a particular experiment to be described below an 
error of 25 wm in the determination of L.,(p) would 
cause an error in the change of area, AA, of more than 
30 percent. 


4. Determination of Actual Sample Length 
L.x{p) from Ultrasonic Measurements 


In equations (11) and (13) the sample lengths com- 
puted for the rigid-walled cylinder were combined with 
an actual sample length obtained from dial gage read- 
ings. If, from ultrasonic measurements under hydro- 
static conditions, both the shear and the longitudinal 
wave velocity (v) are known, these, together with the 
resonance frequencies or the times of flight (7) through 
the sample in the solid-filled cylinder, can be used to 
determine the actual length of the sample: 
L.x(p)=== UT. 


l 
y (14) 


‘ 


a. values for L,,(p) are obtained; one from shear 
waves and the other from longitudinal waves. The 
recision required in f is rather high. An uncer- 
tainty of less than 0.1 percent is necessary and 0.01 
percent is desirable at the highest pressures. 

The two values thus obtained for L,,(p) can be 
combined with the two values for L(p) used in eqs 
(11) and (13) in four different ways to yield: 
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Q6(p) = Loft 
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dp 
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The expressions in brackets are very close to unity 
and, especially with materials of low compressibility, 
do not contribute very much to the total uncertainty. 
Lo and L8, the initial sample lengths, can be measured 
with sufficient precision. The main source of error is 
in measuring the resonance frequencies. While the 
frequencies in hydrostatic environments can often be 
reproduced to better than 0.01 percent at 20 kbar no 
such claim can be made for our solid-filled cylinder 
assemblies where a scatter of the data of up to 1 per- 
cent was encountered. In planned future work a 





very definite improvement is expected from improved 
ultrasonic techniques and choice of more suitable 
materials. 


5. Description of Pressure Equipment and 
Pressure Determination 


A diagrammatic sketch of 
in figure 2. The cylinders used were made from 
cemented tungsten carbide with 6 percent cobalt. 
Radial support was provided by hardened steel sup- 
port rings with an interference fit. Axial support 
was provided by end-loading applied to the carbide 
cylinder through a bridge plate resting on extensions 
of ‘a thick cylinder wall of a ram with a 15 cm diam 
piston. This ram operated on the high pressure piston 
through a hole in the bridge plate. Rotation of the 
high pressure piston and 15 cm diam piston is accom- 
plished by a handle passing between the cylinder wall 
extensions. Both the end loading force and the 15 em 
diam piston force are supplied by the large press in 
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FIGURE 4. 


which this equipment is inserted. Transition pres- 
sure measurements were made in this equipment 
using a high pressure cylinder 5 cm long with a bore 
of 1.26 cm. Quartz transducers were mounted as 
shown in figure 2. In the 5 cm long cylinder a carbide 
filler rod was used to give the desired sample length 
for an entrance of the movable piston of about 0.1 cm 
at zero ram force. This additianal member in the 
path of the ultrasonic signal was undesirable. It 
was later eliminated by the use of a 1 cm long cylinder 
as shown in figure 3. 

Piston displacement was measured with a dial gage 
indicator mounted on the cylinder rim and bearing 
against an extension of the rotation handle. The dis- 
placement read is the sum of the movement of the 
sample end of the high pressure piston, the compres- 
sion of the high pressure piston, and the compression 
of the carbide pieces which back up the piston. The 
movement of the sample end of the high pressure 
piston is due to compression of the sample, expansion 
of the cylinder, ana, if used, the carbide filler rod. 

In the 5 em long cylinder a sample was compacted 
from tellurium powder with a particle size of 50 wm 
or less. The length after compaction to 40 kbar was 
0.6858 cm. Conical steel rings were used to prevent 
extrusion. The transition pressure with rotation 
(approximately two degree oscillation) was determined 
to be 40 kbar. In the 1 cm long cylinder a nickel 
sample, machined to fit, was used for measurements 
to 20 kbar. For the ultrasonic measurements the 
equipment was operated without rotation. 


The determination of internal pressure requires 
correction for friction and the determination of sample 
length requires correction of sample length for com- 
pression of various members. Figure 4 shows a typi- 
cal curve of uncorrected dial gage reading against 
uncorrected increasing and decreasing pressure. 
The dashed line to A is the calculated change in dial 
gage reading with pressure due to the compressibility 


dial gage reading 








+ 
40 kb 
pressure ——> 


Change of length of tellurium sample as indicated by 
dial gage reading and correction to internal pressure. 





of the carbide high pressure piston and its supporting 
pieces. A modified mean curve for the pressure on 
the tellurium is drawn to point F, the start of a tran- 
sition. The corrected pressure is obtained by follow- 
ing a line with a slope equal to that of line “A” from 
an observed point B or D to the intersection with the 
modified mean curve (points C or E) respectively. 
These corrected pressures are used with the observed 
frequencies for the ultrasonic computations. 


6. Determination of the Measured Sample 
Length, L,,(p) 


The measured sample length, L,,, is obtained by 
applying a correction to the dial gage reading for the 
compression of the various carbide members. 
pressures to 20 kbar a constant value of Young’s 
modulus for the carbide of 5.93 X 103 kbar was used. 
For pressures from 20 to 40 kbar a decreasing Young’s 
modulus was used. 

The effects of uncertainties in L,, on a@ are computed 
for a 5 percent uncertainty in Young’s modulus and 
for a 10 percent uncertainty in the estimate of effective 
length of the compressible carbides including the 
tapered piston and the backup pieces of various cross- 
sections. These effects on a(p) are larger in the 
tellurium measurements than in the nickel measure- 
ments primarily because of the greater length of com- 
pressible carbide in the tellurium measurement and 
secondarily because of the shorter sample length of 
the tellurium. These uncertainties in a(p) due to L,, 
are shown in figures 7 and 8 and will be discussed 
later. 


7. Derivation of Expansion Factor a(p) from 
Elastic Theory 
The expansion of the bore of a long thick-walled 
cylinder due to internal pressure p is 


Ad_p GU—-p)+U+p) 


= Cc 
dE 1 o- 


where E is Young’s modulus, pu is Poisson’s ratio, 
and W is the ratio of the outer diameter to inner diam- 
eter of the cylinder. 

The cylinder used with the tellurium had a wall ratio, 
W, of 4 and a length of 5 cm with only 0.686 cm or 
about one radius of that length exposed to internal 
pressure. 

A proportionality factor, F, to be applied to eq (19) 
for this pressure distribution is computed using values 
of tangential strain at the inner surface of a cylinder 
for a semi-infinite zone of internal pressure from the 
“Thick Walled Cylinder Handbook” [1]. For a long 
cylinder with internal pressure over a length equal to 
one radius starting at one end, the proportional expan- 
sion is computed by assuming a mirror image of the 
cylinder at the pressurized end of the cylinder. In 
Figure 5 curves (a and 6) are drawn for semi-infinite 
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FIGURE 5. Derivation of cylinder expansion from the distortion 
caused by reflected semi-infinite zones of internal pressure. 


at the original and reflected 
pressure discontinuities. Curve ¢ is the product of 
“a” and “b” and is an approximation for the factor 
for a pressurized zone of length 2 R in a long cylinder. 

A mean value of 0.75 is taken to apply to the expan- 
sion calculated from eq (19). a(p) is computed using 
the above factor and a Young’s modulus of 5.79 x 108 
kbar for pressures 0 to 20 kbar and a decreasing 
Young’s modulus for pressures of 20 and 40 kbar. 
The computed a(p) for the cylinder used with the 
tellurium sample is shown in figure 7. The effects 
of a5 percent uncertainty in F and a 33 percent uncer- 
tainty in F are shown. 

The 1 em long carbide cylinder used with the nickel 
sample has a wall ratio VW=3. Equation (19) is used 
neglecting the increase in effective W due to the steel 
rings and the variation of the expansion of this short 
section from that of a long cylinder. For a value for E 
of 5.97 < 108 kbar and a p of 0.25: 


pressure zones starting 


“c= 0.00503 at 20 kbar 
cf 


“i. a—1=0.0101 at 20 kbar. 


An estimated uncertainty of 5 percent in E results in 
a 5 percent uncertainty in a—1 and this is shown in 
the value of a@ in figure 8. No value is included for 
the uncertainty introduced by applying the theory for 
elastic distortion of long cylinders to a short section 
and neglecting the friction effects at the end faces. 





8. Experimental Determination of Cylinder 
Expansion Factors 


For the ultrasonic measurements in the solid-filled 
piston and cylinder assemblies the signal had to be 
generated outside the pressurized zone, radiated 
through the sample and picked up again outside the 
pressurized zone. Since plane and parallel interfaces 
along the ultrasonic path are necessary for all precise 
measurements, the signal could only be radiated along 
the longitudinal axis of the apparatus. Figures 2 
and 3 show how the 10 MHz quartz transducers were 
arranged. One transducer was used for shear, the 
other for longitudinal waves. Part of the signal origi- 
nating from a transducer is reflected at the tungsten 
carbide-sample interface, the rest travels through the 
sample and undergoes another reflection at the sample- 
tungsten carbide interface. The returning echoes plus 
usually several artifacts were picked up by the same 
transducer. The signal was then amplified and dis- 
played on the screen of an oscilloscope without 
demodulation. 

The delay time caused by the sample was measured 
in two ways: At high attenuations a short pulse was 
used and the time between the reflection from the 
front (tungsten carbide-sample) interface and the rear 
(sample-tungsten carbide) interface was measured by 
displaying the rf-pulse on the screen of an oscilloscope 
and using the variable sweep delay to let the pulses 
successively coincide with a mark on the screen. The 
transit times could thus be determined with a sensi- 
tivity of +20 ns. At low values of attenuation in the 
samples the phase comparison method developed by 
McSkimin [3] was used. The pulse length was ad- 


justed to provide sufficient overlap of succeeding pulses 
and the rf-frequency was adjusted for destructive 


interference. At constant pressure the critical fre- 
quencies could usually be reset to +1 kHz at 10 to 


12 MHz. 


t~e-, 


time of flight 





The electronic equipment consisted of an rf-cw- 
generator, a gated amplifier and, on the receiver side, 
a wide-band rf-amplifier and a wide band oscilloscope. 
The transducers were connected to the circuit through 
a tuned transformer bridge. 


The same electronic equipment was also used for 
ultrasonic measurements with a nickel sample under 
hydrostatic pressure. In this case only the phase 
comparison method was used. The sample was ce- 
mented between two fused quartz buffer rods with a 
shear and a longitudinal mode transducer attached to 
either end. 

Figure 6 shows a typical plot of transit time versus 
uncorrected pressure obtained in a solid-filled cylinder. 
The large amount of hysteresis is due to the friction 
between the piston and the cylinder. In section 5 we 
have already described how the uncorrected pressure 
is related to the internal pressure. This relation is 
now used to re-plot the ultrasonic data versus internal 
pressure. The result is, in the case of shear waves 
in tellurium, a single curve for increasing and decreas- 
ing pressure. In some cases the results were less 
satisfactory. Transit time data are read from the 
final curves at intervals of 1 kbar to be used for the 
computation of the cylinder expansion. 

The results of two such computations are shown in 
figures 7 and8. Figure 7 is for the long cylinder shown 
in figure 2 filled with sintered tellurium powder com- 
pacted at 40 kbar. The expansion factor computed 
from elastic theory is shown as a solid line. The un- 
certainty is schematically indicated by the vertically 
hatched area. Data for the expansion factor computed 
from ultrasonic measurements and measurements of 
sample length are entered as black dots. The total 
of the estimated uncertainties of the measured length 
is schematically indicated by the horizontally hatched 
area. The estimated uncertainty introduced by the 
ultrasonically measured L(p) is considerably smaller. 
Also entered in figure 7 are the adiabatic compres- 
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FIGURE 6. Time of flight for shear waves in tellurium as function of 
increasing, decreasing and corrected internal pressure. 
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FIGURE 7. Cylinder expansion with 
cylinder. 
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sibility A,q obtained from the ultrasonic measurement 
and the isothermal compressibility A,, obtained from 
the dial gage reading. The latter includes the effect 
of cylinder expansion and is therefore not a correct 
value. 

The overall uncertainty in this determination is 
extremely large. The greater part of the uncertainty 
is due to the determination of L,,(p) from the dial gage 
measurements. This fact emphasizes the need for an 
improved way to measure the actual length of the 
sample under pressure. 

Figure 8 presents the results obtained with a short 
cylinder filled with nickel. The uncertainties are 
indicated in the same way as in figure 7. The overall 
uncertainties are relatively much smaller. This is 
partly due to the improved design of the cylinder and 
partly to the smaller compressibility of the sample. 
The precision of the ultrasonic measurements was 
improved over that with tellurium. Consequently 
the disagreement between the measured and the com- 
puted values for a is much smaller. Again the com- 
pressibilities are entered. 


9. Conclusion 


The experimentally determined expansions for both 
the long and the short cylinder are larger than the 
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expansions computed from elastic theory. The same 
relation was reported by Bridgman [4] for his experi- 
mental determinations of cylinder distortion under 
various conditions. Therefore, a more rigorous com- 
putation of elastic distortion of short cylinders, includ- 
ing end effects and other modifications, seems to be 
desirable. While this mathematical solution does not 
seem to be very near at hand, the main effort will 
have to be in the direction of improved ultrasonic 
measurements with carefully selected materials in 
various piston-cylinder devices. The application of 
ultrasonic measurements to the determination of 
cylinder distortion as proposed in this paper will then 
lead to more meaningful and more precise data. 
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The absorbed dose in graphite from a cobalt-60 gamma-ray source was measured with a spherical 


air-filled cavity ionization chamber and with two spherical calorimeters. 


constructed from high-purity graphite. 


The instruments were 


The current per unit mass of air and the absorbed power 


per unit mass of graphite were determined with uncertainties of 0.40 and 0.17 percent respectively. 
When the two results are combined the value of W,;, * $n is found to be 33.72 electron volts in graphite 
per ion pair in air, with an uncertainty of + 0.14 electron volts. 


Key Words: Absorbed dose, calorimeter, cavity ionization chamber, 


comparison, graphite, W + Sn. 


1. Introduction 


Ionization techniques are often employed for the 
determination of radiation exposure and absorbed 
dose. Especially for photons in the megavolt region, 
the instrument of choice is a cavity ionization chamber. 
For absolute determinations of exposure with such a 
device the effective stopping power ratio, 5», contrib- 
utes the largest single uncertainty—about 0.5 to 2 
percent [1].!_ For absorbed dose determinations with 
a cavity chamber the principal contributors to the 
uncertainty are the stopping power ratio and the 
value of W, the average energy required to produce 
an ion pair in the gas. For air, the uncertainty in V 
is about 0.4 percent [1]. Higher accuracy for such 
determinations therefore necessitates a reexamination 
of the values of W ands». The present investigation 
in which measurements are made, both with the 
ionization technique and with the calorimetric tech- 
nique of the same beam of cobalt-60 gamma rays, 
provides data on W - Sy». 

When the necessary conditions of homogeneity in 
the material and uniformity of primary radiation are 
satisfied the relation [2' hetween ionization and energy 
deposited is given by the Bragg-Gray equation: 

P AT AE, . 
ee ta (1) 
VM m 


where — is the ionization current produced by the 
m 
) 


radiation per unit mass of cavity gas and = is the power 


M 


deposited by the radiation per unit mass of the ab- 
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' Figures in brackets indicate the literature references at the end of this paper. 


cobalt-60 gamma rays, 


. : - , . 
sorbing material. If — and y are given in amperes 


per gram of gas and watts per gram of absorbing 
material respectively, then W is to be given in elec- 


tron volts per ion pair. In the experiments to be 
2 


; ; ; : I 
described in this paper the values of 71 and 7, Were 
measured with small graphite calorimeters and 
an air-filled graphite wall cavity ionization chamber 
respectively, using cobalt-60 gamma radiation. The 

> 

quotient of the measured values ofa and * constitutes 
an experimental determination of the effective value 
of the product W - 5,» in terms of mass and electrical 
units. Such pairs of measurements have been made 
most recently by Bewley [3] and previously by Bernier 
et al. [4], and by Reid and Johns [5]. The present 
experiment repeats work that was done earlier in this 
laboratory [6] but with improved procedures that 
result in increased accuracy. 


2. Apparatus 


2.1. Radiation Source 


The source of radiation for these measurements 
was a cylindrical assembly of twelve encapsulated 
cobalt-60 pencils, each a half-inch in diameter and 6 in 
long, containing about 1000 Ci (curies) in all. The 
pencils are mounted upright on a steel base plate and 
are spaced equally around a 41s in diam circle of 
centers. This assembly is mounted near the bottom 
of a square concrete cistern 6 ft on a side in the lab- 
oratory floor. The source is covered by water to a 
depth of 10° ft. The ion chamber and calorimeter 
were mounted in a watertight can and, during measure- 
ments, were placed at the geometric center of the 





source. The exposure rate at that point was about 
55 roentgens per second. 

The exposure rate is fairly uniform near the center 
of the source. For example, the average exposure 
rate within a spherical volume 1.5 cm in diameter, 
which is the approximate inside dimension of the 
instruments, is only about 0.01 percent greater than 
at the center. 

The spectrum of photons inside the can was ex- 
amined and found to have approximately the compo- 
sition given in table 1 [7]. 


TABLE 1. 


Spectrum of cobalt-60 source 


| Fraction 


Component of 


| intensity 


Primary photons 0.80 
Singly scattered photons 18 
Multiply scattered photons 02 


Energy 


1.17 MeV and 1.33 MeV 
0.21 MeV to 1.33 MeV 
0.105 MeV (effective energy) 





This description of the source is clearly not unique. 
For example, Ritz and Attix, studying a nearly identical 
source, found that the responses they obtained with a 
set of energy-dependent chambers could be adequately 
explained by assuming a scattered intensity of 1.7 
percent at 169 keV and a primary radiation of 98.3 
percent [8]. 


2.2. Calorimeter 


The two adiabatic calorimeters used in these meas- 
urements each comprise two main parts: a spherical 
calorimetric body, or core, and a spherical enclosing 
shell. They are similar to one previously described 
[6, 9]. Figure 1 shows the essential dimensions and 
details of one of the present instruments. A thermo- 


FiGURE 1. Spherical graphite calorimeter with solid core. 

A Shell: Outer diam, 2.05 cm; inner diam, 1.28 cm; graphite. 

B Core: Diam, 1.00 cm; graphite. 

T;. Tz. Ts, Thermocouples: Constantan versus chrome! P; 0.003 in diam wires imbedded 
in holes with epoxy resin. 

H, Core heater made of electrically conducting resin. 

C Copper spots. A No. 43 copper wire was attached to each spot with epoxy resin. 

H, Shell heater. Karma resistance wire. Diam, 0.003 in. Attached to shell with 
epoxy resin. Space between wires is about 0.2 cm 


couple, 7), is inserted in a shallow hole in the core and 
securely fastened. It is used to measure changes in 
core temperature with respect to reference junctions 
held at a constant temperature. Thermocouples 72 
and 73, which are attached to the core and shell respec- 
tively, are connected in opposition and permit the 
measurement of temperature differences that may 
develop between the core and the shell during 
measurements. 

Two cores differing in thickness were used in order 
to be able to determine whether or not the attenuation 
of the gamma rays within the core was properly evalu- 
ated. One core was made of solid hemispheres that 
were joined together by a film, H,, of electrically con- 
ductive resin. The resistor thus formed measured 
about 100 9. It was used as a source of heat for cali- 
brating changes in potential of 7; in terms of energy 
input to the core. The second core was made of two 
hollow hemispheres having the same outer diameter 
as the first, but only about two-thirds the mass. These 
were joined in the same way as the first. The cores 
were carefully weighed at each step of assembly and 
the composition by weight of each component is given 
in table 2. 


TABLE 2. Composition by weight of calorimeter cores 


Solid core Hollow core 


Component 
Mass Mass 


Fraction Fraction 


Gram 
0.5946 


Gram 


Graphite 0.90585 0.9945 


Conducting 


0.9929 


0004, 0005 0000 0000 
0002, 0002 0001; 0003 
0009 0010 0007. 0013 
00024 0003 0001, 0003 
0034 0031 0052 


resin 
Copper spots 
Epoxy 
Copper wire 
Thermocouples 0031, 





Total 0.9108, 0.9999 0.5988, 1.0000 











The calorimeter shells were also made in halves that 
fit tightly together. The thickness was sufficient to 
stop electrons generated outside the shell by cobalt-60 
photons. Four hollow polystyrene spheres of about 
1.5 mm diam were cemented to the inside surface of 
the shell to hold the core in the central position. An 
insulated wire heater, H2, was wound over the outer 
surface of the shells and attached with epoxy resin. 
An adjustable current generated heat in the winding 
so as to permit control of the shell temperature. 
Control was automatic during the measurements, and 
kept the shell temperature equal to that of the core. 

The graphite of which the calorimeters and ion 
chamber were made was a reactor grade material, 
a sample of which left a residue of 0.014 percent when 
ashed in a muffle furnace. Calcium, vanadium, and 
titanium were the principal constituents detected by 
spectrographic analysis of the residue. This level of 
impurity has no significant effect upon the observed 
gamma-ray heating rates. The density of a sample of 
the graphite was 1.76 g per cm*. 
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FiGuRE 2. Electrical circuits used with calorimeters. 

(1) Core temperature measurement T,, core thermocouple; J, constant temperature 
reference cell; R, 10 Q manganin standard; POT, 0-1.6 V potentiometer; D.B., 5 dial 
decade resistor, 0 to 111,110 Q: AMP, contact modulated d-c amplifier; G, microammeter 
with zero at center of scale 

(2) Shell temperature control. 72, core thermocouple: 7s, shell thermocouple: J, con- 
stant temperature reference cell; AMP, contact modulated d-c amplifier; CONTROL, 
3-mode current adjusting controller; MAG AMP, magnetic amplifier; Hz, shell heater 
winding. 

(3) Electrical calibration. SUPPLY, constant current generator 
resistors: r, copper lead wires; H,, core heater 


R,. Rs, wirewound 


Three electrical circuits were required for operating 
the calorimeters. They are shown in figure 2 and 
functioned as follows: 

(1) Core temperature measurement. The EMF of 
the thermocouple, 7;, attached to the core was meas- 
ured by balancing it with an EMF developed by a 
simple Lindeck circuit consisting of a manganin 10 0 
standard resistor, R, a five dial decade resistor and a 
potentiometer set at 0.3 V. Adjustments in the decade 
resistance were made to provide a range of potentials 
of from about 45 to 175 wV across R. The null detec- 
tion system consisted of a contact modulated d-c 
amplifier and a galvanometer. The connections of the 
thermocouple wires to the coépper lead wires at J 
were kept at a steady reference temperature of 26.8 °C 
by immersing them in oil in the central tube of a glass 
cell containing solidifying diphenyl ether [10]. 

(2) Shell temperature control. The EMF of the 
opposed thermocouple pair, 7, and 73, was amplified 
by a second contact modulated amplifier in order to 
obtain indications of temperature differences between 
the shell and the core. During normal calorimeter 
operation the EMF was kept equal to-zero by the.action 
of a three-mode controller and a magnetic amplifier 
that regulated the current flowing through the heater 
winding, H», of the shell. The calorimeters were 
operated a few degrees C above 
the surroundings, so that the electrical heat supplied 
to the shell could be balanced by radiative and con- 
vective heat loss. 

(3) Electrical calibration. This circuit provided 
power to the core heater, H;, and so permitted. cali- 
bration of the core temperature indications in terms 
of added energy. Since it was observed that the core 
resistance increased slightly as the temperature: rose 


the temperature of 


during a measurement, a power stabilizing arrange- 
ment was employed consisting of a constant current 
source, and a calibrated stabilizing resistance, R;, 
shunting the heater. A calibrated resistor, R2, was 
also included as shown. A calibrated potentiometer 
was used to measure the voltages, V; and V2, that 
developed across R; and R». To give the desired 
stabilizing effect the resistance of R; was adjusted 
to be equal to that of the core heater, H,; for con- 
venience in measurement, the resistance of R» was 
adjusted to be one-half that of R,, so that V; was 
nearly equal to V2. 


2.3. Graphite Cavity lonization Chamber 


The important dimensions and details of the cavity 
ionization chamber are shown in figure 3. The wall 
thickness is equal to that of the calorimeter shells. 
The inner and outer surfaces of the chamber are 
spherical except at the end where the instrument 
is supported. The cylindrical collecting rod has a 
spherical enlargement at the end to reduce the elec- 
tric field strength in that vicinity and thus make the 
field strength more uniform. throughout the cavity. 

To obtain the volume of the cavity in the graphite 
shell it was weighed both when empty and when 
filled with distilled water. The volume of the cavity 
was then determined as the quotient of the observed 
mass of added water by the density of water at the 
corresponding temperature. Water was excluded 


























FIGURE 3. Spherical graphite wall cavity ionization chamber. 
Graphite wall assembly. Outer diam, 2.07 cm; inner diam, 1.27 cm. 

Graphite collecting electrode. Diam, 0.1 cm; diam of ball, 0.2 cm; length, 1.1 cm. 
Polystyrene high voltage insulator. 

Aluminum supporting stem with insulating plug to hold collecting electrode. Evacuated. 
Current lead to electrometer. 

Air vent. 





from the air vent with a plug during the weighings. 
Corrections were made for atmospheric buoyancy 
acting on the added water. We believe there was no 
significant absorption of water into the graphite be- 
cause the weight of the shell immediately after being 
emptied of water and wiped dry was not significantly 
different from the weight after the shell was kept 
overnight in vacuum. 

Five determinations of the cavity volume were 
made, including one in which the cavity was filled with 
mercury rather than water. The results are given in 
table 3. 


TABLE 3. Volume of ionization chamber cavity 


Cavity filling Computed volume 
cm? 
] 
4 


Water 1.0839 
Mercury 1.0773 
Water 1.0770 
Water 1.0771 
Water 1.0800 


+ 0.001, 4 


Average... 1.07, 


* Standard error of the average 


The volume of the collecting electrode was calcu- 
lated from its measured dimensions to be 0.0128 cm®. 
This was subtracted from the average determined 
above to give the volume of air within the assembled 
chamber. 

The ionization current from the chamber was deter- 
mined by measuring the corresponding voltage drop 
that it produced in flowing through a calibrated, wire- 
wound resistor. A calibrated potentiometer was used 
for the voltage measurement, with a vibrating reed 
electrometer for null detection. The collecting 
voltage was supplied by dry batteries. An aneroid 
barometer was used to determine the air pressure. 
A mercury thermometer was used to measure the 
temperature of the surrounding water. The humidity 
of the air was shown to have no significant effect on 
the chamber current [11]. 


3. Measuring Procedures 


3.1. Calorimeter 


In both gamma-ray and electrical calibration runs 
the length of time required for the temperature of the 
calorimeter core to rise between two arbitrarily 
selected temperatures was determined. The temper- 
atures corresponded to 160 and 60 uV respectively, 
as indicated by the core thermocouple. Since the 
thermoelectric power of the measuring thermocouple 
was about 60 uV per degree C, the temperature inter- 
val amounted to about 17/3 °C. The time required 
for the gamma rays to produce this change in temper- 
ature was about 40 min. 

The final determination of gamma-ray heating power 
by this method does not require that the end points 
of the temperature interval be known, but only that 
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they remain fixed during the set of runs. It was 
convenient, however, to divide the 100 yuV interval 
into fourths so that during electrical calibration runs 
five power observations could be made at nearly equal 
intervals of time. 

A calibration run was made as follows: The tempera- 
ture of the calorimeter was raised until a temperature 
indication of 175 wV was attained, and at that time the 
automatic shell temperature control was turned on. 
Several observations, usually seven, of microvoltage 
were made at intervals of 5 min to determine an initial 
zero-power drift rate. The electrical calibrating power 
was then turned on. When a temperature indication 
of 160 nV was attained (the controller having eliminated 
the transient shell-core temperature differential), an 
electronic timer was started, and measurements of 
V, and V2 were made. V, and V2 were also measured 
at 135, 110, 85, and 60 wV. The timer was stopped at 
the 60 uV indication. Calibrating power was turned 
off at 45 wV. Observations of microvoltage were then 
made at 5-min intervals as before, to determine a final 
zero-power drift rate. In reducing the data obtained 
in each run, the gross drift rate determined by the time 
interval between the 160 and 60 pV indications was 
diminished by the average of the initial and final zero- 
power drift rates so as to obtain a corrected net rate. 

For gamma-ray heating runs a similar schedule was 
followed. The calorimeter, placed near the corner 
of the pool where the radiation intensity was negligible, 
was brought to 175 wV and the shell temperature was 
put under automatic control. Measurements of micro- 
voltage were made for 30 min to determine the initial 
zero-power drift rate. The calorimeter was then 
inserted into the source. The timer was started when 
the temperature indication attained 160 wV, and was 
stopped when it attained 60 nV. At 45 nV the calorim- 
eter was removed from the radiation source to the 
corner of the pool and another set of measurements 
was made of microvoltage to determine the final zero- 
power drift rate. This measurement procedure was 
carried out with the calorimeter at different positions 
along the vertical axis of the source. Repeated 
measurements were then made at the position of 
maximum rate of temperature rise. 

In this method of calorimetry a systematic error can 
result if the thermocouple pair used for automatic 
shell temperature control do not measure the average 
temperature of the core and shell surfaces to which 
they are attached. To evaluate this error the calibra- 
tions and gamma-ray measurements were done both 
with the calorimeters evacuated and open to normal 
atmospheric air. In addition, measurements of the 
time constants for the relaxation of core-shell tempera- 
ture differentials with no shell heat were made under 
conditions of vacuum and atmospheric pressure. 


3.2. Cavity Chamber 


Preliminary measurements were made of chamber 
current versus chamber position along the vertical 
axis of the source to determine the position of maxi- 





mum response. At that position the current was 
measured over a range of collecting potentials up to 
540 V to provide the data required for extrapolation to 
infinite field strength. To eliminate extra-cameral 
currents each current measurement was made with 
both positive and negative collecting potentials. 

At the time of the measurement the atmospheric 
pressure was read from an aneroid barometer in the 
laboratory, and the temperature of the water was read 
from a mercury thermometer. The thermometer 
was brought to the surface for reading inside a 1-gal 
bottle of water. The effect of scattering and attenu- 
ation in the chamber neck, insulator and stem was 
determined by measuring the chamber current with 
and without an equivalent dummy installed on the 
opposite side of the chamber. 


4. Results 
4.1. Calorimetric 


The calibrations of the two instruments and the 
gamma-ray measurements obtained with them are 
given in table 4. The gamma-ray measurements, 
which were made during a period extending from 
January to July 1965, have all been corrected to a 
common date, March 3], 1965, using 5.24 yr as the 
half-life of the radiation. The standard error shown 
with each mean value was computed from the devia- 
tions of that set of readings from its mean. It can be 
seen that the response of these instruments was, in 
each instance, less when it was operated at atmos- 
pheric pressure than in vacuum. 


TABLE 4. Calorimeter measurements 


1 milliwatt 


8.096 
8.120 
8.164 
8.138 
8.158 


Mean 5.37 


y measurements, in microvolts per minute 


Hollo 


509 
503 


2.520 


2.507 2.548 
2.504 2.531 


2.491 2.544 


acuum | Air Vacuum 
| 


2.518 | 
2.522 2.509 


2.492 


2.493 2.547 


2548 


2.499 + 0.003 4 Mean 2.544 + 0.003 * 





Mean 2.515 + 0.004 2.526 + 0.011 4 


* Standard error of the mean 


The results of the temperature relaxation measure- 
ments are given in the four curves of figure 4, and 








FIGURE 4. Relaxation of temperature differences between core and 


shell of calorimeter with no electrical shell heat. 
4/4 is the fraction of the initial temperature difference. 


Initial temperature difference 
was produced by a pulse of heat in the core 


Scale of ordinates is logarithmic. 1, solid 


ore in vacuum: 2, hollow core in vacuum: 3, solid core in air: 4, hollow core in air 


show clearly the importance of air as a heat conduction 
path between core and shell. The relaxation rates 
(reciprocal of time constant) determined from these 
data are 0.168 and 0.237 per minute respectively, for 
the solid and hollow core instruments in vacuum; 
0.910 and 1.30 per minute respectively, for the solid 
and hollow core instruments in air. Since for each 
; . : , l l ; 
calorimeter this relaxation rate is A=k _ za) min@! 
41 42 
where /= core-shell heat transfer coefficient, joules 
per minute and °C, C,;=heat capacity of core, joules 
per °C and C.=heat capacity of shell, joules per °C, 
the relaxation rates obtained above are proportional 
to the corresponding core-shell heat transfer 
coefficients. 

As was noted from the mean values in table 4, both 
the electrical and gamma-ray heating rates were 
affected by the addition of air with its consequent 
alteration of the core-shell heat transfer coefficient. 
The correct heating rate that is desired in both deter- 
minations is that which would be obtained with no 
heat transfer to the core. That rate may be arrived 
at by extrapolating linearly to zero, on a heat transfer 
scale, the two heating rates measured in air and in 
vacuum. The extrapolations are shown in figures 
5a and 5b and give the values shown in table 5 for the 
calorimetric results. 


TABLE 5. Extrapolated calorimeter measurements 


Gamma-ray Calibration 


heating rate 


Calorimeter core 


wv min“ pV min-'mW- 
Solid 19 5.38 


2.5 
Hollow 2.548 8.213 
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FIGURE 5. Variation of calorimeter response with core-shell heat 
transfer coefficient. 

b. hollow core. The abscissae 

coefficients 


a, solid core: are proportional to the core-shell heat 


transfer 


From the values of the calorimetric results, shown 
in table 5, and the masses given in table 2, the gamma- 
ray heating powers per gram of the cores are obtained 
as follows: 

1) Solid C $s 519 10-3 
(i) Sond Gore: —=—Z.019 X ——— 
VM 5.381 


l 


x ——_ = 0.5140 x 10-3W 2g" 
0.9108? !* . 


2) Hollow C = sec 
(2) Hollow Core: va 0! 8.213 


] e ae 
< ——— = (1.5180 x 10-3 Wee!. 
x 05089 91 8¢ 1( zg 


4.2. lonization Measurements 


Measurements of ionization current were made at 
the position of maximum observed response in the 
source on December 31, 1964, before the calorimetric 
measurements and on July 16, 20, 21, and 23, 1965, 
after their completion. The chamber currents ob- 
served at a series of collecting potentials were extrap- 
olated to infinite field strength according to the method 
of Mie [12, 13] to determine the true ionization current 
in the absence of recombination. From extrapola- 
tions similar to those illustrated in figure 6 it was 
found that the collection efficiency at 540 V_ was 
about 0.993. The measurements with and without 
the dummy stem showed that absorption in the stem 
reduced the observed current by 0.13 percent. There- 
fore, the observed currents were increased by the 

1.0013 


0.993 © 
chamber was computed using the expression 


factor The density of the air in the cavity 


> 9729 


2410.4 


a * ' 
Pair = 0.001293 60) TT g/cm? 


where 


P=observed barometric pressure, mm Hg 
T= observed collecting volume air temperature °K. 
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FIGURE 6. Extrapolation of observed chamber current to infinite 
field strength. 


1, observed chamber current, amperes: V, collecting potential, volts 





In the air density calculation, 0.4 mm Hg was added 
to the observed pressure to allow for the difference in 
altitude between the barometer and the chamber, 
and 0.2 °C was added to the observed water tempera- 
ture to allow for the gamma-ray heating of the chamber 
as measured in a separate experiment. 

The mass of air was computed as the product of 
the volume of the cavity, 1.0663 cm*, and the density 
of the air at the time of measurement. The ioniza- 
tion current per unit mass of air was then determined. 
The observed currents and the results of the ¢omputa- 
tions are shown in table 6. The standard error of the 
average was computed from the deviations of the in- 
dividual determinations from the average. 


TABLE 6. Jonization measurements 


lonization current 
Current 
per unit 


Date of 
measure Air densityMass of ain 
Referred t nass of air 


3/31/65 


Observed* Extrapo 


lated ‘ 


ment 








1 
1.178, 
1.1785 





1.170 


1.168: 














1.551, 


+ 0.000, “ 


Average 


* At 540 V collecting potential 
» Corrected for stem absorption 


© Extrapolated to infinite collecting potential 


‘ Standard error of the average 


5. Analysis of Results 


5.1. Correction of Calorimetric Measurements for 
Impurities and Core Attenuation 


In the relation 
P 
] 
W: Sm >= 


m 


— refers to the cavity chamber’s response in amperes 
m 
2 


per gram of air in the cavity while refers to the 


f 
) 


calorimeter’s response in watts per gram of carbon in 


the core. For this equation to hold accurately, cor- 
rections to the calorimetric results are required to 
allow for the radiation absorption of atoms other than 
carbon in the cores and for the effect of gamma-ray 
attenuation in the cores. No correction is needed for 
attentuation in the calorimeter shells as they were 


equal in thickness to the chamber wall, and the addi- 
tional attenuation in the shell heater windings was 
negligible. ’ 

The effect of impurities is small and results in only 
a slight increase in gamma-ray heating relative to 
pure carbon. We calculate the excess, using the data 
of tables 1 and 2 and the mass energy transfer coefh- 
cients of Berger [14] to be 0.028 percent and 0.041 
percent for the solid and hollow cores respectively. 


In evaluating the core attenuation we have made use 
of data obtained by Loftus and Weaver [15] in this 
laboratory that show the variation with wall thickness 
of the response of a cavity ionization chamber. Three 
spherical chambers made of the same graphite as 
that used for the calorimeters were exposed in turn 
at the same distance from a cobalt-60 source similar 
in spectrum to the NBS water-shielded source. The 
resulting ionization current per gram of air within 
the cavity was determined for each chamber. The 
data of table 7 are the result of a series of such com- 
parisons, all at source-chamber distances greater 
than 1 m. 


TABLE 7. Attenuation in graphite 


Current per gram of air 


) 


Chamber ‘ 
relative to chamber 2 


Diameter | 


of cavity of wall 


Thickness | 





1.0065 
1.0000 
0.9900 


These data indicate a fractional decrease in chamber 
response of 0.0051 per mm of wall thickness.  Al- 
though the average wall thickness for each spherical 
chamber is larger than the radial thickness shown in 
table 7, the differences of these averages were found 
not to be significantly larger than the radial difference. 
The linear attenuation is therefore correctly given as 
0.0051 per mm. 

The length within the cores to be used for calcu- 
lating the attenuation was obtained from a theorem 
of geometry due to Tomkeieff [16] that the average 
length L of the intercept in a spherical body is 


b 
L=4-¢ 


where 


V=volume of sphere and 


S= surface area of sphere. 


If the body contains voids and V is taken to be the 
volume of the sphere less the volume of the voids, 
then the expression gives the average length of in- 
tercept not in a void. We find average gamma-ray 
paths in graphite from this theorem of 0.667 cm and 
0.438 cm for the solid and hollow cores, respectively. 

With no appreciable error of approximation the 
quotient of the average intensity of the radiation 





within the cores by the intensity outside can be taken 
to be 


1—0.0051 x 6.67 _9 9830 for the solid core and 


4.3 : 
1—0.0051 x +58 _ 0.9888 for the hollow core. 


The absorbed gamma-ray power in watts per gram 
of pure carbon in the absence of core attenuation is 


therefore found to be by the two calorimeter 
determinations: ' 


] i ] 


Solid core: 0.5140 x 7 
olid core: 0.514 1.00028 “0.9830 





x 10-8 
= 0.5227 x 10-° watts per gram. 


0.5180 x 


x 


l 
a IS 
0.9888 ad 


Hollow core: 1.00041 7 


= ().5237 X 10-3 watts per gram. 
Average: 0.5232 < 10-° watts per gram. 


re a 
5.2. Value of — = — 
Mem 


© 
The average values of M and 
work give for their quotient 
P 


C- T1551, x 10> 


] : . 
— determined in this 
m 


m 
in carbon per ion pair in air. 


5.3. Accuracy of Result 
P 


M 

T° 

m 
have considered the effect of uncertainties in the 
determinations of mass, ion chamber current, elec- 
trical calibration power, the attenuation correction 
that was applied to the calorimetry, the field strength 
extrapolation of the ion chamber current and in the 
extrapolated values of the electrical calibration and 
gamma-ray rates. 

The measurements of voliage in this experiment are 
subject to a systematic error of + 0.02 percent and a 
random reading error (standard error) of 0.01 percent 
which have been combined to give an overall uncer- 
tainty of 0.05 percent (three times 0.01 percent plus 
0.02 percent). 

The error in the ionization current measurements 
caused by an inability to null precisely is estimated 
to be no more than 0.04 percent. 


In estimating the uncertainty in the quotient 


The resistance value used in the ionization current 
measurements was considered to be accurate to within 
0.01 percent. The values of R; and Rz (fig. 2), used 
in the power calibration each have an uncertainty of 
not more than 0.05 percent. 

An overall uncertainty of 0.13 percent has been 
assumed for temperature measurements, based on an 
estimated systematic error of 0.07 percent and a stand- 
ard error of 0.02 percent, due to reading the instrument. 
No significant error is introduced by the method of 
measurement outlined above. 

Pressure measurements are believed to be accurate 
to one part in 10,000 and are subject to reading errors 
of no more than 0.02 percent. These have been com- 
bined to give an overall uncertainty of 0.07 percent. 

The uncertainty in the collection efficiency is esti- 
mated to be 0.96 percent, based on three times a 
standard error of 0.02 percent. 

It has been estimated that the core attenuation 
correction introduces a systematic error of not more 
than 0.05 percent. 

The mass of the calorimeter cores was believed to 
be known to two parts in 10,000. 

The uncertainty in the collecting volume of the 
cavity chamber has been estimated to be no more 
than three times the standard error of the mean value 
of the cavity chamber volume, from table 3, since 
the error contributed by the collecting electrode vol- 
ume measurement is n-gligible. Thus, an overall 
uncertainty of 0.36 percent has been assigned to the 
collecting volume. 

Measurements before and after each of the two 
series of ionization measurements indicated that the 
chamber was off axis by only a small fraction of a 
millimeter. We estimate that this introduced no 
significant error in the ionization measurements. 

The errors produced by the extrapolation process 
were computed to be no more than 0.10 percent and 
0.04 percent for the gamma-ray heating and electrical 
calibration rates, respectively. 

Tables 8 and 9 show the resulting uncertainties in 
a l 
y and —. 
The combined uncertainty in each ratio is the square 
root of the sum of squares of the individual uncer- 
tainties listed. The uncertainty in the value of 
reg . . , 
7 — has been computed to be (0.17 + 0.40)!? = 0.43 


percent. 


respectively, due to the different variables. 


P 
Uncertainty in 


TABLE 8. M 


T 


Resulting uncer 
> 


Variable —total uncertainty tainty in 
Vv 


Percent 
0.02 percent 0.02 
0.05 percent .10 
0.05 percent 05 
0.05 percent 05 
0.05 percent 05 
0.04 percent J . 04 
0.10 percent , 10 


Mass of core 

Power: Voltage measurement, V2, fig. 2; 
Resistance measurement, R,, fig. 2: 
Resistance measurement, Ro, fig. 2; 

Core attenuation factor; 

Electrical calibration rate extrapolation; 

Gamma-ray rate extrapolation: 





Square root of sum of squares 


26 





TABLE 9. Uncertainty in — 
m 


Resulting uncer 


Variable —total uncertainty tainty in 


m 


Percent 
Mass of air: Volume 0.36 percent 

Pressure: 0.07 percent 
0.13 percent 
0.01 percent 
0.05 percent 
0.04 perce nt 
0.06 percent 


Temperature 

Current: Resistance: 
Voltage measurement 
Null determination 


Collection efficiency: 





Square root of sum of squares 


5.4. Comparison With Other Measurements 


The previous measurements of this quotient that 
were mentioned in the introduction are shown in com- 
parison to the present result in table 10. These 
values should be directly comparable except for the 
slight differences in effective stopping power, Sin, 
that may arise from the spectral dissimilarities of the 
different cobalt sources used. The values of Sy listed 
in table 10 are those given by the authors as appli- 
cable to their sources. The value for the NBS source 
was computed from the expression: 


WE )i(E\dE + | - a ME) E)dE 
m(E) 


where i(£) is the distribution in energy of photon in- 
tensity table 1, w(E) is the mass energy 
transfer coefficient of carbon given by Berger [14] and 
Sm(E) is interpolated from the mean mass stopping 
power ratios for graphite relative to air given in table 
[A7 of NBS Handbook 85 [1]. While the methods 
used for obtaining 5, are not entirely consistent with 
one another, the resulting discrepancies are evidently 
small and probably negligible in comparison to the 
experimental uncertainties. 

The value of Q included in table 10 from the earlier 
NBS measurements [6] differs from the other entries 
in that the value of ionization per unit mass of air was 
derived from an exposure rate measurement of the 
source with a calibrated chamber. ° The original data 
have been recomputed to remove a correction pre- 
viously made for humidity and to change the attenua- 
tion correction so as to conform to that of the present 
experiment using the data of Loftus and Weaver. The 
uncertainty shown is largely that assigned to the NBS 
exposure standard in 1958. 

The present measurements give the value of Q and 
hence the value of the product W - 5, with an uncer- 
tainty of about 0.43 percent. Since the uncertainty 
in each of the factors is about one-half percent, the 
effect of the present determination is to reduce the 
product’s uncertainty by about a factor of two. Our 
estimates of error for the presently reported experi- 
ment are two-thirds and one-third those of Bewley and 
of Reid and Johns respectively. 


given in 


Table 10 includes the quotient of Q by Sm for each of 
the experiments. These should be equal, provided 
the values of 5 have been selected in a consistent way. 
The quotients show a reasonable degree of consistency, 
except for the early measurement of Bernier et al. 
If the data reduction computations of measurements 
3 and 4 were to be altered to eliminate the humidity 
correction that was included, the quotients of Q by 
5m would be increased, although probably by no more 
than 0.1 percent. 

The consistency of the values of Q+5m does not, 
however, imply a corresponding degree of accuracy of 
W, since each quotient carries with it the uncertainty in 
Sm. The uncertainty in 5, does not appear to be easily 
reduced by direct experiment or computation. 

\ more accurate experimental determination of W, 
on the other hand, does appear to be possible, and 
would, if available, also provide a more accurate value 
[OP Sm: 


TABLE 10. Comparison of results 


Measure 
Authors 


ment 
number 


Bernier et al. [4] 
Hart et al. [6] 
Reid and Johns [5] 
Bewley [3] 
Present 





| 
| 
! 
| 
| 
| 
| 
| 








‘Electron volts in carbon per ion pair in air 
Effective mass stopping power of carbon relative to air 


Electron volts in air per ion pair in air. Uncertainty does not include uncertainty in 5». 


The authors are indebted to Harold O. Wyckoff for 
his continued interest in and many valuable discussions 
concerning this measurement. 
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The resistances of three encapsulated, hermetically sealed, arsenic-doped germanium resistors 
from a commerical source rave been measured at every 0.1 °K from 2.1 to 5.0 °K in highly stabilized 
liquid helium baths; values of temperature were derived from vapor pressure measurements associ- 


ated with the liquid bath. 


Ten similar germanium resistors of common origin have been calibrated from 2.1 to 5.0 °K ina 


calibration comparator apparatus that is similar to a calorimeter’in design. 


Calibrations were per- 


formed against a resistor that had been previously referenced to temperature through a helium-4 vapor 


pressure calibration. 


Next, three of the resistors that were calibrated in the comparator apparatus were recalibrated 
in a liquid helium bath (at 2.2, 3.0, and 4.2 °K) and results from the two methods agree within + 1 milli- 


degree. 


The resistance-temperature data from both methods have been fitted to the polynomial 


m 
equation logio R= 3 1, -(logio T)" and the results of a computer program, which evaluates the co- 


n=0 


efficients, are presented. 


1. Introduction 


We have a group of impurity-doped germanium re- 
sistors that possess desirable thermometric properties 
at low temperatures — negative coefficient of tempera- 
ture, good sensitivity and excellent reproducibility 
under thermal cycling [1].2_ This paper describes the 
calibration (in the range 2.1 to 5.0 °K) of the above- 
mentioned resistor group. 

Several resistors have been calibrated directly 
against liquid helium-4 vapor pressures thus giving 
reference to temperatures from the “1958 Het Scale 
of Temperatures” [2]. Since this method of calibra- 
tion is costly in both time and refrigerant materials, 
a different procedure was utilized for extending the 
calibration over a greater number of resistors. An 
apparatus was constructed for comparison calibra- 
tions between “unknowns” and resistors that had 
been previously calibrated against helium-4 vapor 
pressures. 


2. Calibration in a Liquid Helium Bath 


2.1. Liquid Helium Bath 


The apparatus, used for calibrating germanium re- 
sistors in a highly stabilized liquid helium bath, has 
been described in an earlier paper [1]. While the 


'M. H. Edlow is presently associated with the United States Patent Office. 
2 Figures in brackets indicate the literature references at the end of this paper. 


previous report was based upon calibrations in the 
vicinity of 4.2 °K only, no modification of the apparatus 
was necessary for expanding the calibration range to 
include 2.1 to 5.0 °K. 

The bath consisted of about 8 liters of liquid helium 
which filled the spherical liquid helium storage Dewar 
(25 liters) to a height of 14cm. This provided a liquid 
helium surface level that did not extend above the 
thermally insulating vacuum jacket of the vapor pres- 
sure bulb. The storage Dewar was then connected 
to a vapor pressure regulator [3] and the surface vapor 
pressure monitored by a mercury manometer. Small 
amounts of helium gas from a 1-liter filling can were 
bled into the jacketed vapor pressure thermometer. 
After the helium gas pressure in the thermometer 
equaled the controlled bath surface pressure, 1/3 of 
an atmosphere of helium gas was slowly metered into 
the bulb to condense approximately 0.5 cm of liquid. 
Since the thermometer bulb volume was 2 cm*, the 
fluid level was well below the top of the bulb and the 
remaining volume was sufficiently large to contain 
some additional condensate-condensation of gas from 
the vapor pressure thermometer tubing occurred when 
the vapor pressure of the helium bath was intentionally 
decreased during the course of a calibration. 

The calibration from 2.1 to 5.0 °K was procedurally 
divided into three stages. Obvious checks were made 
to ascertain that there was appropriate electrical 
continuity of resistor leads and no undesired ground- 
ing in the resistance measuring circuit before the 
resistors were lowered into helium bath. Then the 
Dewar bath was cooled to appropriate temperatures 





by pumping through the vapor pressure regulating 
manostat and calibrations were performed at every 
0.1 °K in the range 4.2 to2.6°K. Next, the Dewar bath 
was warmed to the helium boiling point (by bleeding 
helium gas into the Dewar) and a second liquid helium 
transfer made. The second transfer was necessary 
since a height of only 3.2 cm of liquid helium remained 
in the Dewar after the calibration at 2.6 °K was finished. 
The bath was then cooled directly to 2.5 °K and subse- 
‘quent calibrations performed at the five lower tempera- 
ture points — 2.5, 2.4, 2.3, 2.2, and 2.1 °K. Thirdly, 
the Dewar bath was warmed to 4.2 °K, a necessary 
third liquid helium transfer made and the calibration 
performed in the range 4.2 to5.0°K. (Vapor pressures 
that exceed atmospheric pressure were obtained by 
closing off the manostat and allowing the Dewar vapor 
pressure to build up to a desired value at which time 
it was again manostatically controlled. The normal 
evaporation rate within the Dewar provided the gas 
for the increased pressure.) 


2.2. Vapor Pressure Thermometry 


The vapor pressure measuring techniques, used 
in these measurements, were consonant with pro- 
‘cedures that we have previously reported [4]. Since 
a temperature gradient exists at the surface of the 
“constant temperature” liquid bath, the vapor pres- 
sure thermometer tubing must be thermally insulated 
from the colder surface to avoid spurious effects caused 
by “cold spots” [4,5]. Consequently, the thermometer 
tubing wes surrounded by a vacuum jacket — the tubing 
length that otherwise would be in contact with the 
liquid surface —thus insuring that the vapor pressure 
thermometer indicated a vapor pressure that was 
associated with the bulk liquid. Vapor pressure 
measurements of both the helium surface and the 
thermometer bulb were routinely performed and re- 
corded. In a previous paper [4], we have listed the 
observed differences in these two measurements for 
tenth-degree intervals from 2.1 to 4.2 °K. Table | 
lists typical comparisons derived from the calibration 
data that are the basis for this paper. The differences 
between the two vapor pressures, P,; and Ps, are not 
unique since the pressure and temperature gradients, 
at the liquid helium surface, appear to be dependent 
upon at least two conditions—the geometry of the 
contained liquid helium and the Dewar evaporation 
rate. Nevertheless the qualitative change in magni- 
tude of the surface temperature gradient (AT) as a 
function of the helium bath temperature is supported 
by the results of many experimental “runs” in our 
laboratory. Additionally, as we have previously re- 
ported [4], calibrations of germanium resistors across 
the helium—A point have been “smooth” when tem- 
peratures in the helium-I region were derived from 
pressures measured in conjunction with the jacketed 
vapor pressure thermometer. We have been unable 
to obtain a “smooth” calibration across the A-point 
when temperatures are derived from surface vapor 
pressure measurements; as might be suspected, tem- 


peratures derived from either method are satisfactory 
in the temperature region of helium-II provided the 
effects of helium film creep do not enter as a variable. 

In the range 4.2 to 5.0 °K, the vapor pressure bulb 
required periodic additions of helium gas since, as 
the vapor pressure was increased, some of the pre- 
viously condensed liquid helium must be evaporated 
to supply the required helium gas in the vapor pressure 
thermometer tubing. The inverse operation was 
necessary as temperatures were reduced below 4.2 
°K i.e., liquid helium was occasionally removed from 
the bulb to avoid an overfilling of the bulb. 


TABLE 1. Vapor pressure measurements in a “constant” tempera- 


ture liquid helium-4 bath 


Temperature 5! AT* Temperature a— Pp," AT‘ 


mdeg 
2.0 
1.9 
2.0 
2.9 


Con uar-oS 


WNNNNNNNN 




















* P, is the helium vapor pressure of the jacketed thermometer 
» Ps; is the helium vapor pressure of the bath surface 
AT (the temperature equivalent of the 
(P,;—Ps) dP 
is obtained from the 


dP dT 
dT 


pressure difference, P;—Ps) is equal to 


‘1958 He* Scale of Temperatures” [2 


In the range 4.8 to 5.0 °K there were a few instances 
when the helium gas oscillated in the jacketed vapor 
pressure tube. These oscillations [6] were indicated 


by: (1) Thermometer vapor pressures greatly exceeded 
surface vapor pressures of the bath after surface pres- 
sure control had been established for several hours; 
(2) the evaporation rate of the helium bath, as measured 
with a flow meter, was greatly increased (by as much 
as a factor of three) but returned to its normal value 
when the bulb was evacuated: and (3) large oscilla- 


tions were observed in the mercury manometer. Re- 
liable measurements were obtained in these instances 
by altering the liquid bath temperature to a value 
where oscillations were not encountered. 

Resistance determinations, at each stabilized tem- 
perature, were measured potentiometrically; the 
electrical power applied to a resistor generally did 
not exceed 0.02 uW, but at temperatures above 
4.6 °K, it was possible to apply 0.1 wW without causing 
appreciable joule heating in the resistors. 


2.3. Discussion 


The calibration at each temperature in the range 2.1 
to 4.2 °K required approximately 3 hr—¥2 to 1 hr for 
cooling the bath to the desired surface vapor pressure, 
2 hrs for the bath to stabilize and 2 hr for performing 
electrical measurements. In one range above the 





boiling point, 4.7 to 5.0 °K, it required about 24 hrs 
to increase the bath temperature by 0.1 °K. There- 
fore, it was not unusual to expend three weeks in 
performing a calibration at every 0.1 °K from 2.1 to 
5.0 °K. The fact that at least three liquid helium 
transfers were necessary to cover the range made this 
calibration expensive as well as time consuming. 
For these reasons, an isothermal comparator was 
constructed for calibrating a large number of resistors 
more efficiently. 


3. Calibration in an Isothermal Comparator 


3.1. Comparator 


The comparator design and construction is shown 
in figure 1. The copper top A of the outside can B is 
6.04 cm in diameter and has an internal extended 
sleeve A’ (1.90 cm O.D., 2.22 cm long, and 1.27 cm 
wall thickness). The outside can B is mainly brass 
and has the dimensions of 6.67 cm diam, 21 cm long, 
and 0.159 cm wall thickness; the upper 1.27 cm of the 
can is a piece that was machined to receive the cylin- 
drical tongue of top A. Wood’s metal is used to solder 
the copper top to the-brass can. The bottom of the 
brass can C is a brass disk that has been silver sol- 
dered to the outside cylinder. During a calibration, 
the can B is contained in a liquid bath whose liquid 
level is well above the top A. The thermal anchor 
post T, 0.477 cm copper rod, is screwed into A. Three 
Teflon rods, each of 0.32 cm diam and 3.2 cm long, are 
equally spaced on the end of cylinder A’ and screwed 
into the wall thickness. A 0.64 cm perforated stain- 
less steel tube E was silver soldered to the copper top 
A; its opposite end was soft soldered to the top of the 
2.5 cm solid copper cylinder J. Twelve 0.36 cm diam 
holes N were drilled at an angle of 45 degrees to the 
axis of J and to a depth of 1.12 cm, which is sufficiently 
deep to contain the length of a germanium resistor 
encapsulation. The inner copper can K, 6.4 cm out- 
side diameter, 9.5 cm long, screws onto the threaded 
surface of the shoulder M. K serves as a radiation 
shield for the portion of the cylinder that contains the 
resistors. Several 0.16 cm holes O were drilled 
through K to enhance evacuation of the volume within 
K. F and F’ (0.64 cm, thin wall, stainless steel tub- 
ings), were soft soldered through A and serve respec- 
tively, as a helium gas filling or evacyation line and 
conduit for lead wires to the resistors and heater coil. 
A total of 14 leads (#38 A. W. G. formvar coated, 
double nylon insulated, copper wire) are used as con- 
necting wires—4 pairs as potential leads to the ger- 
manium resistors that are involved in a calibration, 1 
pair for current leads to the series connected resistors 
and 1 pair each for the sensor and heater coil. For 
the heater coil, #36 A.W.G. manganin wire was non- 
inductively wound on the bobbin H to provide a total 
resistance of 38,000 2 at room temperature. 

The seven pairs of connecting wires were conducted 
down from the Dewar cap through F’ (see fig. 2) and 
wound around the thermal anchor T — prior to wrapping 
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FIGURE 1. Design of comparator—the low temperature component 
of apparatus employed in calibrating germanium thermometers 
from 2.1° to 5.0 °K. 


around 7’, a small loop was made to allow for thermal 
contraction of the leads within the tube F’. Next, 
about 60 cm of the leads were spirally wound about the 
thermal anchor A’ and afterwards, coatings of G.E. 
varnish #7031 were applied to the wires and surfaces 
of T and A’ to enhance the thermal contact between 





FIGURE 2. A view of the upper portion of the cryostat that contains 


the calibration comparator. 


the lead wires and thermal anchors. This treatment 
will greatly reduce the heat leaks along the wires from 
room temperatures to the resistors. Additional 
lengths of wire were then spiralled around the three 
Teflon posts D and strung through holes D’ drilled 
near the end of the Teflon rods. One pair of leads 
was separated from the group and connected to the 
heater coil. The remaining 6 pairs were conducted 
through I, 0.159 cm holes that had been previously 
lined with Teflon, and thence spiralled around J 
where #7031 varnish was also applied. Next, the 
wires were then led through M, 0.159 cm holes also 
lined with Teflon, wrapped once around L and thence 
soldered to the germanium resistor leads with Wood’s 
metal. The soldered connections are enclosed by 
spun glass spaghetti to avoid electrical shorting or 
grounding of the leads. 

In use, the vacuum can B and its top A are exposed 
to a bath temperature that is less than 2 °K, while the 
cylindrical block J is maintained at calibration tem- 
peratures that exceed 2 °K. Thus, without proper 
precautions, large undesired heat-flows could occur 
and cause: (i) a temperature instability of the block J 
with resultant gradients; (ii) abnormally large evapora- 
tion rates of the surrounding liquid helium bath; and 
(iii) a lack of thermal equilibrium between the resistors 
that are to be calibrated. The comparator has been 
designed to reduce and control the flow of heat toward 


and away from J. (1). Radiant transfer from J or K 
to the outer can B is small. (2). The continuous evacu- 
ation of B minimizes heat transfer through gas con- 
duction and convection. (3) Heat conduction through 
the leads and mechanically supporting members has 
been reduced. 


There are several other component systems that 
are used in conjunction with the comparator. Tem- 
perature control of L or J (fig. 1) is effected by an auto- 
matic electronic system. A germanium sensor, 
mounted in a cavity N of L (fig. 1) serves as one leg of 
a Wheatstone bridge. Any unbalance of the bridge 
is amplified and fed into a controller which in turn 
regulates current flow through the heater H (fig. 1). 

Shown in figure 2, is the upper portion of the cryo- 
stat that consists of a glass Dewar A (7.0 cm I.D. and 
100 cm long), which contains the liquid helium in which 
the comparator is immersed, and an external Dewar 
that contains liquid nitrogen for thermally shielding 
the liquid helium Dewar. The figure shows a vacuum 
tight Dewar cap that accommodates the two tubes, 
which lead to the comparator (tubes F and F’ of fig. 1), 
and permits filling of the Dewar with liquid helium. 
Also shown, is the Dewar’s horizontal sidearm D that 
connects to a vacuum pumping line E—by “‘pumping” 
on the vapor pressure through this line, helium bath 
temperatures of less than 2 °K are achieved. The 
vertical tubing. which includes the vacuum valve, 
leads from an external “high vacuum” pumping sys- 
tem to the insulating vacuum space of the comparator; 
electrical lead ingress to the comparator is afforded 
by the tubing on the left. All of the apparatus, with 


the exception of the helium vapor pressure pumping 
system, is located within a shielded room that isolates 
the equipment from local rf interference [7]. 


3.2. Comparator Operating Procedure 


The resistors that were involved in a calibration 
were mounted in the copper block J (fig. 1) in a medium 
of Apiezon stopcock grease to enhance thermal 
contact. The outer can B was soldered in place with 
Wood’s metal and the 14 electrical leads then checked 
for electrical continuity. After appropriate leak de- 
tection the comparator was mounted in the cryostat, 
a small amount of helium gas was bled into the vacuum 
space (within can B) and the apparatus precooled to 
liquid nitrogen temperatures. Next, liquid helium 
was transferred into the inner Dewar, resistor lead 
continuity checks performed again, and then the liquid 
helium bath cooled to 1.5 °K. After the helium ex- 
change gas was sufficiently evacuated, the resistance 
of the Wheatstone bridge was set to correspond to the 
2.1 °K resistance value (determined in an initial cali- 
bration) of the temperature controlling sensor. The 
three control dials— proportional band, rate time and 
reset—of the servocontrol unit were set to obtain 
adequate temperature control of the comparator. 
Proper settings of the three controls result in a tem- 
perature control as fine as 1 X 10-4 °K at 2.1°K. This 
has been evidenced by repetitive resistance deter- 





minations for resistor #1 over a period of 30 min. 
Generally, once a good setting has been obtained at 
the lowest calibration temperature, the dials need not 
be changed for the remainder of the calibration. 

In establishing any given calibration control point, 
a value of resistance was dialed on the Wheatstone 
bridge and sufficient time was allowed for the recorder 
pen of the temperature controlling unit to reach the 
control or zero position. (Generally about 10 min was 
required to attain temperature stability after an in- 
crease of 0.1 °K in the temperature of the comparator.) 
Initial measurements on resistor #1 often indicated 
a deviation of approximately several millidegrees from 
a desired calibration temperature. In this event 
appropriate corrections were made on the Wheatstone 
bridge settings, five additional minutes allowed for 
equilibrium and measurements were again made across 
#1. This process was repeated until the resistance 
of #1 indicated a value that was within +0.5 mdeg 
of the desired calibration temperature. 


4. Germanium Resistance Thermometers 


For all of the germanium resistors, whose calibra- 
tions are being reported in this paper, we have a 
thermal cycling history. That is, the resistors have 
been thermally cycled from room temperature to 4.2 
°K and systematic calibrations of resistance and tem- 
perature performed at 4.2 °K. The methods of 
cycling and calibration have been previously described 
[1], and the reproducibilities ,of particular resistors 
were stated. The previous resistor identification 
symbols will apply to precisely the same resistors in 
this report. All of the resistors are four lead models 
and ‘are hermetically encapsulated in platinum con- 
tainers filled with helium gas. However, the cali- 
brator design and method of mounting resistors in 
the comparator are such that a leaking encapsulation 
will probably not affect the general validity of a cali- 
bration. Most of the resistors are similar in that their 
values of resistance are comparable at 4.2 °K—this 
probably means that their ‘““dopings” are similar and 
consequently, that similar equations will represent 
the resistance temperature calibration data. 


5. Experimental Results 


5.1. From the Stabilized Helium Bath 


Resistors 1, 2, and 3 (reproducible to approximately 
1 mdeg under thermal cycling [1]) have been calibrated 
simultaneously in the “constant temperature” helium 
bath from 2.1 to 5.0 °K at intervals of 0.1 °K; the results 


are listed in table 2. A log-log plot of resistance 
versus temperature for #1, #2, and #3, shown in 
figure 3, suggests the possibility of representing the 
calibration data by a polynomial of the form 


m 
logio R= yy Ay (logiol’)". 


n=0 


TABLE 2. Calibrated resistances obtained in the highly stabilized 
liquid helium baths 








Resistor #1 #2 Tempera- 


ture * 





2 °K 
12435. 2.0993 
10982. 2.2013 

9821.0 2.2995 
8839.0 2.3984 
7979.0 2.5003 


11749. 

10528. 
9491.6 
8588.0 


7841.0 
7106.4 
6611.8 
6119.1 
5681.2 


7271.7 
6572.9 
6108.0 
5639.4 
5227.1 


2.5994 
2.7130 
2.8005 
2.8996 
2.9986 


5287.2 
4938.8 
4625.8 
4345.1 
4090.2 


4855.8 
4527.8 
4232.8 
3968.7 
3729.9 


3.0997 
3.1997 
3.2997 
3.3996 
3.4999 


3859.6 
3650.4 
3457.4 
3279.7 
3118.8 


3513.4 
3317.0 
3136.9 
2970.6 
2820.7 


3.5996 
3.6996 
3.7999 
3.9007 
4.0000 


3660.6 
3448.2 


2972.0 
2829.8 
2699.2 
2579.6 


2468.6 


2683.4 
2551.8 
2430.2 
2318.9 


2216.2 


4.0981 
4.2000 
4.3014 
4.4011 
4.5002 


2498.7 
2381.1 


2272.1 


2363.2 
2266.3 
2173.4 
2086.6 
2005.4 


2119.0 
2029.0 
1942.8 
1862.7 
1787.9 


2169.7 
2075.1 
1985.4 
1901.4 
1823.1 


4.6009 
4.6998 
4.7999 
4.9004 
5.0005 














* Temperatures are derived from the “1958 He* Scale of Temperatures” [2] using values 
of vapor pressure that were obtained with a jacketed vapor pressure thermometer. 
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FIGURE 3. A plot of the resistance-temperature calibration data for 
resistors 1, 2 and 3. 


Temperatures were derived from liquid helium-4 vapor pressures. 





Consequently, the resistance temperature calibration 
data (obtained from the bath calibration) were pro- 
grammed on a computer for a least squares fitting to 
the above polynomial. The program was “‘cascaded” 
so that a succession of least squares fitting was per- 
formed for upper limits of m that extended consecu- 
tively from 3, . . ., 8. The results require discussion 
since there are at least two ways of considering the 
representation of the calibration data. 

(a) Except for the effect of spurious errors, the 
limiting accuracy of a calibration could possibly ap- 
proach +0.001 °K. This would probably occur at 
the lowest temperatures and would be mainly associ- 
ated with the accuracy of our helium vapor pressure 
measurements. 

(b) If, however, one feels that the 0.001 °K assess- 
ment of error is high or might be smoothed out in the 
polynomial fitting, it is not obvious as to how far one 
should proceed with the polynomial fitting. 


Because of these two considerations we felt com- 
pelled to present, in this paper, results which would 
be in accordance with either (a) or (b) of the pre- 
ceding paragraph. In the computer analyses for the 


3 
three resistors, the polynomial logio R= py A,(logio T)" 
n=0 
could possibly represent the data within the limits of 
condition (a). However there exists an obvious 
systematic deviation between the fitted function 
(Reaics T) and the input calibration data (Rgata. Taata)- 
This is apparent in figure 4 where points of inflection 
occur near 2.8, 3.6, and 4.4 °K. While the general 
conclusion from figure 4 might be the indication that 
the input data has been adequately fitted within the 
limits of experimental errors, another possibility is 
suggested—the data is more consistent than + 0.001 
°K and the function that would represent the data is 
3 


adequately R=» 

n=0 
A,(logio T)". Consequently we have continued the 
analysis to higher order polynomials. Table 3 con- 
tains the evaluated polynomial coefficients for resistor 
#1; a double precision computer program prepared by 


not approximated by  logio 


TABLE 3. 


Polynomial 


order 


Coefficient Coef ac Coef 8C * Coef. 6C * 


gg 


5.202 
— 4.249 
3.304 
— 1.593 


(0.0023) 
(.0145) 
(.0287) 
(.0185) 


5.180 

4.062 

2.732 
8355 
3687 


(0.0112) 
(.0922) 
(.279) 
(.369) 
(.179) 


5.380 

— 6.160 
11.348 
— 18.206 
16.846 
—6.715 


(.411) 
(1.68) 
(3.37) 
(3.33) 
(1.30) 





rrr? 


>> > > 


Std. Dev. 1.12 x 10-4 1.06 x 10-4 
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FIGURE 4. Differences between the fitted functions, log R= > A, 
é‘ : n=0 
(log T)" and the calibration data for resistors 1, 2 and 3. 

(a) The residuals AT = (Roen. — Roatald7/dR is the temperature equivalent of (Rogen — Rpata) 
which is the difference between the fitted polynomial function evaluated at Tpaia (Roen) and 
the resistance value of the calibration (Rpata. Tata) 
of resistance. 

(b) The residual plots apply to the three resistors [#1 ( . ), #2 (0) and #3 (+)]; the calibra- 
tion data for each resistor have been individually fitted to third order polynomials. 


dR/dT is the temperature coefficient 


the National Bureau of Standards Statistical Engineer- 
ing Laboratory was used to overcome the roundoff 
error characteristic of fitting high order polynomials, 
and in this case the roundoff error was several orders 
of magnitude smaller than other experimental 
uncertainties. 

As higher order polynomials are used, the syste- 
matic differences—between the fitted polynomial 
and the input data— gradually diminish up to the 7th 
order where apparently random residuals result. 
Figure 5 is the plot of these residuals resulting from 
fitting the seventh order polynomial for resistors 1, 
2, and 3. For n=7 of table 3, it is apparent that the 


m 
Coefficients of the polynomial, logiy R= Ss A, (logio T)", fitted to the calibration data of Resistor #1 


n=0 


Coef. 6C* 


(0.0396) 


5.685 

— 9,995 
31.136 

— 71.870 
> [97.546 
(— 70.554 





[20.767 }* 


(0.198) 

(2.48) 
(12.74) 
(34.41) 
(51.61) 
(40.77) 
(13.26) 


2.219 
40.873 
- 284.78 
1004.6 
- 2076.8 
2533.5 
— 1692.1 
477.6 





(0.735) 
(10.72) 
(66.3) 

(225) 
(454) 
(542) 
(356) 
(99.3) 





Je 


2435] 
[897] 
{— 103] 


(2.63) 
(43.7) 
(316) 
(1291) 
(3268) 
(5248) 
(5221) 
(2943) 
(720) 


5.2 x 10-5 





* 8C is the standard error of a coefficient 


» The brackets indicate that the enclosed coefficient is poorly defined; that is, the ratio of the coefficient to its standard error is less than 2. 
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FIGURE 5. Differences between the fitted functions, log R= + 7. 
=0 


(log T)" and the calibration data for resistors 1, 2 and 3. 


(a) See (a) of figure 4. 
(b) The residual plots apply to the three resistors [#1 (. ), #2 (0) and #3 (+)}; the calibra- 
tion data for each resistor have been individually fitted to seventh order polynomials. 


first coefficient is the least defined coefficient of the 
column. Since random residuals, exhibited in figure 5, 
were achieved by the seventh degree polynomial and 
the coefficients of the eighth degree are poorly deter- 
mined, we did not attempt to extend the fitting to 
higher order polynomials. (It does seem, for resistor 1 
in table 3, that odd order polynomials provide a better 
functional fit than even orders. This observation is 
generally applicable to other germanium resistor cal- 
ibrations between 2 and 5 °K and will be discussed 
later in the paper.) 

In comparing figures 4 and 5, one observes that not 
only has the systematic deviation of figure 4 been 
removed by progressing to the polynomial fitting por- 
trayed in figure 5, but also, the values of standard 
deviation that are listed in the last row of table 3 have 
been reduced by a factor of two. Despite this, one 
should be cautious about inferring: either (a), the cali- 
bration data were valid to + 0.0005 °K; or (b), the 
polynomial fitting is completely satisfactory; or (c), 
conclusions concerning the He* vapor pressure scale 
against which calibrations were performed. These 
items will require further clarification. 

Since resistor #1 was to serve as a basis for cali- 
brating additional resistors, a decision had to be made 
regarding its calibration representation. Although it 
is somewhat arbitrary, we decided to use the 7th 
order polynomial because it might have the effect of 
smoothing out experimental! errors in resistor #1’s 
calibration. The seventh order coefficients (table 3) 


dR ER . ‘ 
aT and TE for resistor #1 at 
intervals of 0.0001 °K from 2.1 to 5.0 °K; this generated 


table is the basis for calibrating additional resistors 


were used to generate R, 


It is pertinent to state that both 
PR 


a” ° Le e 
and TP? did not exhibit a sign change 
dT? 


against resistor #1. 
dR 
dT 
in the range of generation; that is, the fitted seventh 
order polynomial is a smooth function in the range of 
use. 


the generated 


the five calibration runs. 


In general the polynomial fittings of resistors #2 
and #3 produce results that parallel the findings for 
resistor #1. This is indicated by the plots of residuals 
for resistors #2 and #3 in figures 4 and 5. 


5.2. Comparison Calibrations 
Resistor #1 


With Reference to 


Ten resistors have been calibrated in the isothermal 
comparator with reference to the calibrated resistor #1 
and the data are listed in tables 4 to 8. In each of the 
indicated runs, three “unknowns” were calibrated 
against resistor #1 (the data for the third resistor is 
not included in table 4). As a control on the stability 
of the calibrations, resistor I was measured in each of 
Thus there appears in each 
table the measured resistances of resistors #1 and I 
and two unknowns (one unknown in table 4). Tem- 
peratures, which are assigned to the resistance values 
of resistor #1, are taken from the Resistance-Tem- 
perature generation that was derived from a seventh 
order polynomial fitting (sec. 5.1). 

It is also readily noted that particular calibration 
temperatures, in the 5 separate runs, are generally 
realized within a few millidegrees. This was inten- 
tional and eliminates interpolation inaccuracies that 


TABLE 4 


Run 3 


Calibration data for Resistor I and Resistor L 
R, re R, R, 


2 : 2 2 
13241 17117 15178 
11754 14905 13304 
10584 13192 11849 
9515.4 11652 10530 
8604.7 10365 9420.9 


7842.2 
7107.7 
6610.6 
6119.9 
5681.7 


9303.5 
8295.7 
7624.3 
6969.6 


8499.4 
7620.1 
7030.5 
6453.0 
5941.1 


5285.5 
4940.6 
4625.3 
4343.9 
4090.7 


3.1001 
3.1990 
3.2999 
3.4000 
3.4996 


3859.0 
3457.2 
3279.3 
3118.6 


2972.2 


3.6000 
3.7999 
3.9009 
4.0000 
4.0976 


3595.4 
3383.8 
3194.4 
3023.1 


2830.4 
2698.7 


4.1997 
4.3019 
4.4010 
4.5010 
4.6013 


2859.0 
2707.5 
2572.0 
2445.6 
2327.3 


2737.3 
2596.1 
2469.2 
2350.3 
2239.0 


4.6993 
4.8001 
4.9002 
5.0002 


2220.2 
2117.5 
2021.9 
1933.6 


2137.7 
2040.3 
1949.9 
1865.6 





2005.5 











®T is the temperature associated with the 
adjacent value of resistance for R,, which was 
previously calibrated against 73s [2]. 





TABLE 5 TABLE 7 


Run 5 Run 7 


Calibration data for Resistor H, Resistor M, and Resistor I Calibration data for Resistor F, Resistor J, and Resistor | 


. Seeeiecet 
R, t Ry Ry R, ye Ry R, R, 
+0] n 0] +] 0] A 2 0] 2 
13246.0 2. 14596 13819 17121 13238 2.1019 15290 15228 17092.5 
11742.0 2.20 12801 12165 14882 11745 2.2016 13364 13293.5 14873.5 
10579.0 2.29% 11425 10894 13182 10579.5 2.2949 11880.5 11807.0 13173.5 
9507.0 9% 10173 9725.7 11639 9509.8 2.3965 10532.0 10455.4 11632.6 
8609.8 z 9130.1 8754.7 10371 8608.0 2.4981 9408.7 9336.2 10355.8 
7840.1 2.5 8245.1 7924.5 9298.5 7843.9 2.5989 72 8398.6 9297.2 
7107.0 2.712 7402.4 7137.8 8292.8 7108.4 2.7124 757 | 7509.6 8286.5 
6609.6 2.8010 6845.5 6607.1 7620.4 6609.3 2.8010 79.5 | 6912.6 7615.6 
6121.2 2.8990 6293.9 6088.5 6969.8 6119.5 2.8993 6395.9 | 6332.4 6960.7 
5682.7 2.9985 5804.0 5623.7 6392.8 5680.5 2.9990 5880.2 5820.4 6384.4 
5287.8 3.0995 5366.1 5208.4 5879.1 
4941.5 3.1988 4981.9 4843.3 5433.9 4938.4 3.1997 5018.7 4965.9 5427.2 
4626.7 3.2995 4636.0 4513.9 5034.3 4624.8 3.3001 4661.1 4611.4 5029.8 
4344.7 3.3997 4329.6 4219.9 4681.1 4345.4 3.3994 4342.6 4296.0 4680.1 
4091.9 3.499] 4056.7 3956.7 4367.1 4089.7 3.5000 4055.4 4011.6 4362.6 


5287.2 3.0996 5422.9 5366.3 5875.4 


3859.6 3.5998 3806.8 3716.1 4082.2 3859.9 3.5996 3799.7 3760.0 4081.6 
3650.4 3.6993 3584.0 3500.1 3828.0 3649.8 3.6996 5567.2 3530.2 3825.6 
3458.5 3.7992 3379.7 3302.4 3597.3 3457.3 3.7998 2355.6 3322.3 3594.6 
3280.0 3.9005 3190.9 3118.9 3384.7 3279.7 3.9006 5162.5 3131.4 3382.9 
3118.9 3.9998 3020.9 2954.1 3194.5 3118.7 3.9999 2987.0 2959.6 3192.4 
2973.4 4.0968 2869.1 2805.7 3024.8 2972.1 4.0977 2830 8 2805.2 3022.7 
2830.1 4.1999 2718.8 2659.8 2858.5 2830.1 4.1999 2678.2 2655.6 2856.5 
2699.1 4.3015 2582.8 2526.6 2707.9 2699.4 4.3013 2539.8 
2579.1 4.4015 2459.8 2406.0 2571.5 2579.6 4.4011 2414.4 
2468.3 4.5004 2345.7 2294.1 2446.7 2467.9 4.5008 2296.9 





2519.1 2707.0 
2396.0 2571.2 
2445.1 
2363.9 4.5999 2238.3 2188.3 2328.9 
2266.7 4.6988 2139.6 2091.0 2221.1 
2173.6 4.7997 2044.6 1897.4 2117.8 
2086.4 4.9006 1956.5 1910.2 2022.0 
2005.1 5.0007 1874.4 1829.3 1933.3 


2327.3 
2266.4 4.699] ame 

2173.0 4.8004 1992.0 
- 2086.8 4.9001 1903.5 
2005.1 5.0007 1820.7 


919.9 
2116.6 
2022.1 
1932.9 














| 
| 
2363.1 4.6007 2187 8} 2 ] 1327 
| 
| 
| 
| 











T is the temperature associated with the adjacent *T is the temperature associated with the adjacent 
value of resistance for R,, which was previously cali 


value of resistance for R,, which was previously cali 
brated against 73. [2]. brated against 73, [2] 


TABLE 


might otherwise arise in comparing the five individual 
calibrations of resistor I. Table 6 possesses no cali- 
bration data above 4.2 °K because obvious experi- 
mental difficulties? which rendered the data invalid, 
9 occurred at calibration temperatures above 4.2 °K. 

| The calibration data for resistor I will be treated 
2 kK mM | Qo 0 separately in a later section of the paper. The data 


Run 6 


Calibration data for Resistor C, Resistor G, and Resistor I 


| ] 


ae | Seu fae 120495 hd for the other nine resistors has been subjected to 
10580.5 | 2.2049 | 11754 | 112765 | 13185.5 computer fitting and analysis for a polynomial of the 
9513.5 2.3962 10471 | 10031] 11649 . ; F 

8611.0 | 2.4977 | 9395.1 | 8986.7 | 10372.3 form 


7841.3 2.5993 8486.1 | 8105.1 9300.4 
7105.7 2.7129 7623.0 | 7268.8 8291.2 
6612.7 2.8004 7048.8 6713.2 7626.1 


m 
—_ : 
6121.6 | 2.9989 | 6481.7 | 6165.1 6971.2 log R= {1,(log yo7')". 
5681.5 2.9988 5976.7 5676.4 6390.3 t="0 


5287.1 3.0997 55 7.7 5 9243.6 5879.7 

4628.7 | 32008 | 47777 | 49001 3033.1 The determined coefficients and their uncertainties 
a) fae | oo. si gee are given in table 9. Again it is to be noted that: the 
ee eee ee calibrations of resistors in run VI cover a shorter tem- 
3650.0 | 3.6995 | 36925 | 34768 | 38279 perature interval (2.1 to 4.2 °K); and the polynomial 
ac) Ge | eae eee | soe fitting of this data would be expected to produce coefh- 
3119.1 | 3.9997 | 31122 | 2921.2 | 3194.9 cients that are significantly different from the coeff- 
2972.1 | 4.0077 | 2953. 2769.9 3023.3 cients representing the calibration of other resistors. 
2829.8 4.2001 27 2623.6 2858.0 


= = 4 














*T is the temperature associated with the adjacent ’ The experime -_ difficulty applies to run 6 (table 6) only, and was connected with the 
value of resistance for R,, which was previously cali- determination of the current that flowed through the resistor circuit during the 4.3 to 5.0 
brated against 73. [2]. °K calibration. 





TABLE 8 
Run 9 


Calibration data for Resistor K, Resistor E, and Resistor | 


R, fi Ry Ry R, 


2 
13251 
11750 


K 
2.1011 
2.2013 
10603 2.2928 

9519.2 2.3956 

8605.8 2.4984 


0] 
4099 
2377 
1073 
9845.9 
8827.4 


2 
16802 
14658 
13043 
11541 
10291 


0] 
17134 
14889 
13215 
11645 
10366 


7845.3 
7108.5 
6615.9 
6124.4 
5683.6 


2.5987 
2.7124 

.7998 
2.8983 
2.9983 


7979.3 
7167.2 
6628.7 
6094.3 
5616.8 


9302.0 
8293.0 
7626.2 
6971.4 
6411.0 6392.9 
5291.4 3.0985 
4946.7 
4631.6 
4350.9 
1096.8 


5198.0 
4827.1 
4492.5 
4195.7 
3928.4 


5908.8 
5470.8 
5075.0 
4725.7 


4411.8 


5883.6 
5440.2 
5040.1 
1688.3 
4372.5 





3865.3 
3656.9 
3465.5 
3281.9 


3120.3 3 


3686.7 
3469.9 
3271.4 

3082.2 

2916.5 


4129.0 
3876.4 
3646.0 
3426.8 


3235.7 


4088.7 
3835.6 
3605.0 
3386.4 


3195.7 


3.6961 
3.7954 
3.8993 


9989 





2971.1 
2827.7 
2700.5 
2580.8 
2469.4 


.0984 

2017 
3004 
4.4001 
4.4994 


2763.9 
2618.0 
2489.1 
2368.1 


2255.3 


3060.6 
2892.8 
2745.6 
2607.5 
2479.6 


3021.5 
2855.2 
2709.4 
2572.8 
2446.8 


4 
4 
4 


4.6000 
4.6988 
4.7995 
4.8999 
4.9990 


2150.2 
2052.7 
1960.1 
1874.0 
1794.5 


2360.2 
2250.1 

2145.3 
2048.9 
1958.9 


2329.3 
2220.6 
2117.9 
2022.6 


1934.7 








temperature 
for R 


associated with the adjacent 
value 
brated agains » (2 


which was previously cali 


Similarities are to be detected in comparing the general 
characteristics of the coefficients in table 9. The 
third and fifth order polynomial fittings yield well deter- 
mined coefficients while the fourth and sixth order 
coefficients, in some instances, leave much to be 
desired. The Ao’s of the seventh order are generally 
not statistically significant. 

The plots, for resistors K and E, in figures 6, 7, and 
8 are representative of corresponding plots for resistors 
C, F, G, H, J, L, and M. That is, the systematic 
deviations (between the calibration data points [Rgata, 
Taatal and the fitted functions [Reaic. Taatal) portrayed 
in figures 6 and 7 are most similar to residual plots of 
identical polynomial orders for the other seven resis- 
tors. The randomness of deviations for the seventh 
order polynomial (fig. 8) is also representative. 

It is interesting to note that, with the exception of 
run VI, the corresponding coefficients for all of the 
resistors are generally of comparable magnitudes and 
identical signs. This observation is most pertinent 
to the coefficients of the third and fifth order poly- 
nomials, where the coefficients are properly significant; 
and, excepting Ao’s, to the coefficients of the seventh 
order as well. 
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FIGURE 6. Differences between the fitted functions, 


3 
= > 
— 
n=( 


log R A, (log T)" 
) 


and the calibration data for resistors K, E, and I. 

(a) See (a) of figure 4 

(b) The residual plots apply to the three resistors [K (.), E (+) and I (0)]; the calibration 
data for each resistor have been individually fitted tothird order polynomials. 





AT(m®K) 








FIGURE 7. Differences between the fitted functions, 


log R= 


> A, (log T)" 
n= 0 


n= 


and the calibration data for resistors K, E, and I. 

(a) See (a) of figure 4. 

(b) The residual plots apply to the three resistors [K (.), E (+) and I (0)]; the calibration 
data for each resistor have been individually fitted to fifth order polynomials. 
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FIGURE 8. Differences between the fitted functions, 


7 


> An (log TP, 


log R 
; n= 0, * 7 
and the calibration data for resistors K, E, and I. 
(a) See (a) of figure 4 


(b) The residual plots apply to the three resistors [K (.), E (+) and I (0)]; the calibration 
data for each resistor have been individually fitted to seventh order polynomials 


5.3. Analysis of the Control Resistor Data 


The five individual calibrations of resistor I have 
been afforded a separate consideration since this 


TABLE 10. 


resistor served as a “control” through all of the re- 
ported calibrations. Its values of resistance, when 
compared with the corresponding values for resistor 
#1 (or more appropriately temperatures derived from 
resistor #1), indicated that all calibrations of resistor I 
were consistant within +0.0007 °K. While no pro- 
gressive drifting of the calibrations was observed (from 
Run III to IX), Run VII accounted for the largest 
portion of the maximum calibration deviation within 
the group of five runs; we are unable to ascertain 
specific causes for this occurrence. 

The calibration data for resistor I have been func- 
tionally fitted in the same manner as was described 
in 5.1. Each of the five calibrations was treated inde- 
pendently and the resulting coefficients are listed in 
table 10. (The statistical uncertainty of each coefh- 
cient is within the adjacent parenthesis.) Since the 
coefficients of Run VI (table 9) are based upon calibra- 
tion data that terminated at 4.2 °K, one must be 
cautious in comparing these coefficients with those of 
the other four calibration runs. Consequently, the 
general comparison between corresponding coefficients 
will be limited to Runs III, V, VII, and IX. Several 
conclusions are readily apparent. 

(a) Third order fittings—the four coefficients are in 
excellent agreement (within, of course, the limits of 


m 


Coefficients determined from fitting the polynomial logy R= S\ A,(logiwT)" to the calibration data of resistor I 
a 


n=0 


Pd dl a a al 


g 


Run Ill 
Resistor 18¢ 


0.0036 
02 


O4 


03 


15.14 
25.89 
24.94 
10.04 


{0.193} 
74.11 
195 
1722 
3526 
4266 
2828 
793 





x 10-5 





5 x 


».42 
+42 


Resistor 18¢ 


0.003) 
020) 
039) 


025) 


6 


5] 
0008 


6 


x 10-4 


10-5 





Run VI 
Resistor 


5.459 
4.713 
3.289 


1.416 
1.0 10-4 


5) 921 
5.265 
5.089 
3.985 
1.354 
x 10 
9.714 


7.407 


x10 


4.069 
14.536 


[828. ] 
[— 1628. | 
(41870. ] 
[— 1160. ] 
[299.] 


5.2 x 10-5 


16C 


(0.004) 
(.025) 
(.054) 
(.038) 


Run VII 
Resistor 
5.474 
4.815 


3.508 
1.569 


2.049 
45.98 
314.8 





“gis the standard deviation computed from the residuals resulting from the polynomial least squares fitting 
» 6C is the standard error of the adjacent coefficient 
] Statistically, the bracketed coefficients are poorly determined. 
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18c ” 


(0.0027) 
(.017) 
(.033) 
(.021) 





Run IX 
Resistor 


5.476 
4.823 
3.522 


1.577 
1.6 x 10-4 


5.454 
4.639 
2.962 
[— .835] 
{ 361 | 


1.6 x 10-4 
5.780 


~ 8.045 
16.936 


18c " 


(0.0034) 
(.0206) 
(.0409) 
(.026) 





Combined data 


Resistor Isc 
(0.002) 
(.011) 
(.022) 


(.014) 


9.473 
4.806 
3.489 
1.556 


8x 10-4 


(.008) 
(.070) 
(.21) 
(.28) 
(.14) 


(.035) 

(.37) 
(1.5) 
(3.0) 
(3.0) 
(1.2) 


7.409 
14.50 
24.43 
23.31 

9.34 


+x 10-4 
5.89 
9.61 
25.86 


55.27] 


[ 
[69.8] 


6. | 
[12.] 


4x 10-4 


1.857 
419.60 


(0.90) 
(13.1) 
342.1 (81.2) 
1200 (276.) 
2467. (557.) 
2995. (666. ) 
1990. (438. ) 


559. (1.22) 


1.3 x 





the stated statistical uncertainties). Furthermore, 
the ratio of an individual coefficient to its respective 
uncertainty indicates rather well defined coefficients, 
and the standard deviations of the individual least 
squares fittings are comparable. 


(b) Fourth order fittings—while similarities exist 
between some coefficients, two features clearly indicate 
that this order of polynomial fitting presents no ad- 
vantage. The A;’s are not defined, and the standard 
deviations, a, show little improvement over the third 
order of fittings. 

(c) Fifth order fittings —the coefficients, with respect 
to their statistical uncertainties, are clearly not as well 
defined as those of the third order polynomial: but the 
ratios are generally greater than six to one. More 
important, the standard deviations in polynomial fitting 
have materially improved over the third order thus 
indicating the superiority of the fifth degree for rep- 
resenting the calibration data. Here also, most of the 
corresponding coefficients i.e., Ao’s, Ass, are 
nearly identical within the limits of their statistical 
uncertainties. Those from run VII exhibit noticeable 
divergencies: several paragraphs earlier we pointed 
out that the calibration data fer run VII exhibited 
maximum variance with the other calibrations. 

(d) Sixth order fittings—clearly, a sixth order poly- 
nomial produces no beneficial result. 

(e) Seventh order fitting—while we should like to 
employ a seventh order polynomial for the data rep- 
resentation, the inconsistencies of the coefficients 
(Ay and A, in particular) for run III clearly imply the 
possibility of its impropriety. 

The polynomial coefficients for run VI, where the 
data is terminated at the upper limit of 4.2 °K, generally 
exhibit different qualitative trends, but are in quantita- 
tive agreement with the corresponding coefficients of 
‘the other runs for the fifth order polynomial. The 
general conclusion, which may be drawn from the 
comparison of the determined coefficients, strongly 
suggests that a fifth order polynomial offers a consistant 
representation for all of the calibration data (for resistor 
I) and that the standard deviation of such a fitting is 
about 9X 10-5. Deviations of the individual data 
points (Resistor I) from the fitted third order polynomial 
(run IX) are shown in figure 6. The curves’ significant 
features are duplicated, within a spread of + 0.0002 °K, 
for the corresponding representations of runs III, V, 
and VII. In the case of run VI, the general departures 
from the base line are reduced in size. 

Figure 7, which portrays the deviations of individual 
data points (Resistor I) from the fitted fifth order poly- 
nomial (run IX), is also representative of the other com- 
plete calibrations of resistor I. While it exhibits 
systematic deviations, they are much less than those 
exhibited in figure 6. A comparable plot (figure 8), 
associated with the seventh order polynomial, exhibits 
a random distribution: although this condition is desir- 
able, the indeterminacy of run III’s seventh order 
coefficients causes us to hesitate in advocating use of 
the seventh order polynomial for representing this 
calibration data. 








FIGURE 9. Differences between the fitted function, 


log R S A, (log T)", 


n 0 


and all of the calibration data for resistor I. 

(a) See (a) of figure 4 

(b) The residual plots apply to all of resistor I's 
IX) which have been fitted to a fifth order 


cally identihed 


bration data (Runs III, V, VI, VIL, and 


polynor Residuals for Run VII are specifi- 











FIGURE 10. Differences between the fitted function, 


z 


log R Ss A, (log T)" 


n=0 


and all of the calibration data for resistor 1. 
(a) See (a) of figure 4 


(b) The residual plots apply to all of resistor I's calibrated data (Runs III, VI, VII, and IX) 
which have been fitted to a seventh order polynomial. Residuals for Run VII are specifi- 
cally identified 


In summary, the preceeding paragraphs have in- 
dicated that through the five calibration runs: the 
reference resistor #1 and the control resistor I have 
exhibited relative calibration constancies that are 
within a temperature equivalent of + 0.0007 °K; addi- 
tionally, for I, the third and fifth order polynomial 
coefficients are consistantly determined within the 
limits of their computed uncertainties. 





Recognizing that the polynomial fittings of individual 
calibrations for I are somewhat restrictive in intercom- 
parisons between different calibrations, we have 
““polyfitted” all of the resistor I data. The resulting 
coefficients, listed in the last column of table 10, are 
generally not significantly different (within the limits 
of the coefficients’ uncertainties) from the coefficients 
of runs III, V, VII, and IX. Figures 9 and 10 are plots 
of deviations between the calibration data (Rgata. 
Taata) and the fitted polynomials (Reaics Taata) for the 
Sth and 7th order polynomials respectively. It is 
obvious in both plots that the data points from run VII 
consistently deviate from the rest of the calibration 
data and this’ accounts for the increased standard 
deviation that appears in the polynomial fittings of 
the last column. For each calibration temperature in 
figures 9 and 10 it is apparent that all of the calibra- 
tions agree within + 0.0007 °K. 

As a final criterion on the general validity of calibra- 
tions within the calibrator we have spot-checked the 
calibrations of run IX. Resistors I, K, and E were 
calibrated directly against helium-4 vapor pressures, 
in a constant temperature bath, at 2.2, 3.0, and 4.2 °K. 
Comparisons were made with the previous calibrations 
of run IX and the calibration differences (in tempera- 
ture equivalents) are shown in table 11. Overall, the 
agreement between the calibrations is reasonable 
although the discrepancies at 4.2 and 3.0 °K for resistor 
I are somewhat large. 


TABLE 11. Jntercomparison of calibrations from two methods * 


Temperature difference 


mdeg 





> | 
Resistor1® | Resistor K 
i 


Resistor E 





Method A—the resistors were calibrated within the cali 
brator against resistor #1 

Method B—the resistors were calibrated within a “con 
stant temperature” liquid helium bath against 7s (2) by 
means of vapor pressure-measurements 

» The positive sign (+) means, for the same value of resist 

ance, T (K) of method B is greater than T CK) of method A 

Previous thermal histories indicated that, under 86 
thermal cyclings between ambient temperatures and 4.2 °K, 
the 4.2 °K calibrations of resistors I, K, and E were repro 
ducible to 0.8, 1.1 and 1.0 mdeg respectively 


6. Conclusions 


We have shown that the resistance-temperature 
calibration (from 2.1 to 5.0 °K) for a particular group of 
ten similar germanium resistors can be represented by 
a polynomial function within the approximate limits 


of the calibration errors. In each case, based upon 

appropriate polynomial coefficients (fourth degree 

polynomials for resistors C and G, and fifth degree for 
dR 


, d?R ; ; 
the rest), generations of R, — and —3 as a function 


dT dT? 


of T (every 0.010 °K) have demonstrated that the deter- 
mined polynomials are smooth (within the range 2.1 
to 5.0 °K); the first and second derivatives do not 
change sign. While these conclusions are encourag- 
ing for the prospects of employing a fitted polynomial 
function to interpolate between calibration points, it 
should be stressed that the results apply to a select 
group of resistors. Also, caution should be employed 
in any attempts to use the polynomials for extrapola- 
tions beyond the calibration region. 

The type of polynomial that has been used in this 
work, 


m 


log oR = ‘ A, (log 102)", 


n=0 


has been selected somewhat arbitrarily. There is no 
theoretical basis for selecting this polynomial to repre- 
sent the resistors’ calibration data, and the fitting of 
the polynomial has been an empirical procedure that 
is applicable to the limited range of calibration data. 
The values of the determined Ao’s for any given resistor 
in table 10, indicate the inappropriateness of extrapo- 
lations based upon the determined polynomials. 
Additionally, when we have occasionally performed 
dR © , GR 
op and are 
mials) that extended beyond the calibration range, we 
found that the functions became erratic at external 
temperatures. 

We believe that the polynomial, which we have 
employed, is not necessarily the most appropriate 
function for representing the calibration data of these 
resistors. However it does serve as an interpolation 
aid that can be cautiously used until a more meaningful 
representation is achieved. In the immediate future 
we shall seek to find improved functions not only for 
the data presented in this paper but for more extensive 
calibration data (2.3 to 20 °K) that we have previously 
published [1]. 


generations (based on fitted polyno- 
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\ least squares technique has been applied to periodic earth temperature data for the purpose 
of determining basic characteristics of earth temperature cycles, such as thermal diffusivity, average 


temperature, amplitude, and phase angle of the temperature cycle. 


A new procedure was developed 


for obtaining a single thermal diffusivity which represents an average over time and depth at a particu- 


lar temperature site. 


This thermal diffusivity was obtained as a nonlinear part of least squares con- 
stants which yielded a best-fit harmonic curve to a given set of observed earth temperatures. 


The 


thermal diffusivity thus calculated and the calculation method developed are preferred to those obtained 
by current practice, which yields two thermal diffusivities, one based on amplitude decay and another 


on phase angle shift. 


Key Words: Earth temperature, least squares technique, thermal diffusivity of earth. 


Nomenclature 


2 {y= Least squares constants or best 
unbiased estimates of ao. ai, Q2. ay (°C) 

ao, a1, Qe. ay= Cosine coefficients of expected 
earth temperature function (°C) 

B,. Be. Bs, . By=Least squares constants or best 
unbiased estimates of 6,, bo. bs, . by (°C) 

b1, be, bs, by=Sine coefficient of expected 

, earth temperature function (°C) 

B,(x)=ith harmonic amplitude of earth temperature 
cycle at depth x (°C) 

C; or C;.x= Linearized variables (dimensionless) 

D=A least squares estimate of average thermal diffu- 
sivity, a (cm?/sec) 

F(0)= Surface temperature function 

Gi= Amplitude of ith harmonic coefficients = V a? + 6? 
(°C) 

M = Total number of temperature data (number of time 
data multiplied by number of depth data) (dimension- 
less) 

V,=Total number of cycles to be considered (dimen- 
sionless) 

M,=Total number of depths to be considered (dimen- 
sionless) 

N=Highest order of harmonics to be included for 

, expected earth temperature function (dimensionless) 

Pi(x)= Phase angle of ith harmonic of earth tempera- 
ture cycle at depth x, (radians) 

SS =Sum of squares of all of the differences between 
observed and expected earth temperatures 

S; or S;..= Linearized variables (dimensionless) 
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T= Period of cyclic data (31,557,600 sec for annual 
cycle) 

t= Earth temperature (°C) 

t,.= Earth temperature at x > & (°C) 

t;, = Observed earth temperature (°C) 

x == Depth in earth (cm) 

x; = Observed depth in earth (cm) 


Greek Characters 
a= Expected thermal diffusivity of earth (cm?/sec) 


iT ss 
= (ina) 


aT 


Bi= Temperature decay factor = 
= tee : li 
Bi = Estimate of Bj= yee (em ~?) 
aie DT 
_ OSS 
aD 
@=Time coordinate used in the text, measured from 
January 1 (sec) 
0,,= Observed time with datum point at January 1] (sec) 
®;= Phase angle of ith harmonic of earth temperature, 


-1 bi 
aj 
A=Root mean squared deviation (defined in the text). 


(radians) ®;=tan- 


1. Objective 


The purpose of this paper is to discuss a harmonic 
analysis of cyclic temperatures at selected depths of 
earth. The existing techniques employing a simple 





harmonic equation [Kusuda, 1965; Penrod, 1960] or 
Fourier analysis [Carson, 1963] have the disadvantage 
of yielding different thermal diffusivities for a given 
set of data, depending upon the choice of data reduc- 
tion method, choice of cyclic period, choice of two 
depths, and the choice of computational method. The 
technique introduced in this paper eliminates this 
disadvantage and it is statistically more sound than 
previous methods for the analysis of periodic depth 
temperatures of earth. Furthermore, the method 
derived herein may be applied for the accurate calcula- 
tion of thermal diffusivity of the surface layer of other 
celestial bodies in the:event cyclic patterns of their 
depth temperatures are obtained in the future. 


2. Introduction and General Background 


Numerous investigators in the past have studied 
earth temperatures with reference to the earth’s ther- 
mal diffusivity. A recent publication of Kusuda and 
Achenback [1965] includes an extensive bibliography 
on the subject. Most of the analyses in the papers 
listed in this bibliography apply the well-known peri- 
odic heat conduction formula for either diurnal or 
annual (or both) cycles of observed earth temperatures. 

These analyses have been made to establish the 
average, the amplitude, and the phase angle of cyclic 
temperature patterns at various depth levels, so that 
earth temperatures may be predicted by a computa- 
tional method. The amplitude decays and phase 
angle shifts of temperature cycles with respect to 
depths have been used to estimate the thermal dif- 
fusivity of the station where the earth observations 
were made. The analysis of periodic earth tempera- 
tures, therefore, eventually ends up with four basic 
parameters: thermal diffusivity, average temperature, 
amplitude, and phase angle of the surface temperature 
cycle from a given datum time point. 


The investigators have encountered the following 
problems during the course of their studies. (1) The 
earth surface was difficult to define and earth surface 
temperatures were subject to disturbances from cli- 
matic changes above the surface, such as rainfall, 
cloud cover, evaporation of water, and wind velocity. 
(2) The actual earth temperatures observed did not 
follow an ideal simple-harmonic pattern. Moreover, 
their cyclic patterns varied from year to year and from 
day to day. Except in the work of Carson and in that 
of Kusuda and Achenbach, some degree of arbitrary 
human judgment was exercised to establish meaning- 
ful periodic patterns. For example, when several 
years of monthly average earth temperature records 
were available, the annual thermal diffusivity was 
calculated for each set of annual data. Such calcu- 
lations imply an annual change of thermal diffusivity, 
which is a misapplication of basic heat conduction 
theory, as explained later. (3) For a given station, 
thermal diffusivities were computed both from the 
difference of logarithmic amplitudes and from the 
difference of phase angles for two different depths. 
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Unfortunately, the diffusivities calculated by these 
two methods usually did not agree, differing by a 
factor of two in some instances [Kusuda, 1965]. (4) 
When the thermal diffusivities were computed at two 
different pairs of two-depth levels at a given tempera- 
ture station, two different estimates of the diffusivity 
were obtained, even when a consistent method of 
diffusivity computation was employed. 

These four problems resulted in uncertainties of 
thermal diffusivities due to time, due to depth, due to 
method of computation, and due to the method of 
data reduction for a given set of earth conditions at a 
given locality. 

Most of the investigators also attempted to predict 
earth temperatures based on an average of thermal 
diffusivities thus computed, and on other basic param- 
eters; namely, the earth surface temperature ampli- 
tude, the phase angle of the earth surface temperature 
cycle measured from a certain reference datum point, 
and the average earth temperature during the cycle. 

Mathematical description of an earth temperature 
cycle by a Fourier series technique, similar in nature 
with the least squares technique, has been tried by 
several authors. The most recent and comprehensive 
description was made by Carson [1963] for the soil 
of Lemont, Ill. The Fourier series technique is basic- 
ally applicable only to a single complete cycle where 
data points consist of a unique set of temperatures 
for equally incremented time coordinates. Carson 
applied the technique to each individual one-year set 
of monthly average earth temperature cycles over a 
period of three years, resulting in three sets of Fourier 
coefficients, which did not agree with each other. 

Kusuda and Achenbach [1965] applied a least 
squares technique in fitting observed earth tempera- 
ture data to simple harmonic expressions for the 
purpose of computing the four basic parameters of 
earth temperature cycles. 

In their least squares technique, Kusuda and Achen- 
bach [1965] incorporated several years’ records simul- 
taneously, thus treating the temperature data as 
statistical variables or variates rather than functional 
variables, which is the case in a Fourier analysis. The 
least squares technique fits a sinusoidal curve to the 
observed data in such a manner that the variance 
calculated from the residual deviations of the observed 
temperature from the fitted curve is a minimum. The 
number of data points, the time intervals for the tem- 
perature observations, and the number of data for a 
given time are not restricted as they are in a Fourier 
analysis. 

Kusuda and Achenbach [1965], however, began with 
the premise that the earth temperature cycles follow 
a simple harmonic mode. This premise may be 
reasonably applied to annual cycles of monthly average 
earth temperatures, but it is not quite valid for diurnal 
temperature cycles. They, however, still obtained 
two thermal diffusivities based upon ,the amplitude 
decay and phase angle shift techniques, same as those 
used by other investigators. 

The extension of the least squares technique to cover 
the higher harmonic terms is also desirable if the 





method is to be applicable to diurnal earth temperature 
cycles of the earth theoretically differing from the 
simple harmonic pattern because of a complex non- 
linear surface heat exchange with the sun and the 
surrounding space. 

In the following section the basic heat conduction 
relation of the periodic temperature cycle is first re- 
viewed. The least squares technique is then applied 
to annual cycles of monthly average earth temperatures 
as an illustrative example. 


3. Basic Heat Conduction Theory on a Semi- 
Infinite Solid 


The surface layer of earth may be assumed to be a 
semi-infinite solid having a homogeneous thermal 
property, as far as the heat conduction theory is con- 
cerned. This assumption, however, may be inade- 
quate in some cases, and a more refined model, such 
as composite solids (an upper layer of finite thickness 
over an infinitely thick solid bed of different thermal 
properties), is desirable. As a theory, this composite 
model probably represents more closely the actual 
surface of earth than the homogeneous model. But 
the thickness of the top layer, and the variation of 
properties between the top layer and the bottom bed 
are usually unknown, which makes the composite 
layer nodel less meaningful than the homogeneous 
model. Also, the mathematical complexity involved 
for the composite solid model may not be justified in 
finding an accurate thermal diffusivity of a surface 
layer from the observed depth temperature cycles. 
Hence, the homogeneous semi-infinite solid model is 
used to calculate average thermal diffusivity repre- 
sented by the surface region of the earth. Or, con- 
versely, temperature cycles at various depths are 
approximated by the equation derived for semi-infinite 
heat conduction models of homogeneous thermal 
properties. 

The heat conduction of the semi-infinite homoge- 
neous solid can be mathematically expressed as 
follows: 

=o (1) 
Ox” O06 


at x > % t— t,, 
at x= 0 t=F(@) 
and F(6+ T)= F(6) (2) 

In actuality it is difficult to define a unique F(@) 
that represents daily or annual cyclic observations of 
earth surface temperatures, since the character of 
these cycles varies from one period to the next. This 
variation is due to the varying nature of heat exchange 
over the earth surface with respect to climatic condi- 
tions. It is postulated that other celestial bodies 
without atmosphere may experience a much more 
consistent and stable F(0@) than the earth experiences. 
In this paper it is assumed that there exists a unique 
F (6), which may be a statistical average of many 


periodic cycles. Since F (0) is a periodic function of 
period 7, it is possible to expand it as follows: 


F(0)=ay+ > aj COs (2 0)+ bs b; sin (+ 6) (3) 


The values of ao, a1, a2, . . . bi, bs, b3 are dependent 
upon the actual nature of the assumed surface tem- 
perature cyclic pattern. 

The cyclic earth temperature [Eckert, 1959] at depth 
x that satisfies (1), (2), and (3) is 


Jin = im 
al cos T V2 T x) 
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(It should be noted from eqs (3) and (4) that Ao=t,, 
Formula (4) can also be expressed as 


“a T 2710 F 
t=t, + S' Gie Yi+cos (“2_@,— / x) 
V a 7 4 


jaa aT 
where 
Gi=Vae+ Bb 


b 
tan ® =— 
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By letting 


Gie & ‘= Bi(x)= depth amplitude 


®,—./ x =Pi(x)= se angle 
D; Vor Pi(x)= depth phase angle, 


the thermal diffusivity @ can be readily calculated by 
knowing Bix) and P(x) at two different depths x; and 
x2 by the following formulas: 


log = 


= Bilxs) 


api = (9) 


itr X1— X2 2 
= A ; p= 15.2. 4: 2 ee 
7 | P(x) “FI 

In formulas (8) and (9), notation ag and a@p are used 
instead of @ to signify that these two thermal diffusi- 
vities are computed by amplitudes and by phase angles 
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of ith harmonic constants, respectively. Theoreti- 


cally, of course 


Qpi= apj= aa for i= 1, 2, ae aS 


(10) 


Unfortunately, however, these relations have seldom 
been found in practice. Most previous researchers 
have made extensive studies on ag; and ap, only, ig- 
noring higher harmonics. Even for the first harmonic 
constants, agreement between ag and ap; was poor. 
Not only that, the values of ag; and ap; varied, as 
mentioned before, depending upon the choice of x; and 
x2 of formulas (8) and (9). These difficulties have been 
chiefly attributed to the fact that the earth surface is 
not quite homogeneous and the heat conduction char- 
acteristics and surface temperature cycle of the earth 
surface region vary with respect to physical condition 
and time. The previous difficulties are, however, 
mostly due to the misuse of formula (4). It is not 
correct to use formula (4) either for analysis of one 
given periodic cycle of observation, or for analysis 
of depth variation of thermal diffusitivity. Instead, 
formula (4) should be used for the entire temperature 
data, covering the entire depth and the entire time 
period. For instance, if ten years of monthly average 
earth temperature cycles were complete with five 
depth readings, all of the 600 data points should be 
considered simultaneously for formula (4) to yield 
the best fit values of t,, (ai, i=1, 2,3, .. .) and (bj, j=1, 2, 
3...) as well as the best fit value of a. 

It is important here to mention that the foregoing 
statements are based upon the following statistical 
hypothesis: 

Time data (including all the cyclic periods) follow 
basically a given harmonic function of a period T 
(8766 hr for an annual earth cycle). Any deviation 
of actual cycles from a prescribed harmonic function 
are statistical rather than functional. 


4. The Least Squares Technique 


Although the technique described in this paper 
applies to any kind of periodic temperature data ob- 
served at several depths of any semi-solid system, 
the specific treatment will be given for annual cycles 
comprising monthly average earth temperatures taken 
at selected depths. 

Assume that the number of annual cycles of earth 
temperature is M, and each annual cycle consists of 
twelve monthly average readings at each of M,; depths. 
If there are no missing data, although this is not a 
requirement, the total number of observations M 


should be 


M=M, X M; X 12. 


Equation (4), however, should be modified to include 
only the finite number of harmonic terms up to some 
larger number N. In other words, it is assumed that 
a; and }; will be numerically insignificant for all i > N. 
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And SS= > 


The best unbiased estimates of ao, ai, .. . ay, 
bi, bg . . . by, and D, that of @ are given by the solu- 
tions of the following equations 


aSS _ 


“a 


for i=1,2. 





dSS 
6= aD & (11—a) 


where Ao, A:, 42, . . . An, Bi, Bo, . . . By and D are 
the best estimates of do. a;. a2, . . . Gv, bo. b1, bo, . . - 
by, and a. 
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. “ae L =u 
Si, k= e7 Bm sin( T 6: ~Bixs) 


and by noting that 

aSS : _ ISS 

—=() is also satisfied by —"=0 as long as D¥ 0, 
aD 0p 


i 


the equations in (11) can be expressed simply by the 
linear matrix equation shown in table 1. The left- 
hand side of the equal sign of table 1 is actually the 
product of two matrices, one square matrix of (2N 
+1)x(2N+1) and one column matrix of dimension 
(2N+1). Equation (1l-a) is nonlinear and is ex- 
pressed in terms of linearized variables, C;,; and Sj, x 
as is shown in table 2. 

Determination of Ao, A. Ao, ° ° Ay, B,, B:, i etn 
By as a best estimate of ao, ai, d2,. . . av, bi, be, . 
by is straightforward if the C;,;,’s and S;,,’s are known. 
However, these variables include an unknown D which 
is to be solved by (11—a). 

The method employed in this paper is to calculate 
a number of sets of Ao, Ai, Bi; i= 1, 2,. . . N, for sev- 
eral assumed D’s. These calculated sets are then 
applied to the equation in table 2. By an iterative 
procedure, a set of Ao, Ai, Bi; i=1, 2,. . . N, and D 





TABLE 1. Matrix expression of least squares equations (normal equations) 











Least Squares Equations 
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TABLE 2. Mathematical formula for the partial « 


lerivative of the sum of squares of all of the differences between observed and expected earth 
temperatures with respect to thermal diffusivity 
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5=0 


can be found to satisfy 5=0 (see table 2). These co- 
efficients and the assumed D are then substituted 
into formula (12) for the calculation of the mean 
squared deviation of the fitted least square values from 
the observed data by the following equation: 


The value of A, or the square root of the average of 
the sum of squares of deviations is called hereafter 
the root mean squared deviations (RMSD) of the ob- 
served temperatures from those calculated by eq (4). 
The term “standard deviation” is purposely avoided 
mr here since the properties of these deviations resulting 
55 from the iterated nonlinear fitting procedure are not 


A?=—- > 
V (15 known. 
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TABLE 3. Observed earth temperature data of Lemont, Ill., by 
r Carson 


Depth (cm) 
Month a 
100 





4.0 
3.0 
3.5 
6.0 
10.2 


15.5 
18.7 
19.8 
18.8 
15.9 


12.1 
7.8 
4.6 
3.5 


3.9 


6.7 
10.6 
15.4 
19.0 
19.7 


18.9 























Sample Calculations: The normal equations of 
tables 1 and 2 have been solved for annual earth tem- 
perature data of Lexington, Ky., and Lemont, IIl., 
which had been obtained by Penrod [1960] and Carson 
[1963], respectively. The Kentucky data covered five 
years, from 1952 through 1956, at six depths, 0, 61, 
122, 183, 244; and 305 cm. Except for the data at 61 
em and 122 cm during July and August of 1953, the 
data are complete. Lemont data shown in table 3 are 
complete for the period of three years, from 1953 
through 1955, at seven depths, 1, 10, 20, 50, 100, 305. 
and 884 cm, except at the depth of 884 cm during May. 
June, and July of 1953. 

Figures 1 and 2 illustrate the trend of 6 (table 2) 
and corresponding values of A calculated for several 
values of thermal diffusivities for the earth temperature 
data of Carson and Penrod, respectively. The calcu- 
lations in figures 1 and 2 were carried out for V= 4. 

It is seen in figure 2 that the value of A, the RMSD 
of calculated temperatures from the observed data, 
decreases rapidly to the value of 1.07 °C as the assumed 
thermal diffusivity increases to a value at 0.005 
cm?/sec. .On the other hand it is seen from figure 2 
that the diffusivity values that satisfy 5=0 for the 
Lemont data is a=0.0058 cm2/sec. Also it will be 
interesting to compare these thermal diffusivities with 
those calculated by conventional techniques by means 
of amplitude decay and phase angle shift by Penrod 
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from his Lexington data and with the data deduced 
from Carson’s data for Lemont, IIl. 


Lexington data a,; = 0.0088 cm?/sec 


ap: = 0.0056 cm2/sec 


Lemont data ag,;=0.0048 cm2/sec 


ap:;=0.0045 cm2/sec. 


Although the thermal diffusivities calculated by con- 
ventional amplitude and phase angle techniques differ 
from one another considerably, the values of A cor- 
responding to their diffusivities are all less than 1.3 °C 
in their ranges, as seen in figures 1 and 2. (This small 
change of A with respect to a wide range of thermal 
diffusivity may be the very reason that the two methods 
can yield quite different thermal diffusivities.) 

The mean squared deviation A? may not be as in- 
sensitive to the thermal diffusivity for the other types 
of periodic temperature data as it has been found to 
be in the analysis on annual cycles of monthly average 
earth temperatures. 

The Lexington earth temperature data were used 
to examine the influence of including higher harmonic 
terms on the values of A and it was found that the 
least squares analysis, including the higher harmonic 
terms up to degrees 2, 4, and 8, do not yield values 
of A which are significantly different from the results 
of the analysis ignoring all of the higher harmonics 
for all the thermal diffusivities assumed. 

Sample sets of least squares constants calculated 
for these two stations are shown in table 4. 


TABLE 4. Least squares constants for earth temperature data of 


Lexington, Ky., and of Lemont, Ill. 


n, Ky emont, Ill 


4 249 
} } 
0.0050 0.0058 
1.05 72 

14.48 


10.58 


11.39 
10.63 
0.08 
16 
x0 





It may be interesting to compare the calculated earth 
temperatures based upon the least squares constants 
with the original observed data. Table 5 shows the 
calculated earth temperatures when V=4 is employed 
for formulas of table 1 on the observed data of Lemont, 
Ill. (table 3). Similar calculations such as those for 
table 4 for N=1 were also made to evaluate the im- 
provement of taking higher harmonics over using the 


simple harmonic expression for the annual cycle of 


Lemont data. According to such comparison changes 
in the calculated results due to the inclusion of higher 


harmonic terms are very small. Table 5 also may be 
compared with table 6, which is the arithmetic average 
(or norm) of the 3-year cyclic data of table 3. Sig- 
nificant discrepancy of the arithmetic average from the 
least squares result exists for the region near the 
surface. 


TABLE 5. Predicted earth temperatures for Lemont, Illinois, using 


least squares constants up to fourth harmonics 


Depth (cm) 
Month }-— —————— —- 
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earth temperature data in table 3 
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Figure 3 compares the least squares calculation of 
the first harmonic constants with those derived from 
the observed three annual cycles of earth temperature 
in Lemont, Ill. 

The coordinates a; and 6b, of the figure represent 
first harmonic constants describing the observed earth 
temperature at a depth x from the ground surface. 
Since the Lemont data consists of seven depth observa- 
tions for a period of three years, each depth point is 
indicated with three numerals representing the first, 
second, and third year. 

The solid curve was obtained by connecting a; and 
by calculated by theoretical relation (depicted in the 
upper right corner of fig. 3) for several depths. The 
calculations of a; and b, were made with the least 
squares constants indicated in table 5. The curve 
representing the computed values of az and by fit 
extremely well for all of seven groups of three points. 
A similar comparison was made for the second har- 
monic constants of the same earth temperature data. 
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Graphical representation of harmonic constants for the earth temperature cycles 


of Lemont, Ill. 


Although the scatter of these second harmonic con- 
stants derived from observed data is much more than 
those indicated in figure 3, the computed curve did 
closely follow the centroid of the three points. 


5. Conclusions 


With the use of annual cyclic earth temperature 
data of Lexington, Ky., and of Lemont, Ill., a new 
technique is illustrated to yield single thermal diffu- 
sivity that will satisfy the least squares relation from 
the entirety of observed data. The technique requires 
an iterative solution of a nonlinear part of normal 
equations. The calculated earth temperatures based 
upon the theoretical heat conduction relation agrees 
very well with the average annual cyclic data of the 
monthly average earth temperatures when the least 
squares constants are properly evaluated. 

It has also been found in this analysis that the devia- 
tion of the observed earth temperatures from the 
calculated temperatures is fairly insensitive to the 
thermal diffusivity as long as the thermal diffusivity 
is in the neighborhood of 0.005 cm?2/sec, at least for 
the annual cycles of monthly average earth tempera- 
tures for two stations studied in this paper. This 
insensitivity may be the very reason that the current 
and existing techniques resulted in two quite different 
thermal diffusivities depending upon the method of 
calculation used, the amplitude method or the phase 
angle method. 


The method developed in this paper is recommended 
as the preferred method for the determination of 
thermal diffusivity from the periodic data of depth 
temperatures, not only of the earth surface by of other 
celestial bodies as well, in the event such data becomes 
available. 


This paper is a modification of a report prepared 
under the contracted research for the Office of Civil 
Defense. The permission from OCD to release the 
paper for publication is appreciated. 

The author greatly appreciates much critical review 
given by A. McNish of the National Bureau of 
Standards. 
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A New Near-Zone Electric-Field-Strength Meter 
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The National Bureau of Standards has recently completed the development of prototype instrumen- 
tation for measuring the electric-field components of complex, high-level, near-zone electromagnetic 


fields. 


The instrumentation is intended for use in evaluating hazards of high-level electromagnetic 
radiation to electroexplosive ordnance devices at military installations. 


The measuring range is from 


0.1 to 1000 V per meter, at frequencies from 150 kHz to 30 MHz, with a present uncertainty of less than 


+2 dB. 


The design of the NBS meters is based on the use of a novel form of telemetry, which apparently 


has not been fully exploited heretofore. 


This involves the use of a completely nonmetallic electrical 


transmission line over which the field information is transmitted from the measuring antenna to a 


remote readout unit. 


The line is essentially “transparent” to the field being measured, and reduces 


the perturbation of the field two orders of magnitude below that normally experienced when using a 


metallic line. 


The high r-f line loss involved necessitates miniaturizing the r-f portions of the receiving 


and calibrating instrumentation and placing them and their associated battery supplies inside the 


measuring antenna. 


The design and performance of the meters are discussed in some detail. 


Key Words: Device, electroexplosive ordnance; field, near-zone electromagnetic; hazards, elec- 
tromagnetic radiation; line, nonmetallic electrical transmission; line, semiconducting 
plastic transmission; meter, electric field-strength; telemetry, novel form of. 


1. Introduction 


The National Bureau of Standards, under the spon- 
sorship of the Defense Atomic Support Agency, has 
recently completed the development of prototype 
instrumentation for measuring the electric-field com- 
ponents of complex, high-level, near-zone electro- 
magnetic fields. The instrumentation is capable of 
measuring both the magnitude and direction of ellip- 
tically polarized CW electric fields having strengths in 
the range from 0.1 to 1000 V per meter at frequencies 
from 150 kHz to 30 MHz. This represents one step 
in the effort to develop accurate and meaningful 
measuring and calibrating instrumentation, for use 
in evaluating hazards of electromagnetic radiation to 
electroexplosive ordnance devices. It is well known 
that high-level fields of powerful nearby radio trans- 
mitters may cause premature detonation of missile and 
rocket type weapons during storage or loading opera- 
tions on shipboard, or at other military installations. 
It was from the standpoint of improving both the 
weapons reliability, as well as the safety of operating 
personnel and equipment, that the program was 
undertaken. 

The design of the NBS meters is based on the use 
of a novel form of telemetry, employing a completely 
nonmetallic electrical transmission line to avoid per- 


*For complete design details the reader is referred to “A New Near-Zone Electric-Field 


Strength Meter” — Frank M. Greene, National Bureau of Standards Technical Note No. 345, 


November 15, 1966. For sale by the Superintendent of Documents, U.S. Government 
Printing Office, Washington, D.C. 20402. Price 35 cents. 

**Radio Standards Engineering Division, National Bureau of Standards, Boulder, Colo. 
80302. 
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turbing the field being measured. This approach has 
apparently not been fully exploited, heretofore. The 
extremely high r-f loss of the line used attenuates any 
r-f currents induced on the line by the surrounding 
field and essentially eliminates reradiation from the 
line, or unwanted coupling between the line and the 
measuring antenna. The high line-loss, however, 
necessitates miniaturizing the r-f portions of the re- 
ceiving and calibrating instrumentation and placing 
them and the associated battery supplies inside the 
measuring antenna. The field information contained 
in the detected DC-AF output of the receiver is trans- 
mitted over the line to a remote readout unit, where 
the strength of the electric-field component parallel 
to the axis of the antenna is read directly. 

Interim electric-field standards and calibration tech- 
niques were also develped at NBS to evaluate the 
performance of the field-strength meters during their 
development, as well as to provide a tentative calibra- 
tion of the completed instruments. The present 
uncertainty of the standards is believed to be less than 
+2 dB, but further development effort is expected to 
reduce this to less than +1 dB. 


2. The New “Semiconducting” Plastic 
Transmission Line 


2.1. Errors Caused by Metallic Lines 


In the past, electric field-strength meters have’ 
usually made use of a long metallic r-f transmission 





line to connect the measuring antenna with its receiver, 
usually located at a point remote from the antenna. 
Such metallic lines often cause large measurement 
errors; especially when measuring near-zone fields 
having complex spatial distributions. in these cases, 
a neutral-field path, along which the line can be placed, 
either seldom exists, or the orientation of such a path, 
if it does exist is usually not known in advance of 
making the measurements. Under these conditions, 
it is difficult or impossible to orient the line so that it 
is everywhere normal to existing electric-field compo- 
nents. Thus, the line not only may perturb the field 
being measured, but unwanted r-f currents induced 
on the line can be coupled into the antenna and con- 
tribute to the total response of the field-strength meter. 


2.2. The Use of a Nonmetallic Line 


The NBS near-zone meters avoid this difficulty by 
making use of a special nonmetallic balanced trans- 
mission line in which the conductor r-f loss has been 
purposely made extremely high compared to that of 
the usual copper line. If the conductors are made of 
sufficiently high resistance material, the line can be 
made essentially “transparent” to the surrounding 
field. This is to say that the r-f currents induced on 
the transmission line by the field will be negligibly 
small, resulting in insignificant reradiation, or that any 
r-f energy propagating along the line will be heavily 
attenuated because of the extremely high loss of the 
line. It has been found that this can be achieved if 
the volume resistivity of the conductor material used 
is of the order of a million times or more higher than 
that of copper. This is roughly midway (on a loga- 
rithmic scale) between the volume resistivity of metals, 
on the one hand, and that of insulators such as glass or 
mica on the other. Such materials can therefore be 
said to be “semiconducting.” 


2.3. Description of the “Semiconducting” Line 


The material used in the NBS lines is basically 
polytetrafluoroethylene (PTFE), rendered “‘semicon- 
ducting” by uniformly dispersing finely divided carbon 
black (approximately 30 percent by weight) throughout 
the plastic while it is in the semifluid state during its 
manufacture. Tests at NBS have indicated the 
PTFE material to be more: stable electrically than 
other types tried, such as silicone rubber, polyethylene, 
carbon-coated multifiber glass, etc. Other types of 
conductor materials (including high-resistance metallic 
alloy wire and thin deposited metallic films) were 
examined, but were not found to be practical for this 
use, either because of unsuitable electrical charac- 
teristics or excessive cost. 

Several parallel-conductor balanced transmission 
lines were made at NBS using the semiconducting 
PTFE material. The line conductors are in the form 
of 0.03 in diam monofilaments in place of the usual 
copper conductors. They are coated with a thin nylon 


film (approximately 0.005 in thick) to further improve 
their strength and electrical stability. Each of the 
filaments is in turn enclosed in a vinyl-coated woven 
fiber glass sleeve approximately 0.1 inO.D. Twosuch 
insulated conductors are encased in a_ polyvinyl 
chloride outer jacket, for added strength and protec- 
tion, to form the completed transmission line. Military 
Type JJ-048 or PJ-291 twin-conductor plugs are used 
at the two ends for connecting to the dipole antenna, 
and to the remote indicator unit, respectively. The 
overall length of the lines is approximately 30 ft. 


2.4. Electrical Characteristics of the Line 


The semiconducting PTFE material used in the line 
has a volume resistivity of approximately 3.0 O-cm 
(compared to 1.7 10-® Q-cm for copper), giving a 
resistance of approximately 20,000 © per lineal foot, 
or a loop resistance of approximately 1.2 MQ for a line 
30 ft in length. The measured mutual capacitance 
between the two line conductors is approximately 
10 pF per foot, or roughly 300 pF for the entire line. 
The transmission loss, to the desired field-strength 
information at dc and low af, is negligible if a high- 
impedance d-c a-f vacuum-tube voltmeter is used 
as the readout device at the receiving end of 
the line. The line attenuation increases as_ the 
square-root of. frequency, reaching a theoretical as 
well as a measured value of approximately 53 dB per 
foot at a frequency of 30 MHz, for the differential mode 
of transmission. The common-mode attenuation with 
the line lying flat on a metal ground plane turns out to 
have essentially the same value. The use of a high- 
loss, parallel-conductor, balanced transmission line of 
this type has been found to reduce the perturbation 
of the surrounding field by more than two orders of 
magnitude below that existing in the case of a copper 
line. This renders any effect of the line on the field 
wholly negligible in most instances. A type of “‘flex- 
ural’”” noise is generated on the transmission line 
when it is suddenly flexed or moved physically. The 
common-mode component of this “noise” may be 
relatively high, having a peak value of several volts 
on open circuit. However, the differential-mode 
component is appreciably smaller and does not present 
a problem in the NBS meters at the level of d-c 
transmission used (see table 1). 


3. Description of the Near-Zone Field-Strength 
Meters 


3.1. Principles of Operation 


Two prototype near-zone field-strength meters were 
developed for the Defense Atomic Support Agency, and 
cover the frequency ranges 150 to 250 kHz, and 18 to 
30 MHz, respectively. These meters make use of 
a tuned r-f type of receiver to amplify the voltage in- 
duced in the receiving antenna by the component of 
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the electric field being measured. A tunable CW 
oscillator is included to permit calibrating the overall 
receiver gain periodically during a sequence of meas- 
urements in the field. A block diagram of the basic 
circuitry is shown in figure 1. The meters employ 
the same basic method of receiver-gain calibration 
that has been used in CW-type field-strength meters 
in this country for many years. The receiver portion 
of the field-strength meter consists basically of an r-f 
input capacitive step attenuator, a fixed-tuned band- 
pass r-f amplifier (with manual gain control), and a 
diode detector. The r-f output level of the oscillator 
is monitored by its own detector. When calibrating, 
the receiver gain is simply adjusted until the d-c out- 
puts of the two detectors are identical. In this 
method, the exact level of the calibrating signal is not 
important, and does not have to be known. However, 
it is important that the oscillator monitoring detector 
remain stable. 


TABLE 1. Distribution of nominal operating rf-voltage levels 


(in the TRF field-strength meter) 


E(V/m) | adB) Vi volts) Volvolts) V (volts) V.(volts) V,(volts) 


0.1-1.0 0 
1.0-10 20 
10.—10? 40 
102-108 60 


0.015-0.15 
0.15-1.50 
1.50-15.0 
15.0-150. 


0.003-0.03 0.003-0.03 
0.03-0.30 0.003-0.03 
0.30-3.00 0.003-0.03 


3.00-30.0 | 0.003-0.03 


0.3-3.0 
0.3-3.0 
0.3-3.0 
0.3-3.0 


0.75-7.5 
0.75-7.5 
0.75-7.5 


0.75-7.5 


E(V/m) = Electric field strength, volts per meter 
a(dB) Measuring dipole attenuator setting 

Vi Induced voltage in measuring dipole 

‘0 Output voltage of measuring dipole 

4 Input voltage to rf amplifier 

A Output voltage of rf amplifier 

3 Output voltage of amplifier detector (DC). 


y 
U 
U 
U 


A list of the various ranges of r-f voltage over which 
the principal components operate is shown in table 1 
for the range of electric field strength from 0.1 to 1000 
V per meter. An explanation of how the ranges were 
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Block diagram of the NBS field-strength meter. 


calculated is given in the following: 
Effective length of measuring dipole..., Leg = 0.15 m 
Induced voltage in dipole Vi=E + Lets 
Dipole voltage-transfer ratio, Vo/Vi=Cq/(Cat+C,)=0.20 
C, = Internal capacitance of dipole Ca =3 pF 
C, = Residual gap and circuit capacitance, 
C, = 12 pF 
, V2/V; = 100. 


RF amplifier nominal voltage gain 


3.2. Electrical Design 


The Dipole Antenna. A cutaway view of the dipole 
assembly is shown in figure 2, and a view of an assem- 
bled unit is shown in figure 3. The dipole is made of 
1g in wall copper tubing 1°/s in O.D. The length is 
approximately 13 in with a ¥/s2 in insulated gap at the 
center. The miniaturized solid-state receiver and its 
battery supply are contained in the left half, and the 
CW calibrating oscillator and its battery supply are 
in the right half. The capacitive step attenuator, the 
DPDT line-transfer switch, and the oscillator ON-OFF 
control comprise the five rotary switch wafers seen in 
the central portion of the dipole assembly. The 
balanced transmission-line plug connects to two 
insulated stainless steel pins at the center of the dipole. 
The line enters the interior of the dipole through a two- 
section balanced RC filter to help preserve the elec- 
trical symmetry or balance of the dipole with respect 
to the transmission line, and to attenuate any differen- 
tial, as well as common-mode r-f pickup that might 
exist on the line. The symmetry is important if the 
antenna is to respond only to the component of the 
electric field parallel to its axis. The symmetry of 
the NBS dipole is such that its response to the cross- 
polarized component of the electric field is 40 dB or 
more below the principal response. 

It can be shown that if the length of the measuring 
dipole does not exceed 0.03 wavelength, its effective 
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FIGURE 2. Cutaway drawing of the measuring dipole. 











FIGURE 3. View of the complete measuring dipole. 
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length when measuring a field having a highly complex 
spatial distribution will not differ by more than 2 or 3 
percent (in the worst case) from the effective length 
when immersed in a uniform plane wave. Likewise, 
the interaction error, resulting from coupling between 
such a short antenna and its image, will be small pro- 
vided the loading is light and that measurements are 
not made closer than 2 or 3 ft to the ground or large 
metallic objects. 

RF Capacitive Step Attenuator. This attenuator has 
four principal steps, 0, 20, 40, and 60 dB, with two 
additional 0 dB positions, calibrate A(OSC) and cali- 
brate B(AMP).._ In the latter two positions, the calibrat- 
ing oscillator is automatically turned ON, and the 
“semiconducting” transmission line can then be 
switched between the d-c outputs of the oscillator and 
amplifier detectors, while calibrating the amplifier 
gain. The attenuator is constructed using special 
subminiature rotary wafer switches. Miniature low- 
temperature coefficient ceramic capacitors are used 
in either an “L” or a “Pi” configuration for the various 
attenuator pads. The attenuator steps are accurate 
to +0.1 dB. The attenuator module is approximately 
1 in O.D. by 1 in long. 


Band-Pass Amplifier and Detector. In order to 
provide a high input impedance, so as not to appre- 
ciably load the dipole antenna or the capacitive step 
attenuator, a metal-oxide-silicon field-effect transistor 
(MOS-FET) is used in the first stage as an untuned 
input-impedance converter. This is followed by a 
six-pole, fixed-tuned, bandpass filter having a 3 dB 
pass-band of 150 to 250 kHz, and 18 to 30 MHz respec- 
tively in the two units and a 40 dB pass-band of 100 to 
300 kHz and 12 to 36 MHz, respectively. The filter 
is followed by a 3-stage r-f amplifier. In the low- 
frequency unit, this is broadband resistance-coupled, 
using overall negative r-f feedback to improve gain 
stability with respect to temperature and _ supply- 
voltage variations. In the high-frequency unit, the 
r-f amplifier has 3 stages using wide-band synchro- 
nously tuned interstage transformers. The overall 
gain of both amplifiers can be manually adjusted to the 
required operating value of approximately 40 dB. 
The output detectors of both amplifiers employ two 
stable silicon diodes in a voltage-doubling circuit, to 
ease voltage-swing requirements on the final stage and 
consequently to reduce battery drain. A 12 V silver- 
oxide battery supply is used, having a useful life of 
from 10 to 20 hr. The complete amplifier module 
measures 11/4 in O.D. by 5 in long including its filter 
and battery supply. 


Calibrating Oscillator and Detector. This unit is 
also transistorized. It consists of a permeability-tuned 
Clapp-Colpitts oscillator followed by a two-stage auto- 
matic-level-control amplifier to maintain the r-f output 
essentially constant over its tuning range. The tuning 
range of each of the two oscillators is the same as the 
3 dB bandwidth of the band-pass amplifier with which 
itis used. The r-f output of the oscillator is adjustable 
to roughly 3.0 V and fed into a fixed 30 dB resistive 
attenuator pad. When calibrating the gain of the 
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band-pass amplifier, the output of approximately 0.1 V 
from the pad is applied to the amplifier input through a 
small variable capacitor (adjusted to approximately 
3 pF) which simulates the dipole probe capacitance. 
The detector circuit is identical to that used with the 
band-pass amplifiers and is used to monitor the 3.0 V 
oscillator level at the input terminals of the 30 dB pad. 
A 13.5 V silver-oxide battery supply is used, having a 
useful life of about 20 hr. A portion of the supply is 
Zener regulated and applied to the transistor base 
circuits for improved amplitude stability. The com- 
plete calibrating oscillator module measures 11/4 in 
O.D. by 5 in long including the battery supply. 

The Remote Indicator Unit. This unit consists 
basically of a battery-operated, balanced d-c vacuum- 
tube voltmeter, with its d-c microammeter, on which 
the magnitude of the electric field strength is indicated. 
This VTVM has exceptionally good stability and 
linearity of response because of the large amount of 
negative d-c feedback employed in its design. It is 
used over a single 20 dB input range of approximately 
0.75 to 7.5 V. The balanced input circuitry provides 
from 40 to 60 dB rejection of common mode input in the 
60 to 400 Hz power-frequency range. In addition, a 
two-section balanced r-f filter provides at least 40 dB 
suppression of common and differential-mode input at 
the higher r-f frequencies. An output circuit in the 
VTVM is used to drive a transistorized audio amplifier 
and loud speaker. These units are used as a zero-beat 
indicator when adjusting the frequency of the calibrat- 
ing oscillator to that of the field being measured. 
Separate, self-contained battery supplies are used 
with the VTVM and a-f amplifier. The complete 
remote indicator is housed in a drawn-aluminum case, 
using a tight-fitting panel to minimize r-f leakage into 
the unit. A small drawn-aluminum case is used to 
enclose the rear side of the d-c microammeter to reduce 
r-f entry through the meter face. 

The Assembled Field-Strength Meter. Tine meas- 
uring dipole of the near-zone field-strength meter is 
mounted on a hollow fiber-glass shaft which is perpen- 
dicular to its axis. This shaft is in turn supported in 
plastic bearings in a circular reinforced polyfoam 
gimbal, approximately 21 in. in diameter, as shown in 
figure 4. This assembly is mounted on a sturdy 
wooden surveyor’s tripod, in such a manner that 
manual adjustments in both azimuth and elevation 
can be readily made in the orientation of the dipole. 
A view of the complete field-strength meter is shown 
in figure 5. The balanced transmission line connects 
to the dipole through the hollow shaft on the right. 
The other end of the 30 ft line connects to the remote 
indicator unit, a closeup view of which is shown in 


figure 6. 


3.3. Performance 


The two near-zone field-strength meters underwent 
extensive field tests at the U.S. Naval Weapons Labo- 


ratory at Dahlgren, Va. Electric field strengths were 
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FIGURE 4. View of the measuring dipole and mounting gimbal. FIGURE 6. View of the remote indicator unit. 


measured without difficulty over the entire range from 
0.1 to 1000 V per meter (CW-RMS). There was no 
observable susceptibility to out-of-band CW signals 
having levels as high as 200 V per meter, or to pulsed 
signals from nearby high-powered radar transmitters. 
No interaction was observed between the semicon- 
ducting transmission line and the field being measured. 
In general, the agreement was about + 1 dB between 
the instrument readings as calibrated in the NBS 
standard field at Boulder and those obtained in the 
NWL standard field at Dahlgren. The uncertainty 
in the calibration accuracy of the instruments is 
believed to be less than+2 dB. Further development 
effort on the standard fields is expected to reduce this 
to less than + 1 dB. 





The work described in this paper was performed for 
the Field Command, Defense Atomic Support Agency, 
under DASA Order No. HD-1102-3171—0001 (June 
20, 1963). 
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assistance of William E. Jessen of NBS-BL in the 
assembly and performance evaluation of the instru- 
mentation described herein. The assistance of Ezra 
B. Larsen in the standard-field development work is 

FIGURE 5. The complete near-zone field-strength meter. also gratefully acknowledged. 
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A general method is developed for determining the inductance and characteristic impedance of 


uniform transmission lines. 


A non-uniform current distribution is allowed in the transverse plane. 
The system is represented by a matrix equation which can be programmed for computer solution. 
correct inductance and impedance are obtained as the result of a simple limiting process. 
is applied to one particular geometry, a four-tape stripline system. 


The 
The method 
Results are given for the induc- 


tance, resistance, and current distribution as functions of frequency and resistivity for a particular 


geometry. 


A method for extending the results to strip lines with proportional dimensions is developed. 


An accuracy of one part in 10° was found to be feasible for the determination of the inductance per 


unit length. 


Key Words: 


1. Introduction 


The determination of the inductance, capacitance, 
or characteristic impedance of coaxial systems has 
been the subject of numerous reports. The system 
considered here is the special case of a strip- 
transmission line consisting of four vanishingly thin 
parallel tapes. 

The approach used in this paper is to determine the 
current distribution in a transverse plane and therefore 
the transverse magnetic field configuration. Knowing 
the current distribution, the effective inductance per 
unit length is found directly. The characteristic 
impedance of a lossless line can then be determined 
from the inductance. 

Almost all previous authors have used the approach 
of calculating the capacitance per unit length, C, after 
determining the configuration of the transverse electric 
field. The result is then used to calculate the charac- 
teristic impedance of the lossless line, using the 
relationship 


Zo=1fvC (1) 


where v is the velocity of propagation in the line. 
Several applications of this approach are listed in the 
references [1—4].!' The results of Cohn [2] and Bates 
[3] give analytic results which are useful design tools. 

The approach used in this paper, while nonanalytic, 
is quite general and without any geometric limitation 

* Microwave Circuits Standards, Radio Standards Engineering, NBS Laboratories, 


Boulder, Colo. 80302. 


' Figures in brackets indicate the literature references at the end of this paper. 


Characteristic impedance, inductance, and strip-transmission lines. 


in the transverse plane. The effect of finite conductor 
losses and frequency dependence can also be included. 
It is not necessary to assume any distribution or magni- 
tude for the interdependent variables as this approach 
determines all necessary information simultaneously. 


2. Consideration of Finite Losses 


If the line has finite losses, the characteristic im- 
pedance and propagation constant, y, can be calculated 
from 


, _ |R+joL 
Zo= G+joC 





y=at+jB= V(R+ joL (G+ jal) (2b) 
where R and G are the series resistance and shunt 
conductance per unit length of the line and @ and B are 
the attenuation and phase constant of the line. For a 
low-loss line, (2a) and (2b) can be closely approximated 


by 
: ae ] 
— olf oe 


RR. GuL 
aoc = (3b) 


(3a) 


B=o ViC=— (3c) 





where Q; =@L/R and Q.=@C/G. Terms involving 


Q? have been neglected. It is also assumed that 


Sn R, 
=a, 8, and ._ 


which must be satisfied for transmission line theory 
to be valid [5] where o, and €, are the conductivity and 
dielectric constant of the conductor, R, is the skin- 
effect surface resistivity of the conductor, and 7 is the 
intrinsic impedance of the dielectric. 

Both R and L are functions of frequency and resis- 
tivity as well as geometry for lines with finite loss. 
The inductance is therefore a function of resistance 
so that R must also be calculated to obtain accurate 
values of L. A computer method for calculating R 
and L as a function of frequency, resistivity, and 
geometry is developed in this paper. Values of v and 
G are readily obtained from properties of the dielectric 
material in the line. These values together with R and 
L can be used to obtain the attenuation constant from 
(3b). If the contribution of Q, and Q, in (3a) can be 
ignored, the characteristic impedance can be obtained 
from the simple relationship 


Zo — vL. 


which is analogous to (1). 


3. D-C Inductance 


Figure 1 shows the cross section of a four-tape strip- 
line. The conductors marked 1 and 2 are to act as a 
single outer conductor while conductors 3 and 4 act as 
a single inner conductor in a go-and-return circuit. 
For conductors with finite resistance the current will 
be uniformly distributed throughout the outer and 
inner conductors at zero frequency. The total induc- 
tance of the system, , can be calculated from 


(4) 


14 4 u 
SFf=— ij 
4 = 2, e 


* as 


where Mj; is the mutual inductance between tape i and 
j, and Mj is the self inductance of tape i. A recent 








FIGURE 1. Cross section of a four-tape stripline. 


report by Hoer and Love [6] gives exact equations for 
the inductances of rectangular conductors carrying a 
uniform current. Using these equations to calculate 
the Mj, the exact d-c inductance of any rectangular 
coaxial system can be obtained. 


4. A-C Resistance and Inductance 


When an alternating voltage is applied to the sys- 
tem, the current is no longer uniformly distributed 
throughout the cross section of the conductors. The 
inductance as well as resistance then differs from 
the d-c values. 

Consider a short length of a long uniform line, such 
as shown in figure 2. Writing V2 in terms of V; gives 

V2=V, cosh yAl—1,Zo sinh yAl. (5) 
Choose Al such 
approximated by 


that |yAl| < Then (5) can be 


V2=V,—IZoyAl 


by. 
~ Al = IZ oy. 


Substituting (1) and (2) into (5) gives 


AV Vi-D, 
at Le + jwL). 


Alo d 
Our object now will be to calculate the total current, 
I, for some arbitrary value of AV/Al and from this 
result obtain R and L from (7). With R and L known, 
the transmission line constants may be obtained 
from (3). 

To determine the total current in the stripline, each 
tape is mathematically divided into a number of 
smaller, parallel sections as shown in figure 3. An 
equivalent circuit of the transmission line then looks 




















FIGURE 2. Short length of a long uniform transmission line. 











---— 


| 
ee 4nj8n,ee* ) | 7nt! 
tt | 





2ntl,2nt2, eee 3n] 7n, eee 





| 
| 
| 
| 
| 


| 
y 


FIGURE 3. Cross section of the stripline showing how the 


sections are labeled. 
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FIGURE 4. Equivalent circuit for the stripline with the tapes 


divided into small parallel sections. 


like that in figure 4. The width of each section is 
chosen small enough so that the current density may 
be considered uniform throughout that section. The 
appendix contains de inductance equations which 
can be used to calculate the self, and the mutual 
inductance per unit length of these sections. The 
resistance per unit length of each section will be the 
dc resistance. 

Although the method of subdivision is almost com- 
pletely arbitrary, one system will be explained, which, 
although strange in appearance, will simplify later 
expressions considerably. Each tape is divided into 
2n sections labeled in such a way that the index 
increases as the y axis is approached from any edge 
as shown in figure 3. The widths of the sections need 
not be equal and indeed will not be in later calculations. 
Since the geometry has symmetry about the origin and 
both axes, this symmetry is retained in subdividing the 
tapes. Symmetry produces the following relations 
between the currents in the different sections. For 
the outer tapes, 


T;=Ten + j= ans j=Tonsj, j= a ae 


For the inner tapes, 

Ins 5 =I 5n4j=15n+j=1in+j,j=1, 2, . - 4 f. 
The total voltage drop per unit length along the length 
of any one section may be written, 


8n 
Vi. = rel + jw bo Mili. (8) 
(=1 


For the purposes of this report the d-c resistance per 
unit length, r;, of each section will be used in the form 
of resistivity/area. There are 8n equations which may 
be written in matrix form as 


V; ry} 





M,,1 My sn ] I, 


(9) 


+ ja) 


Msn,1 . Man, sn I sn 





As a result of the current symmetry conditions, only 
2n of the equations are independent. The four groups 
k=1 to 2n, 2n+1 to 4n, 4n+ 1 to 6n, and 6n+ 1 to 8n 
are all equivalent. 

Choosing the first group, with =1 to 2n as the 2n 
independent equations to be solved, rewrite (9) as 





\ 
. M ,2n ' 


ote Mon,2n 





where 


M = y Mi,2pn + jo 


p=0 


i,j=1,2, . 


Equation (10) can now be written in the form 
V=[R+joM]. 
Since V;, and J from (8) are complex, let 


Vie=ext ifr, T= ay — jor. 





In terms of column matrices, 


V=E+jF, 1=A-jB. (13) 


Substituting (13) into (1]) produces two real matrix 
equations: 


E=RA+MB, 
F=—RB+oMA. 
Solving for A and B produces 


A=R-(1)+(@MR-')]-YEHOMR-)F] (16) 


B=R-"{(1)+(@MR-!P]-"[(oMR-YE-F]. (17) 


In these equations for A and B, E and F have not yet 
been determined. To simplify their determination 
the following matrix definitions are used: 


w= R-'{(1)+(@M@R -1)?]-1 
b=aMR-'! 


T=(]1,, . — ]5,) 


6=(0,, . ee | ee, ea . lon) (18) 


Since the total current in the outer conductor must 
equal the total current in the inner conductor, we have 


n 2 
Sa- ¥ 
k 1 k= 


Ss bi. =— 


k=1 k= 


or in matrix notation, using (16) 


in (18) 


and the definitions 


TA=TH[E+ oF] =0 (19) 


TB=7(hE —F | =0. (20) 


The scalar voltage drops in the inner sections and outer 
sections obey the following conditions: 


ot k=1,2,. 


—=€k 


OP ait. n+2,. 
=e,} 


Cy T Cn +1 +j(fi + fn + 


Cn+1 


1)= Av/Al. 


Since the value of Av/Al is arbitrary, it is set equal to 
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1+j0. This condition applied to the above results 
produces 


ei tensi1=1 


The column matrices E and F can now be written in 
terms of known matrices and one scalar unknown each: 


E=6" + er" 
(21) 
F=f" 


where the superscript T indicates the transposed 
matrix. Substituting (21) into (19) and (20) and solving 
for e; and f, produces 


(rhhr? \(rWhO? ) + (7WO? \(rWr") 
e = a 
(rr? )? + (rr? )? 





(22) 


- TUhO? \(rwr? ) — (rh? \(rWhT" ) 


fi (tusr? )2 + (rub? )? 





(23) 


Every term in parenthesis on the right sides of these 
equations is a calculable scalar. The result obtained 
can be substituted into (16) and (17) and the current 
in each section can be determined. With the currents 
known, the approximate inductance and resistance per 
unit length can be calculated from (7), which in terms 
of the components of the current in each 
becomes 


section 





and 





5. Limiting Value of L, and R,, 


The correct value of inductance and resistance per 
unit length is obtained in the limit as n becomes 
infinite. A good estimate of the correct L and R can 
be obtained from several approximations calculated 
for different values of n. The approximate values of 
L and R are plotted against n as shown in figure 5 and 6. 
The asymptotes to these curves are the values as n 
becomes infinite. 
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FIGURE 5. Approximate inductance versus number of subdivisions, 
showing limiting value of inductance for infinite n. 


The curve of L versus n can be approximated by an 
equation of the form 


Enp=Lot+an— JS 


(26) 


where L, is the approximate inductance obtained from 
(24) for a given n, Lo is the inductance as n becomes 
infinite, and a and.¥ are constants. Equation (26) can 
be solved for Lo if Ly is calculated for four different 
n’s chosen such that 

ny = Nes 


Ne n4 


Then from (26) 


(Ln —Ln (Ln — 
3 4 2 

(Ln ai 
1 


Ln ) 
Ln) — (Ln, — Ln) 





Lo=Ly — 


A plot of (Ly, — Loo) versus n for these 4 values of L, will 
be a straight line on log—log paper only if (26) is a 
valid representation of the curve L, versus n. Thus, 
four points can be used to verify (26) as well as to calcu- 
late Lo. Once it has been shown that (26) is a valid 
representation of the curve, Lo can be calculated from 
three values of n chosen such that 


ny Nz , 
—=— (29) 
Ne Neg 


Then (26) gives 


(Ln, — Ln)? 
(Ln, — Ln,) — (Ln, — Ln.) 





Lo= Ly — (30) 


The value of Loo for the curve shown in figure 5 was 
calculated from three sets of values of n. For values of 
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FIGURE 6. Approximate resistance versus number of subdivision, 
showing limiting value of resistance for infinite n. 


n equal to 4, 6, and 9, (30) gave a value of Lo=0.201352 
#H/m. For values of n equal to 9, 12, and 16, (30) 
gave a value of Lo=0.201344. For values of n equal 
12, 15, 16, and 20, (28) gave a value of Lo=0.201343. 
The maximum difference is 5 parts in 10°. The latter 
value of Loo was used to obtain the curve of (L, — Loo) 
versus n shown in figure 7. The time to run these 
sets of n on a large high speed computer was 10, 35, 
and 75 seconds respectively. 











20 
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FIGURE 7. Difference between approximate and limiting value of 


inductance versus number of subdivisions. 





No simple approximation was found for the curve 
of R versus n. 


6. Subdividing the Tapes 


The exponent, .7, in (26) is an indication of the rate 
of convergence of L, to its limiting value, Lx. Sub- 
dividing the tapes into segments of equal width gave a 
SY of approximately 1. The variation of the current 
density across the width of the inner and outer tapes 
is shown in figure 8. This curve suggests that the 
width of the segments of the inner tapes be made 
smaller near the edges where the curve is steepest. 
One such method @f subdivision is to let the width, 
Wkn+j of the kn+j segment be given by, 


C ‘£=1, 3, 5, 7. 
(n+2—j)" j=1, 2, 


(31) 


Wkn + j = 


where C is determined by setting the sum of the widths 
equal to W/2. Thus 





fet2-s» 


The constant, 7, was chosen such that the rate of 


convergence was near maximum. With the outer tape 
divided into segments of equal width and the inner 
tape divided into segments having widths calculated 
from (31), an 7 of 3 was sufficient to yield a ¥ of 
approximately 3. 


7. Typical Results 


If the line has losses, the R and L will be a function 
of the actual dimensions instead of just the ratios 
tlb, W/b and Wo/b. The R and L of a stripline having 
dimensions b= W=0.5 cm and t=0.025 cm, was cal- 
culated as a function of frequency f, resistivity p, and 
width of outer conductor Wo. The thickness, T, of all 
tapes was assumed to be equal and small enough to 
permit calculation of the Mj; in (10) using inductance 
equations for zero thickness tapes.” Machine storage 
limitations prevented including the variation of current 
density with thickness. This limits the strict applica- 
tion of these results to tapes whose thickness is of the 
order of a skin depth. For the geometry considered 
here the errors introduced even at high frequencies 
should be negligible 


2 See appendix. 
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FIGURE 8. Variation of current density in the inner and outer tapes. 
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FIGURE 9. Variation of inductance per unit length with q. 


The variation of inductance per unit length with 
frequency and resistivity is shown in figure 9. The 
reason for using the variable g=/Tk/p is that one curve 
of L versus q is sufficient for all f, 7, k, and p, provided 
that T meets the requirements just discussed. The 
scale factor k is a dimensionless proportionality con- 
stant for striplines having proportionate dimensions 
and shape. If all dimensions, except the thickness, 





are multiplied by a scale factor, k, the matrix M in 
(16) and (17) remains unchanged. The matrix R-! in 
(16) and (17) is a function of k and may be written 


(33) 


where 


W1 0 
UW = . 
0 W2n 


Equations (16) and (17) can now be written, 


- ae ae 
4=— 9 |(y+(2* ar) | 
p \ p 


[e+(“ar)F] 0 


a i (en MU | 


p 


| (A ar)e- FI. (35) 


Since the matrices .“# and WY are independent of f, T, 
k, and p, one curve of L versus the variable q=/Tk/p 
will be sufficient for all f, 7,4, andp. The same result 
is true for R/Rae versus q. For the strip line under 
consideration, with T=10-4 em, p=2X10-® Q em, 
and k=1, a q of 5X 104 hertz/Q corresponds to a 
frequency of 10* hertz. 

A curve of R/Rac for the stripline versus q is shown 
in figure 10. An explanation of the shape of the curve 
can be had from an examination of the current density 
curves such as shown in figures 8 and 11. The current 
density in the outer tapes remains uniform up to about 
qg=5. From q=5 to about g=5000 the current dis- 
tribution is changing. Above a q of about 5000 the 
current distribution in the outer tapes no longer 
changes. The resistance of the .outer tape would 
therefore remain constant up to a q of 5, increase until 
a q of 5000 is reached and then remain constant again. 
The current density in the inner tapes remains uniform 
up to a q of 50. The current then begins to crowd to 
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FIGURE 10. Variation of resistance with q. 
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FIGURE 11. Normalized current density in outer tapes 
as a function of q. 
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FIGURE 12. Variation of Resistance with q for the outer 
and inner tapes. 
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FIGURE 13. Variation of inductance per unit length with 


width of outer tapes. 


the edges of the inner tapes. The current continues 
to move to the edges as long as the value of g increases, 
never reaching a constant distribution as the current 
in the outer tape does. The curves of the change in 
resistance of the outer and the inner tapes with g would 
look like those shown in figure 12. The sum of these 
two curves gives a curve like that shown in figure 10. 

The variation of inductance per unit length with the 
width of the outer tapes is shown in figure 13 for 
different values of g, for the dimensions given above. 


8. Summary 


The inductance and resistance per unit length of 


strip transmission line are functions of the current 
distribution in the conductors. At frequencies where 
the current distribution throughout the conductors is 
uniform, the inductance and resistance can be calcu- 
lated exactly. For frequencies where the current 
distribution is not uniform, the conductors are divided 
into sections. From these sections the approximate 
current distribution of the conductors is calculated. 
Approximate values of the inductance and resistance 
per unit length are then calculated from the approxi- 
mate current distribution. Using high speed com- 
puters, it is not difficult to divide the conductors into 
sufficiently small sections to calculate the inductance 
to an accuracy of one part in 10° or better. The 
accuracy on resistance is considerably less, depending 
on frequency. Other methods of subdividing the 
tapes are being considered to improve the convergence 
of R, to its limiting value, R o. 

The variation of the inductance per unit length with 
frequency, resistivity, and width of the outer conduc- 
tors is given for a particular geometry. Tables and 
graphs of the inductance and resistance per unit length 
fer other geometries are being prepared for future 
publication. 

The technique of calculating R and L as given in this 
paper for a stripline can be applied to other lines having 
more complicated cross sections provided that the 


66 


conductors can be divided into sections whose d-c 
self and mutual inductances can be calculated. This 
technique is used in a forthcoming paper by Brooke 
and Cruz [7] to calculate the L and hence Zp of a loss- 
less, rectangular line with a rotating center conductor. 


The authors gratefully acknowledge the assistance 
and helpful criticism of C. M. Allred, and W. E. Little. 


9. Appendix 


The mutual inductance per unit length between two 
long, thin, parallel tapes such as shown in figure 14 
can be obtained from equation 8 of reference [6] and is 


E-—a, E+d 
(x) 
E+d-—a,E 


P2— x? = 
——— ln (P?+-x?2)—xP Tan-!— 


q 
1 P 


I 


+ 0.2| In 2/-+ |. MH/m (Al) 


where the limits, which have been retained for com- 
pactness, are substituted as follows: 


$183 4 


[f(x)] (x) = Ss (-~ U°* hie. 


as. 


The origin coincides with the left edge of one tape. 
If the left edge of the second tape is in any quadrant 
other than the first as shown, one or both of the values 
of E and P will be negative. The self-inductance per 
unit length of a long thin tape is 


oY es 
In 2/+<— 


|. H/m. 


] 
L,=0.2 In -+ 0.2| (A2) 
a 


In both of these expressions it is assumed that the 
current density is uniform throughout the conductors. 
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FIGURE 14. Cross sectional view of two long parallel thin tapes. 





Note that the term 


0.2] In 21+5| (A3) 


appears in both M, and L; as the only term involving 
length. However, if M, and L; are substituted into (8), 
Substi- 


the terms involving length will exactly cancel. 
tuting (Al) and (A2) into (8) gives 


sn 

r . — 
V;, = rich, +J@ > 
a 


(=1 


M1) + jel 0 (In 214 -)| Sh (A4) 
kui J | +s - : py l iM 


l=1 


where the M;, are the expressions in (Al) and (A2) 
with the (A3) term excluded. 
because 


But the last term is zero 


3 1, =0. 


t=1 


That is, the total current in the outer conductors is 
equal and opposite to the total current in the inner 
conductors. Therefore the ;, in (8) may be calculated 
from (Al) and (A2) with the (A3) term excluded. 
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The exact solution for the field equations of a cylindrical TMi;9 mode cavity has previously allowed 


accurate measurements of tensor permeability to be obtained at X-band frequencies. 
strated that this method is also applicable at frequencies down to 1 GHz. 
cavity and measurement system for obtaining data at these lower frequencies is given. 


It is demon- 
A brief description of the 
Both intrinsic 


and external permeability results on three commercially available polycrystalline garnets are shown. 
The larger size rods required for measurements at these frequencies result in some sample size effects 


in the data. 
materials was observed. 


In addition, a previously unreported absorption in the external tensor permeability of the 


Key Words: Degenerate mode cavity, L-band cavity, polycrystalline garnets, tensor permeability. 


1. Introduction 


A knowledge of the tensor permeability components 
as described by Polder [1]! of a ferrimagnetic material 
is often required in the theory and design of ferrite 
devices at microwave frequencies. Considerable effort 
has been directed towards the measurement of these 
parameters in recent years. Most of these tensor 
permeability studies have been carried out using cavity 
techniques [2-18]. In particular, degenerate mode 
cavities excited with either a positive or negative pure 
circularly polarized wave have been the most accurate 
means available for obtaining measurements of this 
type. Both TE and TM mode cavities have been used 
in this manner. In comparing the degenerate mode 
cavity methods, the measurement of rod-shaped 
samples placed along the axis of a TMiio mode cylin- 
drical cavity has been particularly attractive since this 
geometry allows an exact solution to be obtained for 
the field equations [19-20]. 

The measurement of tensor permeability using the 
TMi1o mode cavity as well as any of the other cavity 
methods has in general been confined to frequencies 
of 3 GHz or higher. Exceptions to this may be found 
in the work of Ogasawara [21—23] as well as Waldron 
and Maxwell [24-27], who use stripline techniques for 
measurements from 3 GHz down to the UHF range. 
Ogasawara describes measurements of disks in crossed 
stripline cavities. Waldron and Maxwell use thin 
rectangular plates in a single stripline cavity. In 


*Radio and Microwave Materials, Radio Standards Physics, NBS Laboratories, Boulder, 
Colo. 80302. 


' Figures in brackets indicate the literature references at the end of this paper. 


either case, the use of disk or rectangular shaped 
samples in stripline may be subject to errors due to 
sample geometry effects, nonuniform electromagnetic 
fields, etc. [28-29]. 

It appears feasible that some of the errors associated 
with sample measurements in stripline could be 
avoided, at least in the frequency range down to about 
1 GHz, by the use of the previously mentioned TMi 
mode cylindrical cavity at these lower frequencies. 
This cavity offers advantages over other waveguide 
cavities at lower frequencies for two reasons. First, 
this mode allows a flat cylindrical type cavity to be 
designed which is sufficiently thin along the axis, even 
at these lower microwave frequencies, to fit between 
the pole pieces of an electromagnet. Second, the 
rod-shaped samples and exact solution for the field 
equations mentioned above shouid ‘:elp to avoid 
sample shape errors appearing in other methods. 

This paper describes a feasibility study which uses a 
TMiio mode cavity for tensor permeability measure- 
ments in the L-band region. A brief description is 
given of the measurement method, the cavity, and the 
associated circuitry used in these measurements. 
Several improvements in the TMiio mode cavity are 
noted which are also applicable to the use of this 
method at higher frequencies. This is followed by a 
description of typical data obtained on three com- 
mercially available polycrystalline garnets at 1106 MHz. 


2. Measurement Technique 
2.1. General Method 


The general procedure for measuring tensor per- 
meability using degenerate mode cavities is based on 
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the fact that the magnetic flux b in a magnetic material 
is related to a circularly polarized field A; inside the 
material by the relation [4] 

b. =(wta)hi.. (1) 
In this equation, «x and @ are both complex and are the 
diagonal and off diagonal components of the intrinsic 
tensor permeability. The subscripts refer to either a 
positive or a negative circularly polarized wave. 

In a typical measurement, the frequency shift and 
change in Q of a cavity containing a small sample are 
determined as a function of an applied d-c field for 
each of the two senses of circular polarization. In 
most cases, the experimental relationship between 
(2+ a) and the changes in the resonant frequency and 
Q of the cavity are determined from cavity perturbation 
theory. However, in the case of the TMiio mode 
cylindrical cavity containing a rod-shaped sample along 
the axis, the correspondingly more exact relations may 
be obtained from a direct application of Maxwell’s 
equations to the sample and cavity geometry as shown 
by Bussey and Steinert [19-20]. This approach is 
well described in their papers and will not be repeated 
here. It suffices to point out that the intrinsic per- 
meability data in this paper have been calculated with 
the aid of a computer from eqs (10) and (11) of their 
paper [20].? 

To accurately identify the resonant field and the 
maximum value of the imaginary part of (u+ a), it was 
necessary to know the d-c field value at increments of 
one oersted immediately above and below resonance. 
The accuracy of the available equipment for measuring 
the d-c field in the range of values used in these meas- 
urements was not adequate to obtain this resolution. 
The relation between the measured frequency and the 
applied d-c field was very nearly linear. This made it 
possible to plot this relation and interpolate between 
experimentally observed points. The experimental 
data and the data obtained by this interpolation were 
combined to obtain the intrinsic permeability data 
shown in this paper. 

It should be emphasized that the permeability (4 + @) 


defined in eq (1) and calculated from the formulas of 


Bussey and Steinert is the intrinsic permeability which 
does not depend upon the geometry of the sample. 
This situation arises because the field A; in eq (1) refers 
to the rf field inside the sample. In some cases a 
so-called “external” permeability (7+) has been 
referred to in the literature. This is defined in terms 
of the rf field ho existing in the absence of the specimen. 
In this case, we have an analogy to eq (1) 


b.=(ntk)ho.. 


where 7 and « are again complex and are the diagonal 
and off diagonal components of the external permea- 


* The signs preceding the e” and €' terms of eq (11) in ref. 20 should be — and + respec 
tively rather than as given in the paper 
meability m* should have been m* =m 
in their discussion. 


This is because the expression for inverse per- 


-im" rather than m' + im" which the authors used 
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bility. In the case of a rod sample, Waldron [30] has 
shown that the external permeability is related to 
intrinsic permeability by the relation 


icinie 2(u+ a) 


~ 1+(uta) (3) 


In this paper, we are primarily concerned with the 
evaluation of the intrinsic permeability. However, 
the above concept of external permeability is included 
primarily because we have observed 


previously 
unreported absorption phenomena. 


2.2. Cavity Details 


A simplified drawing of the circular TMi:o mode 
cavity developed for making tensor permeability meas- 
urements at L-band frequencies is shown in figure 1. 
A resonant frequency of 1106 MHz was chosen which 
results in a cavity having an inside diameter of about 
13 in. and an inside height along the axis of approxi- 
mately 1.4 in. As previously mentioned, this readily 
allows the cavity to be placed between the pole pieces 
of a typical electromagnet since the magnetic field is 
applied parallel to the axis of the sample and cavity. 


TMyo CAVITY FOR TENSOR 
PERMEABILITY MEASUREMENTS 


rt 
seals 


COUPLING LOOPS 


| 


=— TUNING PROBE 


FIGURE 1. Simplified drawing of TMj,;9 mode cavity for tensor 


permeability measurements at 1106 MHz. 





The cavity itself was fabricated from a solid piece of 
brass with an access lid having a radius of 48 percent 
of the inside radius of the cavity which corresponds 
to a low current region. The solid piece of brass was 
employed to avoid the use of mechanical joints in the 
regions of high current such as at the periphery of the 
cavity. Earlier designs of cylindrical cavities using 
lids that made contact at the corners demonstrated 
that it is difficult to obtain a perfectly symmetrical con- 
tact resistance around the periphery of the cavity. 
Tuning probes are located in the region of high electric 
field on the solid end of the cavity as seen in figure 1. 
It was observed that better circular polarization condi- 
tions could be obtained by using the above techniques. 


Further minor adjustments are still necessary for 
obtaining optimum circular polarization conditions in 
the cavity. For example, as shown in figure 1, a heavy 
ring clamp is placed around the outer circumference 
of the cavity with adjusting screws for stressing the 
cavity in compression. This allows slight adjustments 
to be made in the circular symmetry of the cavity. 
Four stnall tuning probes as shown in figure 1 were 
placed in line with the coupling loops and at the point 
of maximum electric field to provide further fine 
adjustments in the tuning of the cavity. Adjustable 
coaxial coupling loops were also used. These of 
course readily allow coaxial lines, rather than wave- 
guide, to be used external to the cavity which is advan- 
tageous at these lower frequencies. At the same time, 
the coaxial loops easily allow small adjustments in 
the coupling to be made in order to obtain a pure circu- 
larly polarized mode. 

The procedure for obtaining circularly polarized 
waves in the cavity is carried out by first adjusting the 
compression screws, tuning probes, and coupling loops 
of the empty cavity in an iterating procedure to equal- 
ize the magnitude and frequency of the two possible 
modes of the degenerate mode cavity. Further adjust- 
ments are then made by placing a ferrite rod in the 
cavity and observing any resulting mode splitting 
when the ferrite is biased with a d-c field. Final 
adjustments using the coupling loops and line stretch- 
ers, which are external to the cavity, are then made in 
an effort to completely suppress one mode when the 
cavity is excited by a circularly polarized wave. It 
should be noted that it is more difficult to obtain the 
final adjustments with the sample in the cavity at 
L-band frequencies than at the higher frequencies. 
This is partly because the greater low field losses make 
it more difficult to distinctly see the splitting of the 
two modes. However, as will be shown in this paper, 
it is still possible to obtain good data. 


Figure 2 is a photograph of the cavity showing many 
of the features described above. A typical sample 
is shown at the center of the cavity. In an actual 
measurement, the sample would be inserted through 
the small hole in the lid and cavity as shown. The 
lid would of course be on the cavity and would not be 
opened during the measurement. The small rod 
protruding from the stressing ring in the lower right- 
hand corner of the photo is used for hanging the cavity 
between the pole pieces of the electromagnet. 
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FIGURE 2. Photograph of TM,19 mode cavity with lid removed show- 


ing typical sample. 


2.3. Measurement System 


The basic requirements in the circuit for tensor 
permeability measurements are shown in figure 3. 
The equipment used in this system is derived from 
commercially available coaxial components modified 
where necessary to meet the needs of these measure- 
ments. In particular, the microwave source is a 
commercial oscillator whose power supply was re- 


placed by well-regulated direct current power supplies 


to improve the quality and stability of the signal. The 
plate supply was regulated to 0.002 percent, and the 
filament supply regulation was adequate at 0.5 percent. 
The system also makes use of IF detection to increase 
sensitivity and a phase-locked local oscillator system 
was used to assist in tracking the microwave source 
during the measurements [31]. 

Optimum accuracy in the resonant frequency and Q 
measurements of the cavity was assured by using a 
transfer oscillator in conjunction with a frequency 
counter. The crystal detector and galvanometer was 
in turn calibrated using a 30-MHz piston attenuator in 
the IF circuit in order to improve the accuracy of the 
3-dB points on the cavity resonance curve used in the 
Q determinations. 

The equal 90° phase shifted signals needed for cir- 
cular polarization are obtained by the phase and power 
splitter in conjunction with the line stretchers. These 
adjustable length coaxial lines readily allow fine adjust- 
ments to be made on the phase which is necessary for 
obtaining good circular polarization. 


3. Measurement Results 


The permeabilities of three commercially available 
polycrystalline garnet materials of different magnetiza- 
tions were measured at 1106 MHz. These materials 
will be identified as follows: 



































FIGURE 3. Schematic diagram of circuit components used in tensor 
permeability measurements at L-band frequencies. 


Material Composition Saturation magnetization, 
477M in gauss 


Substituted YIG 580 
Substituted YIG 1158 
YIG 1799 








Measurements were made on rods of diameters 
0.100, 0.150, 0.200, 0.250, and 0.300 in. Although data 
could be obtained on all these rods, the smaller diame- 
ters (0.1 in—0.15 in) resulted in poor precision because 
of the resulting small frequency shifts while the large 
rods (0.3 in) proved difficult to measure because of the 
greater loss introduced into the system. The best 
data were obtained on the 0.250 in rods. 

The real and imaginary parts of the intrinsic tensor 
permeability values, (u’+a’) and (uw”+a”"), for the 
above three materials using 0.250 in diameter rods 
whose lengths ranged from 2.850 in to 2.882 in are 
shown in the top graph of figures 4,5, and 6. It is of 
interest that these relatively good tensor permeability 
resonance curves were obtained on rod samples whose 
ordinary Kittel resonance properties could not be 
measured at L-band frequencies using a rod geometry 
due to the low fields resulting in unsaturated conditions 
in the sample. 

It is feasible that low field effects may also be influ- APPLIED FIELD-OERSTEDS 
encing the accuracy of these tensor permeability _ 
measurements. For example, we have shown the FIGURE 4. Intrinsic and external tensor permeabilities for poly- 


“mn . Ath 4 : crystalline substituted yttrium iron garnet. material A M = 
external tensor permeability for positive circular polari- 580 G. garnet, material A, 4M 
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zation conditions corresponding tothe above intrinsic 
tensor results in the bottom graphs of figures 4, 5, and 
6. These curves were calculated from the intrinsic 
tensor values using eq (3). The external permeability 
curves are essentially just the high field side of the 
Kittel resonance curve. Based on data obtained at 
higher frequencies [5], both y'+.«’ and 7n”’+k" might 
be expected to be smoothly decreasing curves. How- 
ever, in all three cases, a slight absorption was ob- 
served in the 7” + x” data in the vicinity of the intrinsic 
tensor permeability resonance. This same phenomena 
was similarly observed in measuring one of the rods 
in a rectangular TE,o2 rectangular cavity at 1107 MHz 
indicating that it probably is not related to the adjust- 
ment of the circular polarization in our cylindrical 
cavity. It was also observed on rods ground into 
prolate ellipsoids indicating that it is not a macroscopic 
geometry effect. It is probably related to the rela- 
tively low fields used in these experiments although 
further studies are necessary in this area. 

It is also of interest that the small absorption located 
at a field of about 200 Oe in the w”—a” curve for 
material A in figure 4 appeared in measurements of 
different size rods of the same material and also in 
a prolate ellipsoid formed from one of the rods. This 
phenomena again may be characteristic of this material 
at low fields. 

Further studies were carried out in an effort to deter- 
mine the precision and accuracy of the system for 
measuring the intrinsic tensor permeability which 
was of primary interest in our program. First of all, 
it was observed that, for any given rod, the field for 
maximum absorption in w”+ a” could be repeated to 
better than | percent while the linewidth and maximum 
value of w”+a” at resonance could be repeated to 
better than 5 percent. However, it was noted that 
data obtained on different rods may show a much 
greater variation than this which is probably due to 
sample inhomogeneity effects, etc. 

The linewidths of the three materials obtained from 
the plotted results in figures, 4,5, and 6 were compared 
with corresponding measurements on disk samples 
in a rectangular TE;o2 mode cavity at 1107 MHz using 
the technique adopted by the ASTM [32]. Disks 
having aspect ratios (diameter/thickness) of the order 
of 83 to 85 were used in an effort to minimize shape 
effects [33]. The results are shown in the following 
table: 





FiGuRE 5. Intrinsic and external tensor permeabilities for polycrys- 
talline substituted yttrium iron garnet, material B, 4nM= 1158 CG. 








Linewidth (oersteds) 





Sample 
Tensor meas. | ASTM meas. 
(rods) (disks) 





APPLIED FIELD-OERSTEDS = 
51 48 


82 78 
FIGURE 6. Intrinsic and external tensor permeabilities for polycrys- _ 66 62 
talline yttrium iron garnet, material C, 4nM= 1799 G. 














The agreement is remarkably good considering the 
variations that are often experienced even from sam- 
ples cut from the same slab of material. These results 
indicate that the tensor results give a very good indica- 
tion of the linewidth in addition to providing information 
on the individual tensor permeability components. 

It is also recognized that the rods used at these 
frequencies do not approximate infinitely long cylin- 
ders, and demagnetizing corrections may be necessary. 
For example, it is well known that the internal field 
H; in a material is related to the external field Ho by 
the relation 

H;= H,— NM, (4) 
where JN is the demagnetizing factor and M is the mag- 
netization. The data shown in figure 7 for the three 
materials indicate that the applied field for resonance, 
H,, does increase slightly with increasing rod diameter. 
Better results are obtained using the corresponding 
data for the internal field H; calculated using eq (4). 
and N values obtained from Bozorth [34], with values 
of M corresponding to the internal field as determined 
independently with a magnetometer. It thus appears 
that small size effects in resonant field values can be 
expected; however, these can be partially corrected 
for in the case of the rod samples. 
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FIGURE 7. Dependence of intrinsic tensor permeability resonant 
field on sample diameter for materials A, B, and C. 
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No definite trend in linewidth on rod diameter 
could be identified within the variations of data ob- 
served on different rods. The linewidth varied con- 
siderably on different size rods that were cut from the 
same slab: therefore, it is possible that some size 
effects may still be present within the total variation. 

It is apparent that some of the above variations in 
resonant field and linewidth could be comparable to 
similar variations on size previously observed in ferro- 
magnetic resonance studies on disks [33]. It is also 
of interest to note that differences in the properties 
of ellipsoids and rods using the TMi cavity, as well 
as ellipsoids and disks measured in ,the previous pro- 
gram, were observed. It thus appears that rod samples 
in the TMi cavity as well as thin disk samples such 
as might be used in a stripline cavity may both result 
in some sample-shape errors in tensor permeability 
measurements. 

Although the 7M,,o cavity method does not com- 
pletely eliminate the problem of sample shape errors 
occurring in other methods, the cylindrical symmetry 
and ability to make several adjustments on the cavity 
for obtaining circular polarization may have advantages 
over stripline methods using disks. Furthermore, 
the cylindrical cavity allows the use of a long rod 
sample which is a distinct advantage when data in the 
vicinity of the intrinsic tensor permeability resonance 
are desired since this resonance coincides with the 
ordinary ferromagnetic resonance in disks but not in 
rods. The complications introduced by this higher 
loss in the disks have been previously recognized at 
higher frequencies [6] and further substantiated by 
our own efforts to measure disks at these lower 
frequencies. 


4. Conclusions 


On the basis of the above study, it is concluded 
that a degenerate TMiio mode cavity is very satis- 
factory for measuring the tensor permeability of rod 
samples at frequencies as low as 1 GHz. It is pos- 
sible to obtain good linewidth measurements on cylin- 
drical samples which ordinarily could not be measured 
at these frequencies using standard rectangular wave- 
guide techniques normally used in studies of the Kittel 
resonance. 

Although the use of rods in this method offers some 
decided advantages over disks in stripline when 
measuring intrinsic tensor permeability near reso- 
nance, the TMiio cavity method is still influenced 
by some sample shape effects. It is possible that all 
low frequency tensor permeability measurements 
may be influenced by sample shape effects that are 
related to the low fields involved. These results along 
with the appearance of a small anomalous absorption 
in the external permeability observed in our rod data 
in the vicinity of the intrinsic tensor permeability 
resonance indicate that there is still a need for further 
work describing ferromagnetic resonance and tensor 
permeability properties at low frequencies and applied 


fields. 





The authors express their appreciation to H. E. 
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Standard x-ray diffraction powder patterns, H. E. Swanson, 
M. C. Morris, and E. H. Evans, Mono. 25, Section 4 (June 28, 1966), 
55 cents. 

Standard x-ray diffraction powder data are presented for 103 sub- 
stances. Thirty-two of these patterns represent experimental data 
and 71 are calculated. Ten experimental patterns replace eleven 
cards already in the X-ray Powder Diffraction File published by 
the American Society for Testing and Materials; twenty-two experi- 
mental patterns and seventy-one calculated patterns are for sub- 
stances not previously included in the File. The experimental 
x-ray powder diffraction patterns were made with a Geiger counter 
x-ray diffractometer, using samples of high purity. All d-values 
were assigned Miller indices determined by comparison with the- 
oretical interplanar spacings and from space group extinctions. 
The densities and lattice constants were calculated, and the refrac- 
tive indices were measured whenever possible. The calculated 
x-ray powder diffraction patterns were obtained from single crystal 
structure data. The reported peak height intensities for calculated 
patterns were converted from integrated intensities. 

Accurate cell determination measurements for the internal standards 
were obtained by using a flat-plate back reflection focusing camera. 


Key Words: Constants, crystal, diffraction, lattice, measurements, 
powder, standard, structure, x-ray. 


Electrical parameters of precision, coaxial, air-dielectric 
transmission lines, R. E. Nelson and M. R. Coryell, Mono. 96 
(June 30, 1966), $1.25. 

Since precision coaxial connectors have become commercially 
available, precision coaxial, air-dielectric transmission lines are 
being widely used as radio frequency immittance standards. The 
evaluation of seven electrical parameters for four different line sizes, 
which are commonly used as standards, are presented in graphic 
form. The seven parameters, evaluated as functions of frequency, 
are inductance per inch, resistance per inch, characteristic im- 
pedance magnitude, characteristic impedance phase angle, attenua- 
tion per inch, phase-shift per inch, and wavelength. Also included 
are graphs showing how these parameters vary with changes in the 
resistivity of the conductors, the relative dielectric constant of air, 
and the diameters of the conductors. 


Key Words: Air-dielectric, coaxial transmission lines, electrical 
parameters, evaluation, functions of frequency, graphic presenta- 
tion, and precision. 


Dimensional metrology subject-classified with abstracts, 
I. H. Fullmer, Misc. Publ. 265 (August 1, 1966), $3.00. 

This bibliography covers two of the three principal divisions of di- 
mensional metrology, namely (1) linear, angular, and geometrical 
measurements of solid bodies, and (2) in-process control of sizes. 
The remaining division is gages, gaging, and inspection as to speci- 
fied limits of size. There are about 3,600 references classified under 
main subject headings in 12 sections, which are further subdivided 
into 102 subsections covering 100 subjects. The 12 main subjects 
are (1) general metrology, physical and dimensional; (2) length and 
diameter measurements by interferometry; (3) length and diameter 
measurements by mechanical methods; (4) calibration of line stand- 
ards of length, including tapes: (5) dimensional measuring instru- 
ment design features; (6) angle measurement: (7) measurements of 
deviations from geometrical regularity; (8) measurement of profile; 
(9) measurement of screw threads; (10) measurement of gears; (11) 
measurement of thickness of thin films: and (12) measurement and 
production techniques for accurate in-process control of size and 
form. The references in each subsection are arranged chronologi- 
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cally and a brief abstract of each is given with few exceptions. The 
bibliography is an attempt to organize and make more readily avail- 
able the extensive existing knowledge in the field covered. 


Measures for progress. A History of the National Bureau 
of Standards, R. C. Cochrane, Misc. Publ. 275 (1966), $5.25. 
The story of the first half century of the Federal Government’s 
primary scientific research body is told in appealing, interest-holding 
literary style. Starts with the founding of the Bureau in 1901, when 
this country was beginning its industrial growth, and _ physical 
weights and measures were in a hodge-podge of disorganization. 
Tells how the Bureau brought order out of this chaos and, through 
scientific research, went on to establish exacting, reliable new stand- 
ards involving incredible sophistication in physical measurement. 
Reflected throughout the book is the remarkable technological, 
industrial, and commerical growth that the country itself was ex- 
periencing, 1901-1951, and the key role played by the National 
Bureau of Standards in this evolving technology. Thorough docu- 
mentation and indexing, plus 15 Appendices containing a fund of 
information not heretofore assembled, make this a valuable refer- 
ence volume. 


Analytical mass spectrometry section: instrumentation and 
procedures for isotopic analysis, ed. W. R. Shields, Tech. Note 
277 (July 25, 1966), 60 cents. 

This report describes the general instrumentation of the Analytical 
Mass Spectrometry Section and the specific analytical techniques 
which have been devised for the measurement of isotopic ratios of 
Ag, Br, Cl, Cr, Cs, Cu, Mg, Pu, and U. Interim procedures for B, Li, 
Rb, and Sr are also given 

In the appendix some general statistical principles used in the 
design and analysis are briefly discussed; an example is given in 
detail illustrating the various steps involved leading from original 
data to the reported uncertainties for the isotopic ratio of bromine. 


Key Words: Instrumentation, isotopic analysis, mass spectrometry, 
procedures. 


Emission spectrometry, M. Margoshes and B. F. Scribner, 
Anal. Chem. 38, No. 5, 297R-310R (Apr. 1966). 

The literature of analytical interest in atomic emission spectrometry, 
flame photometry, and atomic absorption spectrometry is reviewed 
for the years 1964-65. The topics included are books and reviews, 
spectral descriptions and classifications, instrumentation, excita- 
tion, standards, calibration, and analytical applications. Literature 
citations total 351. 


Application of sulfur coatings to integrating spheres, S. T. 
Dunn, Appl. Opt. Letters 4, No. 7, 877 (July 1965). 

This letter describes several methods of preparing infrared inte- 
grating sphere coatings. 


Simple contact type temperature controller, L. A. Marzetta, 
Rev. Sci. Instr. Note 37, No. 6, 789 (June 1966). 

A simple, low cost temperature controller is described. The circuit 
has two independent current paths: One provides an On-Off control 
current, and the other a preselected current for the background 
heat. With the stated components, 250 watts of heater power can 
be controlled. For applications such as laboratory liquid baths, 
1/1000 °C stability is possible. 


Key Words: Mercury-toulene thermometer, temperature controller. 


Some statistical properties of pulses from photomultipliers, 
M. Gadsden, Appl. Opt. 4, No. 11, 1446-1452 (Nov. 1965). 

An examination of the integrated counts from-9502S and 9558A 
photomultipliers suggests that the statistical distribution is closely 





Poissonian. An examination of the time intervals between consecu- 
tive pulses has shown an interesting effect: with a 9502S tube, a 
highly significant, although small, number of intervals are much 
longer than expected on a random hypothesis. This effect is ex- 
plicable in terms of a loss of electrons between cathode and anode. 
Data on the variation of the signal : noise ratio with differing dis- 
criminator bias levels are also presented. 


Gas density stabilizer for flow proportional counters, R. D. 
Deslattes, B. G. Simson, and R. E. La Villa, Rev. Sci. Instr. 37, No. 
5, 596-599 (May 1966). 

A controller for gas density (pycnostat) is described with particular 
view to gain stabilization in flow proportional counters. The device 
is capable of operation at small flow rates and at pressures both above 
and below atmospheric. Constructional details and performance 
measurements are reported. 


Use of an auxiliary sphere with a spectrophotometer to ob- 
tain absolute reflectance, D. G. Goebel, B. P. Caldwell, and H. K. 
Hammond, J. Opt. Soc. Am. 56, No. 6, 783-788 (June 1966). 
Reflectance measurements that are made on a scale that is not 
relative to an arbitrary standard are often called “absolute” meas- 
urements. The method presented here uses an auxiliary sphere 
with a double-beam integrating-sphere spectrophotometer to make 
measurements on an absolute basis. The basic requirements are: 
(1) The auxiliary sphere must be uniformly coated with a highly- 
reflecting, highly-diffusing material. (2) A flat plate must be coated 
in an identical manner to provide a measure of the coating reflec- 
tance. (3) The interior-surface area of the sphere and the area of 
the entrance port must be measured. 

The theory of the method is discussed and an error analysis is made. 
Reflectance data are reported for specimens of smoked MgO and 
pressed powders of MgO and BaSQs. 

The precision of repeatability has been evaluated from measure- 
ments of a Vitrolite-reference standard. More than a dozen meas- 
urements at each of eight wavelengths made over a 3-year period 
exhibited a standard deviation of 0.003 for the spectral reflectance. 


Key Words: Absolute, reflectance, spectral. 


Influence of an external magnetic field on nuclear resonance 
in ferromagnetic materials, L. H. Bennett, J. Appl. Phys. 37, 
No. 3, 1242-1243 (Mar. 1, ! 966). 

The reduction in the amplitude of the nuclear magnetic resonance 
signal, as well as a change in line shape, upon the application of an 
external steady magnetic field has been measured for pure nickel 
at room temperature. Both the transverse and longitudinal con- 
figurations (steady field and rf field perpendicular and parallel, 
respectively) were measured for a sample of nickel powder and for 
bulk nickel. The signal amplitude drops faster in the longitudinal 
configuration than in the transverse for each specimen. The reduc- 
tion in amplitude is considerably less for the pellets than for the 
powder. There is no sign of a critical field as previously reported. 
A comparison of the resonance results with the apparent permeabil- 
ities is presented. 


Key Words: Critical field, field dependence, nickel pellets, nickel 
powder, nuclear resonance, permeability. 


Measurement of absolute wavelength stability of lasers, 
K. D. Mielenz, R. B. Stephens, K. E. Gillilland, and K. F. Nefflen, 
J. Opt. Soc. Am. 56, No. 2, 156-162 (Feb. 1966). 

A technique was devised to measure the absolute wavelength sta- 
bility of a gas laser by direct interferometric comparison with a mer- 
cury 198 standard lamp. The apparatus used is described; its limits 
of precision are discussed. 

The wavelength fluctuations of a free-running, unstabilized helium- 
neon laser were measured and found to be of the magnitude (several 
parts in 107) due to thermal and mechanical instabilities of the laser 
cavity. 

Manual control of the laser was seen to yield a wavelength constancy 
(a few parts in 108) comparable to the limits of accuracy of inter- 
ferometric standard sources. 


Key Words: Fringe, interferometer, laser, optical, photoelectric, 
precision, scanning, stability, standard, wavelength. 
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Measurement of the 633-nm wavelength of helium-neon 
lasers, K. D. Mielenz, R. B. Stephens, K. E. Gillilland, K. F. Nefflen, 
and R. B. Zipin, Appl. Phys. Letters 7, No. 10, 277-279 (Nov. 1965). 
The wavelength of the 3s.,—2p, transition of neon was measured 
by comparison of helium-neon lasers with the Krypton standard of 
length. 


Some applications of the wave front shearing interferometer, 
J. B. Saunders, Japanese J. Appl. Phys. Suppl. 1, 4, 99-105 (1965). 
This paper gives the results of several applications of the wave front 
shearing prism interferometer. The instrument is very compact 
and easy to apply. It is applied to the testing of chromatic aber- 
ration of simple and compound lenses; and to the testing of wave 
forms that characterize the monochromatic aberrations (spherical, 
coma and astigmatism). Results are shown for several different 
type lenses. This interferometer is equally applicable to the testing 
of small lenses and large telescope objectives. 


The measurement of homogeneity of optical materials in 
the visible and near infrared, F. W. Rosberry, Appl. Opt. 5, 
No. 6, 961-966 (June 1966). 

A procedure is described for determining the inhomogeneities of a 
sample of optical material in terms of small changes in index of re- 
fraction at discrete points over the measured surface. The method 
is used in the visible and near infrared regions of the spectrum. A 
He-Ne gas laser source was used for the infrared measurements. 
The variation in index was mapped with contour lines enclosing 
areas of similar index variations. The largest index change in one 
cm distance was noted and recorded as the maximum gradient per 
cm. Over 25 different samples of six different materials were ex- 
amined. The results are presented in a chart indicating the range 
of maximum index gradients of the samples observed. Data are 
shown to illustrate the possibility of determining index values to 
higher precision than the material justifies. 


Key Words: Homogeneity, index gradient, infrared, interferometer, 
laser, measurement, optical, refractive index, shadowgraph. 


Electron paramagnetic resonance and nuclear magnetic 
resonance as analytical tools, T. C. Farrar, (Proc. Conf. Puri- 
fication of Materials, New York Academy {nn. N.Y. 
Acad. Sci. 137, 323-334 (Jan. 1966). 

Nuclear magnetic resonance (NMR) spectroscopy is becoming an 
increasingly important and useful analytical tool. The advantages 
(such as rapid, non-destructive analysis), disadvantages (such as 
large sample requirement and low sensitivity), and limitations of 
high-resolution, broad-line and pulsed NMR techniques as applied 
to analytical problems are pointed out. Recent developments in 
instrumentation which have significantly increased the stability 
and sensitivity of NMR techniques are discussed. 

Examples of the use of high-resolution NMR to determine the struc- 
ture and the purity of samples, of the recently-synthesized molecules 
HBF:. SisBF., and SisBF. are given. The results of a broad-line 


of Sciences), 


NMR study to determine the number and the position(s) of the hydro- 
gen atoms in the molecule Re3(CO),.H,,are presented. 

Although paramagnetic species can be detected in concentrations as 
low as 10-'' M, electron paramagnetic resonance (e.p.r.) has been 


used for quantitative analysis only in a few cases. It is somewhat 
more useful for qualitative analysis. At this time its prospects for 
developing into a useful analytical tool seem limited. 
for this are briefly discussed. 


The reasons 


A cooling system for a laboratory magnet, R. A. Forman and 
T. McKneely, Appl. Spectry. 20, 189 (May-June 1966). 
The cooling system is described for use with a5 kW magnet system 


which enables one to obtain high electrical and thermal stability. 


A method for calibrating volt boxes, with analysis of volt- 
box self-heating characteristics, R. F. Dziuba and T. M. Souders, 
IEEE Intern. Conf. Record, Pt. 10, pp. 17 (Mar. 21-25, 1966). 

A modification of the Julie ratiometric method for calibrating volt 
boxes has been under study at NBS to evaluate its performance as 
a useful calibration technique. The modification features speed of 
operation, high accuracy and rated voltage measurements. Evalu- 
ation of the method accelerated an investigation of self-heating errors 





of volt boxes that was in progress. The results of this study to- 
gether with an analysis of self-heating curves are included. 


Key Words: Calibration, errors, heating, volt-box, voltage divider. 


Anomalous conductance behavior in polymers, A. H. Scott 
(1965 Annual Report Conf. Electrical Insulation, NAS—NRC Publ. 
1356, pp. 98-102 (Natl. Acad. Sci.-Natl. Res. Council, Washington, 
D.C., 1966). 

Conductance measurements on polymer materials are found to be 
affected by discharge currents even though the specimen has not 
had a voltage applied to it. It is necessary to fully discharge a 
specimen (sometimes for several days) before a repeatable conduct- 
ance measurement is obtained. The discharge currents are am- 
plified by the use of conductive rubber electrodes rather than tinfoil 
electrodes. Conductivities that were obtained with conductive 
rubber electrodes were lower in all cases than those obtained with 
tinfoil electrodes. 


Key Words: Anomalous d-c conductance, conductive rubber elec- 
trodes, epoxy, poly¢arbonate, 


polystyrene, 
vinylchloride, tinfoil electrodes. 


polyvinylbutyl, poly- 


Apparent negative impedances and their effect on three- 
terminal dielectric loss measurements, W. P.. Harris (1965 
{nnual Report Conf. Electrical Insulation), NAS-NRC Publ. 1356, 
pp. 83-86 (Natl. Acad. Sci.-Natl. Res. Council, Washington, D.C., 
1966). 
Tee networks consisting of “stray” impedances can cause “errors” 
in the conductance and capacitance of three-terminal capacitance 
measurements. In particular. a tee with two capacitive arms and 
a resistance leg to ground can introduce a negative component of 
conductance. Therefore, a three-terminal capacitor can appear to 
have zero or negative losses. A measurement of the dielectric loss 
of a material can indicate a negative loss unless steps are taken to 
eliminate by carefvl cell design this source of error. 
A component of negative loss can arise from the coaxial cables con- 
necting the cell or capacitor to the bridge. This error can be re- 
duced by the use of coaxial chokes’ or by a feedback circuit. 

7 
Key Words: Dielectric losses, dielectric materials, dielectric meas- 


urements, errors, low-loss, negative impedances, negative losses. 


Flow measurement standardization, T. Filban and M. R. Shafer, 
Proc. ISA Conf. 1964, 19, Pt. 1, Paper No. 12—2—4-64 (1964) 

The availability and types of standards related to the measurement 
of flow are reviewed. The facilities available at the National Bureau 
of Standards for the calibration of flow-measuring instruments are 
described: the inaccuracies or uncertainties associated with these 
measurements are discussed. 


Standards for galvanomagnetic devices, S. Rubin, Solid State 
Electronics 9, 5 559-566 (May 1966). 

A review is given of the Military Standard and the IEC standard for 
Hall effect devices. Points of outstanding interest such as linearity 
definitions are detailed. A review of a proposed standard for mag- 
netoresistive devices is given. Similarities between the two stand- 
ards are noted as are some of the similarities with standards for other 
flux sensitive resistors such as thermistors and photoresistors. 


The NBS instrumentation for use in the measurement of 
spectral irradiances of solar simulators, R. Stair and W. E. 
Schneider, Proc. Inter. Symp. Solar Radiation Simulation, Los 
Angeles, Calif., Jan. 18-20, 1965, pp. 75-94 (Inst. Environmental 
Sciences and Am. Soc. Testing and Materials, 1965). 

Three methods are described for the measurement of the spectral 
irradiances of solar simulators at the National Bureau of Standards. 
The design of acceptable equipment is significantly affected by large 
variations in sensitivity of photoelectric and thermoelectric detectors 
with the wavelength and over their surfaces. The first instrumenta- 
tion consists of a conventional spectroradiometer in which the source 
flux is collected in an integrating sphere which insures accurate 
evaluation of the solar simulator flux in terms of that from a 1000- 
watt lamp standard of spectral irradiance. A MgO-coated sphere is 
used for the ultraviolet, visible and near infrared to 1.6 microns. 
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For longer wavelengths a gold-coated sphere was found more effl- 
cient. The second instrumentation consists of a photoelectric 
filter spectroradiometer in which 36 interference filters mounted on a 
filter wheel are selected and placed in correct optical position 
through the use of a T—V rotor motor. In both of these instrumen- 
tations photoelectric detectors are employed with the flux beam 
chopped at 510 cps and the signal measured ac-wise. The third 
instrumentation is a thermoelectric filter spectroradiometer wherein 
the filters are likewise mounted on a filter wheel. The thermoelec- 
tric output is measured dce-wise with the detector calibrated through 
the use of a lamp standard of total irradiance. The supplementary 
use of the first two instrumentations is recommended. 


The response of flush diaphragm pressure transducers to 
thermal gradients, L. Horn, /S4 20th Annual Conf. and Exhibit, 
Los Angeles, Calif., Oct. 4-7, 1965, 20, Pt. III, ISA Preprint No. 
13.3-4-65 (1965). 

The zero signal output of pressure transducers may alter during 
field use in the presence of thermal gradients. The changes can be 
large enough to cause serious questioning of the validity of test 
data taken without corrections for the zero shift. 

Variations of up to 100% FS (full scale) have been found using 
transducers that indicated less than a 6% shift in the standard 
steady state temperature tests. 

A simple method of dynamic testing is described and experimental 
results are given for a number of tests. 


Key Words: Flush diaphragm, pressure gage, temperature, tempera- 
ture gradient, transducers. 


Color designation and specification, D. B. Judd, Encyclopedia 
of Industrial Analysis 1, 315 (1966). i 
Tristimulus colorimetry is summarized from basic assumptions to 
the 1931 CIE standard observer and the 1964 CIE suppjementary 
observer and coordinate system for large-field colorimetry. Three 
types of tristimulus colorimetry (visual, photoelectric, spectro- 
photometric) are described including operating procedures, methods 
of reducing the data, and a summary of applications. A similar 
treatment is given of approximate colorimetry by systematic collec- 
tions of color standards (Munsell, Ostwald, Maerz & Paul dictionary) 
including determination of the ISCC-NBS color designation. Defi- 
nitions of radiometric, photometric, and colorimetric terms used 
are given in the form of a glossary. 


Electrical testing in chemical technology, A. H. Scott, En- 
cyclopedia of Chemical Technology, 2d ed. 7, 716-726 (1965). 
Electrical testing of one sort or another is being used today to (a) de- 
termine the electrical properties of materials for design purposes, 
(b) determine the purity of the materials, (c) check the uniformity of 
manufacture or the manufacturing process, (d) aid in developing 
an understanding of the molecular structure of materials, (e) char- 
acterize materials, and (f) indicate the influence of environment on 
the properties of materials. The purpose of this article is to briefly 
describe such electrical measurements and to indicate their useful- 
ness and limitations. 


Measurement of DC dielectric conductance (reciprocal 
resistance) at elevated temperatures, A. H. Scott, Proc. 6th 
Electrical Insulation Conf., Sept. 13-16, 1965, pp. 252-254 (1965). 
The new term adopted by Committee D-9 ASTM (DC Dielectric 
Conductance) for the conductance observed in electrical insulating 
materials differentiates this type of conductance from that in metals 
or other conducting materials. Measurements at elevated tempera- 
tures require special precautions and special emphasis on some of 
the usual precautions of measurement. Special cells or holders 
are required for these measurements. 


The calibration of permanent magnet standards, |. L. Cooter, 
ISA 20th Annual Conf. and Exhibit, Part 1. Measurement Stand- 
ards and Instrumentation, Los Angeles, Calif., Oct. 4-7, 1965, 20 Pt. 
1, ISA Preprint 14.1—3—65 (1965). 

The increased demand for the calibration of permanent magnets 
used as reference standards has required a new magnet calibration 
facility at the National Bureau of Standards, Washington, D.C. The 
method procedure and apparatus used in the calibration of reference 
magnets are described. 





The nature of measurement, A. G. McNish, Book, Industrial 
“Metrology (Sponsored by the American Society of Tool and Manu- 
facturing Engineers), chs. 1, 4, No. 2, 24-31 (Mar. 1966). 

General principles for measurement. Accuracy and precision are 
discussed. 


Resolving power reiated to aberration, F. E. Washer, Photog- 
rammetric Eng. XXXII, No. 2, 213-226 (Mar. 1966). 

An heuristic method which accounts for the discrepancy between 
measured and theoretical values of resolving power of lenses is 
presented. The analysis also accounts for the sharp decline in 
contrast -with increasing values of resolving power expressed in 
lines per mm in the image plane. Analyses are given of the probable 
effects on resolution and contrast of the fraction of the lens area 
transmitting image-forming light, longitudinal chromatic aberration, 
and longitudinal spherical aberration. An example is given that 
demonstrates change in focus with change in resolving power and 
shows how it is readily explained on the basis of simple diffraction 
theory. 


Second breakdown and current distributions in transistors, 
H. A. Schafft and J. C. French, Solid-State Electron. 9, No. 7, 
681-688 (July 1966). 

A study of second breakdown in transistors, which included the use 
of temperature-sensitive phosphors to reveal current distributions, 
has emphasized how intimately the susceptibility of the transistor 
to second breakdown is linked to the internal distribution of current 
or of energy dissipation. The effect of such factors as the base 
drive and internal structural irregularities on the internal current 
distribution is described. A better understanding of the role of the 
base drive in second breakdown was achieved and led to a test 
which was used to distinguish second breakdown from other low 
voltage modes that have been recently reported to be different levels 
of second breakdown. 


Key Words: Diodes, failure mechanisms, internal current distribu- 
tions, second breakdown, semiconductors, temperature sensitive 
phosphors, transistors. 
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